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by Julie- Ann Goo diet 

1 define 3 samples of extragalactic radio sources of type FRll, containing 26 
objects in total. The control sample consists of 6C and 7C sources with radio- 
luminosities of around 10^''WHz~^ at 151MHz and redshifts oi z ~ 1. The other 
samples contain 3CRR sources with either comparable redshifts but radio lumi- 
nosities about a decade larger or with comparable radio-luminosities but redshifts 
around z ~ 0.4. These samples are used to investigate the possible evolution of 
depolarisation and rotation measure properties with redshift and radio luminosity 
independently. All sources show evidence for an external Faraday screen. Depo- 
larisation and variations in the rotation measure of a source increase with redshift, 
only the depolarisation measure shows any trend with radio-luminosity. There are 
no strong trends with spectral index. 

1 use the Kaiser et al. (1997) model with 3 fiux measurements to estimate the 
density of the source environment, source age and jet power. Depolarisation and 
variations in the rotation measure were found to be insensitive to changes in the 
density of the source environment. This indicates that depolarisation and variations 
in the rotation measure are not density indicators but map how the disorder in the 
magnetic field changes with redshift. 

The model is then used on a larger set of FRIl sources taken from the complete 
3CRR, 6CE and 7C HI samples. I find that the density of the source environment 
is a strong function of the source radio-luminosity but is insensitive to changes in 
redshift. The jet power is also a strong function of radio-luminosity but also a weak 
function of redshift. This indicates that the mass of the black-hole powering FRII 
sources is epoch dependent. In general, I find that the source environment is domi- 
nated by changes in the radio-luminosity but is relatively unaffected by changes in 
cosmic epoch. 
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Chapter 1 
Introduction 



The formation of all massive galaxies and also understanding their evolution is a vital 
component in understanding how the Universe evolved. Some of the most extreme, 
distant and powerful sources in the Universe have been observed in the radio band. 
AGN and radio galaxies were first solidly detected in radio in the 1950's using large 
aperture radio telescopes such as the Cambridge University and Sydney University 
radio telescopes (Krolik, 1999). The publication of the 3'''^ Cambridge survey (3C) 
in 1962 (Bennett, 1962) represented the first comprehensive list of radio sources 
in the northern hemisphere. Unfortunately it was not the first complete survey 
using the Cambridge telescope. Embarrassingly the 1*^* and 2"^^^ Cambridge surveys 
were found to be highly confusion limited and had to be abandoned (Mills & Slee, 
1957). The Cambridge telescope was later upgraded and the weU known 3C survey 
was produced. The revised 3C sample (3CR) is flux-limited, only sources brighter 
than 12 Jy at 178 MHz are included and hence only contains sources out to z ~ 2. 
Current technology has vastly improved on this limit, the new Texas - Oxford Survey 
(TOOT) (Hill & Rawlings, 2003) has a flux limit 100 times fainter than the 3CR 
sample and reaches as far back as 2; ~ 4. 

Radio sources are visible in all wavebands ranging from X-ray, through the optical 

and into the radio. The Cosmic Microwave Background radiation scatters off the 
radio synchrotron emitting electrons into the X-ray band and 3C 273 was the first 
radio source to be found in X-rays (Bowyer et al., 1970). The Compton scattered 
emission from radio jets emits in the X-ray but AGN also emit X-rays directly. It 
was not until the late 1970's that the Ariel V survey found a sub-class of radio 
galaxies that were strong X-ray emitters (Elvis et al., 1978). Radio sources often 
have optical counterparts that have been observed with the large optical telescopes 
such as Keck and the Hubble Space Telescope. By observing these sources in the 
optical, observers can determine the host galaxy properties. For many sources this 
is the only way to determine where the nucleus of the radio source should lie as it 
may be too faint to be present at radio wavelengths. Recent studies by Stevens et al. 
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(2003) using SCUBA in the sub-mm have found that AGN may be present in many 
if not all large elhptical galaxies. Radio galaxies have been found to contain some 
of the most massive and distant black holes observed (Willott et al., 2003; Almaini 
et al., 1999). 

Observing in the radio offers an unique opportunity that is absent from any other 
waveband. Radio observations are not limited by obscuration of dust such as in the 
UV and soft X-ray. Radio observations also do not depend on the time of day and 
only the most extreme of weather conditions can disrupt observations. Telescopes 
can be built near cities, e.g. Jodrell Bank, and do not need to be built on high 
mountain tops. Radio astronomy is often connected to the large impressive arrays 
of dish antennas that span many kilometres such as the Very Large Array (VLA) or 
the globe with Very Long Baseline Interferometry (VLBI). However, radio astronomy 
can and is being done on much smaller and simpler scales. LOFAR (low frequency 
array) for example is currently being built and will use 10,000 wire antennas, spread 
out across the Netherlands and Germany. 

It is perhaps the high redshift Universe that is best suited to radio observations 
and by studying how the environments of populations of radio galaxies evolve we 
can begin to understand how the Universe itself evolved and also understand how 
the properties of the host galaxy are affected by changes in the source environment. 



1.1 Classifying radio galaxies 

Fanaroff & Riley (1974) were the first to suggest that radio galaxies could be effec- 
tively split into two distinct classifications, the so called FRI's and FRII's. FRI's are 
bright near their centers and fade out towards the edge whereas FRII's are brightest 

at their edges and fainter towards the center. Around Pivsmhz ^ 10^^ Hz~^ there 
is a break between the two classifications, all FRI's, in general, lie below this break 
and all FRII's above. Owen & Ledlow (1994) found that the break in radio lumi- 
nosity extends into the optical and that FRIIs reside in the more optically-luminous 
galaxies. 

The distinction between the two classifications becomes even more complicated 
when their morphologies are considered: FRI's generally have both jet and counter- 
jet structures visible and exhibit complicated non-uniform structure (e.g. 3G 48 
Leahy et al., 2000). FRII's on the other-hand have very rarely been observed to 
have both jet and counter-jet structures. FRII jets are well confined and their 
lobe structure can be seen to extend up to several Mpc (e.g. 3C 236 and 3C 326). 
The luminosity break and the associated change in the source morphology has been 
suggested to be due to magnetohydrodynamical properties of the jet fiow which 
determine the stability of the jet (Rawhngs, 2002; Kaiser & Alexander, 1997). 
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Figure 1.1: High resolution image of Cygnus A by Carilli & Barthel (1996). This 
image was generated with data from telescopes of the National Radio Astronomy 
Observatory, a National Science Foundation Facility, managed by Associated Uni- 
versities. Inc. 



The most famous example of a FRII is 3C 405, otherwise known as Cygnus A 
(see Carilli & Barthel (1996) for a good overview on this source). Figure 1.1 shows 
a high resolution image of Cygnus A, a "Classical Double" source. The low redshift 
{z ~ 0.06, Owen et al., 1997) of this object means that its radio structure was been 
well studied. At the top-right and bottom-left of Figure 1.1 there are two large 
diffuse structures known as the radio lobes. Within each lobe is a bright locale, 
known as the radio hotspot. Cygnus A is atypical in the fact that it shows a strong 
core (near the center) with both jet and counter-jet structures propagating outward. 
Both lobes show extensive filamentary structure, see Fomalont et al. (1989) for a 
discussion on the nature of these filaments. Cygnus A is often said to be over- 
luminous (Carilli & Barthel, 1996) compared to most sources at this redshift. The 
relationship between redshift and radio-luminosity is discussed later in this chapter. 
In this thesis I concentrate on FRII objects only. 



1.2 Environment or Orientation? 



Many authors have used polarisation and flux measurements of powerful radio 
sources to argue that any observed asymmetries in these measurements are caused 
by orientation and/or environmental differences. Observing in different wavebands 
gives an insight into different regions of the source. In the optical, emission line prop- 
erties of a source can give an indication of the richness of the source environment. 
X-ray observations sample the hot gas around a source and can determine how this 
gas influences the radio source. In the following sections I review the evidence for 
the interaction of powerful radio sources with their environment and explain how 
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orientation effects also have a strong impact on what is observed at radio and other 
wavebands. I will then outline my own approach for determining the relationship 
between the observed source properties and the source environment. 



1.2.1 The Laing-Garrington effect 

Garrington et al. (1988) and Laing (1988) defined a sample of FRII sources with well 
defined jets. Both studies found that depolarisation, which is defined as the ratio of 

percentage polarisation at two separate frequencies, was almost always larger on the 
counter-jet side. This was deemed to be a simple orientation effect. The emission 
from the counter-jet, which is pointing away from an observer, has to travel along a 
longer path through the magnetised plasma in the vicinity of the source, causing it 
to appear more depolarised due to an inhomogeneous Faraday rotation (see chapters 
4 and 5). This is known as the Laing-Garrington effect. Gregory & Condon (1991) 
proved that depolarisation by gas in the lobes themselves was unlikely to cause 
the observed asymmetry. This is now the common assumption: the depolarising 
medium in radio galaxies is external to the lobes and any internal component will 
not significantly contribute to the overall depolarisation. 



1.2.2 Liu-Pooley effect 

Liu & Pooley (1991a,b) found that the least depolarised lobe in a source also has 
the flatter spectral index. General arguments for the asymmetry of the spectral 
index in the lobes can be attributed to Doppler beaming of the hotspots. However, 
Liu & Pooley (1991a,b) found that the entire lobe shows a flatter spectral index, 
and thus the observed spectral asymmetry must be due to some intrinsic difference 
of the source environments. Later studies by Ishwar a- Chandra et al. (2001) and 
Dennett-Thorpe et al. (1999) found that the Liu-Pooley effect was preferentially 
stronger for smaller sources and was independent of whether the sample comprised 
of radio galaxies, quasars or a mixture of both, see chapter 4 for a more detailed 
discussion. 

1.2.3 Evidence for a non uniform medium around radio 
galaxies 

Many authors have tried to determine if the so called Faraday screen around powerful 
radio sources is uniform or clumpy. The Faraday screen can cause the intrinsic 
polarisation vector of a source to rotate. It also causes a drop in the observed 
percentage polarisation if the structure of the magnetic fleld, within the Faraday 
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screen, varies on scales smaller than the size of the beam of the telescope used in 
the observations. Both effects are highly wavelength dependent. The angle over 
which the polarisation vector rotates is directly related to the line-of-sight strength 
of the magnetic field and the column density of electrons surrounding a source. 
Carilli & Barthel (1996) found that Cygnus A showed a magnetic field structure 
that was uniform over ~ 10 kpc. 

Pedelty et al. (1989) found that the lobe length directly related to the strength of 
the depolarisation. Almost all of their high redshift, high radio-luminosity sources 
showed a preference for the smaller lobe to be more depolarised. This was explained 
as an environmental difference between the two lobes of an individual source and 
Pedelty et al. (1989) concluded that the Faraday medium around their radio galaxies 
was distinctly clumpy. If one lobe was embedded in a denser environment, it would 
not be able to expand as freely as the other lobe and the denser environment would 
cause a larger degree of depolarisation. This trend with size has also been found 
in many other studies (e.g. Strom, 1973; Strom & Jagers, 1988; Best et al., 1997). 
Strom (1973) also concluded that the depolarisation-size trend could be explained by 
the presence of a large gaseous halo around a source. Small sources would be deeply 
embedded in this Faraday screen and would be expected to show more depolarisation 
than a source that was only partially embedded. This is consistent with observations 
that there is a distinct lack of linear polarisation at long wavelengths (A > 49cm), 
which can be explained as a result of a large halo of gas around the source, generally 
associated with the host galaxy (Strom & Jagers, 1988). This halo then creates a 
large Faraday rotation, vastly reducing the degree of polarisation observed. 

The extent to which the Faraday screen affects the observed radio properties of 
a source was first investigated by Burn (1966). More recent and perhaps realistic 
interpretations have been suggested by Tribble (1991) and EnBlin & Vogt (2003). 



1.2.4 X-ray 

A radio emitting source is generally assumed to be surrounded by a halo of hot, 
diffuse X-ray emitting material. Gregory & Condon (1991) found that the X-ray 
emission from elliptical galaxies was comprised of a symmetric halo, independent of 
whether the galaxy was in a group or a cluster. By combining X-ray observations 
with radio observations it is possible to show that in at least some sources, the 
diffuse halo of material is associated with the Faraday medium causing significant 
depolarisation (Gregory & Condon, 1991; Carilli & Barthel, 1996). 

Observations of radio sources using ROSAT and CHANDRA have shown that 
FRII sources have a strong effect on the material surrounding them. A good example 
of this interaction is Cygnus A where the radio lobes arc seen to forcibly expand 
into the X-ray emitting material (Carilli & Barthel, 1996). Young et al. (2001) also 
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found this effect in their CHANDRA observations of M87. X-ray observations can 
be used to directly determine the gas density in the source environment. In a recent 
study by Hardcastle & Worrall (2000), it was found that low redshift FRII's do 
not reside in dense environments which is consistent with the results using galaxy 
correlation methods. 



1.2.5 Galaxy correlation methods 



Galaxy-correlation methods are an indirect method of sampling the environment 
around powerful radio galaxies. We expect that the density of the gas in-between 
galaxies is higher in 'rich' environments, e.g. galaxy clusters, compared to 'poorer' 
environments, e.g. isolated galaxies. The extent to which galaxies cluster around 
each other can be analysed using cross-correlation methods (Longair & Seldner, 
1979; Prestage & Peacock, 1988; Yates et al., 1989; Yee & Green, 1987). This 
method is independent of any other observational method for sampling the source 
environment and only depends on how accurate galaxies can be correlated with each 
other. 

Longair & Seldner (1979) were among the first to try and determine the extent 
to which radio sources in the SCR and 4C samples (with z < 0.1), associated with 
the Zwicky galaxies (Zwicky et al., 1968), are clustered. Longair & Seldner (1979) 
found that extended powerful radio galaxies were in richer regions of space than 
galaxies selected at random. However, they also found that the sources with the 
strongest emission line spectra were isolated and very massive. This indicated that 
some sources could be massive enough to provide their own atmospheres and did 
not need to be within clusters or groups of galaxies. A later study by Prestage & 
Peacock (1988) using a much larger sample of galaxies out to 2; < 0.25, did not find 
this tendency for strong emission galaxies to be isolated, but did find that the most 
powerful galaxies were in rich clusters. 

The conflicting results lead to the question as to whether the difference in the 
richness of the environment that surrounded a source was due to the differences 
in the redshift or radio-luminosity of the samples used. To try to overcome this 
degeneracy in redshift and radio-luminosity Yates et al. (1989) studied yet a larger 
volume of space with 0.15 < z < 0.82. Yates et al. (1989) found that it was probably 
the differing radio-luminosity of the sources indicating differences in the richness of 
the environment, but they still could not rule out the possibility that there was some 
residual cosmic epoch effect. 
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1.2.6 Alignment effect 

Powerful radio sources are often found to be located in giant elliptical galaxies and 
also to contain extended regions producing strong emission lines (e.g Inskip et al., 
2002; Best et al., 2000a,b; Tadhunter et al, 1998; McCarthy et al., 1991; Ridgway 
et al., 2004; Inskip et al., 2003). These regions producing optical emission lines 
have often been found to ahgn closely to the radio emission. A striking case is the 
aligned radio galaxy 3C 280 (Ridgway et al., 2004), which displays a tight correlation 
between its rest frame UV line emission and its radio emission. An expansion of 
a radio source into the surrounding medium can trigger star-formation (Storrie- 
Lombardi et al., 1996; Dey et al., 1997). 4C 41.17 at z ~ 3.8 is perhaps the strongest 
case for triggered star-formation. Optical measurements of the stellar population 
show a striking correlation with the position of the radio lobes. Chambers et al. 
(1987) also found that higher redshift sources displayed a stronger alignment effect, 
with very little alignment present in lower redshift sources (Inskip et al., 2002). 
Interestingly Inskip et al. (2003) found that the strength of the alignment effect 
was also dependent on the radio-luminosity of the sample, with the strength of the 
correlation weakening in the fainter 6C sample compared to the 3CR sample. 

The observed alignment effect has also been found to depend on size, (Best et al., 
1996) and the wavelength at which the observations were taken (Inskip et al., 2003). 
Inskip et al. (2003) found that the alignment effect was non-existent in their infrared 
observations of the 6C and 3CR sources at z ~ 1. McCarthy et al. (1991) found 
that there was more emission-line gas present in the near side of a source compared 
to the far side which argues for an orientation dependence of the emission-line gas. 



1.2.7 Spectral Ageing 

The spectral index of a source has been found, in general, to be flatter in the vicinity 
of the hotspot and steeper nearer the core regions (e.g Alexander & Leahy, 1987). 
This can be explained through energy losses of the electrons due to synchrotron 
radiation. The electrons towards the core are older than those closer to the hotspot 
and so have undergone more spectral ageing. Spectral ageing thus offers a unique 
way to determine the ages of radio sources. There are three main assumptions in 
this method: 

i. ) The main site of particle acceleration is the hotspot at the end of the jet, 

ii. ) The magnetic field is constant through-out the source, 

iii. ) The distribution of electrons can be described as a power-law. 
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Blundell & Rawlings (2001) argue against spectral ageing in radio sources, suggesting 
that the observed trend could be caused by a gradient in the magnetic field combined 
with a curved power-law distribution. Spectral ageing also assumes that there is no 
mixing between different populations of particles and that the radiative lifetime of 
synchrotron electrons is longer than the spectral age. Blundell & Rawlings (2001) 
suggest that both assumptions may be unrealistic. 

Although the spectral ageing method may be too simplistic it is interesting to 
note that Alexander & Leahy (1987) find that the environmental effects on a source 
only become important around the FRI/FRII break. In other words, more luminous 
sources do not show spectral signatures caused by environmental effects. 



1.2.8 Unified Model 

There are many sub-classes within the FRII definition. Perhaps the most well known 
are the terms, radio galaxy, quasar and blazar. Originally these sources were thought 
to be separate and independent populations. It was proposed (e.g Barthel, 1989a) 
that the difference between these three populations was simply the angle at which 
the source was viewed, see Figure 1.2. At large angles to the line-of-sight an observer 
would simply see a radio galaxy with narrow line emission properties. As the angle 
decreases, the flux ratio of the lobes increases and the emission becomes dominated 
by the lobe pointing towards the observer. The narrow-line emission can no longer 
be seen and only broad line regions are visible, which is known as a broad line quasar. 
At very small angles to the line-of-sight the emission is completely swamped by the 
highly polarised jet emission. The separate FRII sub-classes can be simply stated a 
function of orientation angle (Scheuer, 1987; Peacock, 1987; Barthel, 1989a). 




blazar 



radio galaxy 

radio galaxy 

Figure 1.2: Simple cartoon picture of the unifled model. The line-of-sight orientation 
determines if a source is observed to be a radio galaxy, quasar or blazar. 

If the orientation angle was a factor in whether of not a source was observed to 
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be a quasar or a 'normal' radio galaxy would expect all quasars and blazars to show 
evidence of beaming. Barthel (1989b) did indeed find that all the quasars in his 
sample showed some evidence of beaming. 



1.2.9 Source or environment? 

Sections 1.2.1 to 1.2.8 indicate that it is very hard to disentangle whether or not the 
gaseous environments of FRII radio sources change systematically as a function of 
redshift and/or radio luminosity. This is further complicated by orientation effects 
(e.g. Garrington et al., 1988), Doppler boosting (e.g. Leahy & Perley, 1991; Liu & 
Pooley, 1991a), size (Best et al., 1996) and wavelength (Inskip et al., 2003). 

Studies are further complicated by the use of fiux-limited samples and the result- 
ing degeneracy between redshift and radio-luminosity. 



1.3 Flux limited samples 

A purely observational effect hampering the study of source environments is the 
Malmquist bias inherent in any fiux-limited sample. Figure 1.3 demonstrates how 
severe this problem is in the case of the 3CRR sample by Laing et al. (1983). All 
complete surveys in the radio are flux-limited thus any analysis of the environments 
of these sources with respect to cosmic epoch or radio-luminosity will be greatly 
limited. To overcome the degeneracy between redshift and radio-luminosity it is 
possible to use sources taken from several samples, generally with different flux- 
limits, e.g. combining the 3CRR, 6CE and 7CRS samples (Blundell et al., 1999). 

Figure 1.3 shows that there is a distinct lack of high radio- luminosity sources at 
lower redshifts. In fact, all low redshift sources (albeit a few over-luminous sources 
like Cygnus A) are at much lower luminosities than their high redshift cousins. 
As we observe at higher redshifts we are observing larger volumes of space. At 
low redshift the chance of finding a relatively powerful source is small because the 
radio luminosity function is steep. As the sampled volume of space increases, the 
probability of finding a rare high-luminosity source naturally increases. One of the 
aims of radio astronomy is to 'fill in' the z-P plane and remove the strong bias that 
is present in all surveys, even when using several complete samples, at different flux 
hmits. Hill & Rawlings (2003) gives a concise review of the flux-hmit problem. 



1. Introduction 



-10- 



X >k 5<^'^ 



X X 



Xx ^ 







0.5 



1 .5 



z 



Figure 1.3: A radio luminosity-redshift plot for the 3CRR sample. 



1.4 Modelling FRII's and their environment 



In the previous sections 1 have reviewed observational techniques to test the envi- 
ronments of powerful radio sources. However, we can also attempt to model these 
sources and try to gain insights into their environments from a more theoretical 
point of view. 

Scheuer (1974) and Blandford & Rees (1974) were the first to explain the double- 
lobed structure. In the simplest model, energy is transported from the nucleus of 
the source to the radio lobes by a relativistic beam, more commonly known as the 
radio-jet. Scheuer (1974) predicted that only a small fraction of the energy in the 
jet would be radiated away from the tip of the beam, known as the hotspot. The jet 
terminates at the hotspot and creates a jet shock. The lobe is generally assumed to 
be over-pressured and expands supersonically into the external medium, creating a 
bow shock (Alexander, 2002). A simple picture of the dynamics of a radio source is 
presented in Figure 1.4. 

FRll's are generally assumed to expand with time along the jet axis and the 
pressure exerted by the jet is balanced by the ram pressure of the external gas 
(e.g Scheuer, 1974; Kaiser & Alexander, 1997; Kaiser et al., 1997; Hardcastle & 
Worrall, 2000). The jet must have enough ram pressure to overcome the external 
gas pressure, otherwise the source will not have a well defined hotspot (Prestage 
& Peacock, 1988). The rate of expansion of the lobe is highly dependent on the 
strength of the jet and the density of the source environment (Kaiser & Alexander, 
1997). If a source was strongly constrained perpendicular to the jet axis, then there 
would be a tendency for sources to be long and thin which has not been observed 
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IGM 



Figure 1.4: A simple impression of a FRII source. 



(Miller et al., 1985). 

In this thesis I apply the model of Kaiser Sz Alexander (1997) and Kaiser, Dennett- 
Thorpe & Alexander (1997) to the radio observations of FRII radio galaxies and 
quasars. Using this model I obtain estimates for the gas density in the environ- 
ment of these sources. I then compare these mo del- dependent estimates with my 
observational results on the depolarisation of the radio emission. 



1.5 This work 

It is often difficult to obtain information on the properties of the black holes at the 
centres of active galaxies, their jets and their large scale environments from direct 
observations. As I have demonstrated, many studies to try and disentangle the 
source/environment interaction suffer from the redshift-radio luminosity degeneracy 
present in all flux-limited samples. 

To break the degeneracy between redshift and radio luminosity I chose 3 subsam- 
ples of sources from the 3CRR and 6C/7C catalogues. The control sample consists 
of 6C and 7C sources at 2; ~ 1 with radio luminosities of around 10^^ WHz^^ at 
151 MHz. Another sample at the same redshift consists of 3CRR sources with radio 
luminosities around a magnitude higher at 151 MHz. The final sample consists of 
3CRR sources at 2; ~ 0.4, again with radio luminosities of around 10^'' WHz~^ at 
151 MHz. The observations of these sources can then be used to study the source 
properties and the medium around the source, thus discovering which correlate with 
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redshift and which correlate with radio luminosity enabling the following questions 
to be answered: 

i. ) Does a relationship exist between radio luminosity and the environment in 

which a given radio source lives? 

ii. ) Do the source environments evolve with redshift? 

To arrive at answers to these questions I investigate in detail the depolarisation 
properties of my sample. I also apply the models of source evolution mentioned 
above to my observational results. Finally, I apply the later technique to a much 
larger sample. 

As previously stated depolarisation and Faraday rotation are wavelength depen- 
dent. At low frequencies the degree of polarisation would be too low for an in-depth 
analysis of the environment and at high frequencies some of my sources would have 
too low signal-to-noise. Thus an intermediate frequency range of 1.4 GHz to 4.8 
GHz was chosen. 

The body of this work is dedicated to answering the question as to whether it is 
the redshift or the radio-luminosity of powerful radio sources that is the dominant 
factor in determining their observed characteristics. 



1.6 Synopsis 

In chapter 2 I give a detailed description of the sample selection and subsequent 
methods used in the data reduction. I also present individual notes for each source 
used in the sample. All sources which have pre-published maps between 1.4 GHz 
and 4.8 GHz are detailed in Table 2.6. 

In chapter 3 1 present the polarisation and flux properties for each source present 
in my sample and analyse the bulk trends present between my sub-samples. I 
also present the results of correcting the depolarisation measurements to a common 
redshift, z — 1, using the Burn (1966) law. 

In chapter 4 I present a more detailed analysis of the results presented in chapter 
3. I use the Spearman Rank, Partial Spearman Rank and Principal Component 

Analysis tests to determine how the rotation measure and polarisation parameters 
relate to the fundamental parameters of redshift, radio-luminosity and radio size. I 
also analyse the asymmetries present in the flux and polarisation results to determine 
how the environments of these 26 sources differ and if these differences can be related 
to the fundamental parameters. 
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In chapter 5 I present the results of using the Tribble (1991) models to correct 
the depolarisation measurements to a common redshift, -2=1, and the results are 
compared with those from chapter 3 using the simpler Burn (1966) prescription. 
I then use the Tribble (1991) models to determine if there is any evolution in the 
Faraday screen with respect to redshift or radio-luminosity. 

In chapter 6 I use the Kaiser & Alexander (1997) and Kaiser et al. (1997) models 
to estimate the density of the source environment and lobe pressure for each source 
described in chapter 2. I compare the density of the source environment and the 
lobe pressure, from the model, with the observed depolarisation and variations in 
the rotation measure to determine if the polarisation parameters are affected by 
the source environment. I then apply the models to the FRII sources from the 
combined 3CRR, 6CE and 7C III samples to estimate the density of the source 
environment, lobe pressure, jet-power and lobe age for each source. This allows a 
global approach to determine how the environments of powerful radio sources evolve 
with cosmological epoch and also with radio-luminosity. 

Chapter 7 details the conclusions drawn from my findings and gives possible avenues 
for an extension to this study. 

In appendices A to C I present the radio maps of all 26 sources from the small 
sample. Each source has an associated polarisation intensity and polarisation angle 
map at 4.8 GHz and 1.4 GHz, a spectral index and depolarisation map between 1.4 
GHz and 4.8 GHz and a rotation measure and magnetic field angle map between 1.4 
GHz and 4.8 GHz. 3C 457 is the exception as this source does not have any rotation 
measure or magnetic field maps associated with it. All chapters use a cosmology 
which assumes Ho = 75 kms~^Mpc~^, and Qm = 0.35 (A = 0.65). For full details of 
the cosmology used see appendix D. Finally, appendices E to G detail the archive 
information used in the modeUing of the 3CRR, 6CE and 7C HI samples. 



Chapter 2 



Observations and Data Analysis 



Observations of the polarisation properties of extragalactic radio sources can provide 
information on the relationships between the radio source properties and their en- 
vironments as well as the evolution of both with redshift. Many previous studies of 
variations in polarisation properties have suffered from a degeneracy between radio 
luminosity and redshift due to a Malmquist bias, present in all flux-limited samples. 
A good example of this effect is the depolarisation correlations (i.e. the ratio of per- 
centage polarised flux at two frequencies) found independently by Kronberg et al. 
(1972) and Morris & Tabara (1973). 

Kronberg et al. (1972) found that depolarisation of the radio lobes generally 
increased with redshift whereas Morris & Tabara (1973) found depolarisation to 
increase with radio luminosity. Due to the flux-limited samples (PKS and 3C) used 
by both authors it is difficult to distinguish which is the fundamental correlation, or 
whether some combination of the two occurs. Both suggestions have ready explana- 
tions: (i) If radio sources arc confined by a dense medium then synchrotron losses 
due to adiabatic expansion are reduced, the internal magnetic field is stronger and 
a more luminous radio source results; if this confining medium also acts as a Fara- 
day medium, more luminous sources will tend to be more depolarised, (ii) Sources 
at different cosmological epochs may reside in different environments and/or their 
intrinsic properties may change with redshift. 

Hill & Lilly (1991) observed that galaxy densities around FRII radio sources 
increased with redshift out to 2; ~ 0.5 and beyond, but Wold et al. (2001) did not 
find this trend in a recent study. Welter et al. (1984) argued that the increase in 
rotation measure with redshift is primarily attributable to an increasing contribution 
of intervening matter. However, depolarisation asymmetries within a source, e.g the 
Laing-Garrington effect, increase with redshift which imply an origin local to the 
host galaxy (Garrington & Conway, 1991; Laing, 1980; Laing & Peacock, 1980). 
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2.1 Observational program 
2.1.1 Sample selection 
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Figure 2.1: A radio luminosity-redshift plot showing the 3 subsamples used in the 
observations. Sample A is represented by 'x', sample B by 'o' and sample C by '*'. 
The lines mark the flux limits for the 3CRR and 6C samples. A spectral index of 
-0.75 was used to shift the 3CRR data, originally observed at 178 MHz to 151 MHz. 

Sample A was deflned as a subsample chosen from the 6CE (Eales et al., 1997) 
subregion of the 6C survey (Hales et al., 1990), and the 7C HI subsample (Lacy 
et al., 1999), drawn from the 7C and 8C surveys (Pooley et al., 1998). The selected 
sources have redshifts 0.8 < z < 1.3, and radio luminosities at 151 MHz between 
6.5 X lO^^WHz-i < PisiMHz < 1-35 X lO^^WHz-^. Sample B was deflned as a 
subsample from the revised 3CR survey by Laing et al. (1983) containing sources 
within the same redshift range, but with luminosities in the range 6.5 x 10^^ WHz~^ < 
-Pi5iMHz < 1-35 X lO^^WHz"^. Sample C is also from the 3CRR catalogue; it has 
the same radio luminosity distribution as the control sample, sample A, but with 
0.3 < z < 0.5. I only include sources that are more luminous than the flux limits 
of the original samples at 151 MHz (Figure 2.1). In all samples only sources with 
angular sizes Q > 10" (corresponding to ~ 90 kpc at z=l) were included, see Table 
2.1. This angular size limit is imposed by the depolarisation measurements. I 
required a minimum of ten independent telescope beams (1" per beam) over the 
entire source to ensure sufficient signal-to-noise. The distributions of linear sizes of 
the radio lobes are reasonably matched across all the samples (Figure 2.2), but note 
the two 'giant' sources in sample B. 

Each sample initially contained 9 sources. The source 3C 109 was subsequently 
excluded from sample C as the VLA data for that source was of much poorer quality 
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Figure 2.2: Linear size-redshift plots of the 3 subsamples used in the observations. 
Symbols as in Figure 2.1. Figure (a) assumes Ho = 75 kms~^Mpc~^, and = 0.5, 
fl\ = 0. Figure (b) assumes Ho = 75 kms~^Mpc~^, and Qm = 0.35, = 0.65. 



than that for the rest of the sample. The sources in the 3 subsamples are represen- 
tative of sources with similar redshifts, radio-luminosities and sizes. However, the 
samples are not statistically complete because of observing time limitations. Sources 
that were fairly well observed at 4.8 GHz at B array (and C array if needed) were 
selected from the archives, this ensured that only a minimum of new observations 
were needed. 

The ratio of angular to physical size varies only by a factor 1.3 between z = 0.4 
and z = 1.4. This ensured that all the sources were observed at similar physical 
resolutions. All redshift values for 6C sources were taken from Bales et al. (1997) 
except for 6C 1018+37 which was taken from Rawlings et al. (2001). Redshifts for 
the 7C sources were taken from Lacy et al. (1999). Redshifts for the 3CRR sources 
were taken from Spinrad et al. (1985), 4C 16.49 was taken from Barkhouse & Hall 
(2001), 4C 14.27 was taken from Herbig & Readhead (1992) and 3C 457 was taken 
from Hewitt & Burbidge (1991). 



2.1.2 Very Large Array observations 



Observations of all 26 radio galaxies were made close to 1.4 GHz using the A-array 
configuration and a 25-MHz bandwidth. This bandwidth was used instead of 50 
MHz to reduce the effect of bandwidth depolarisation which occurs when there are 
large changes in the rotation measure on the same scale as the beam size. The 
maximum angular size that can be successfully imaged using A array at 1.4 GHz is 
38". 12 sources were larger than this and were observed additionally with B array. 
Observations were also made at 4.8 GHz using a 50 MHz bandwidth. The maximum 
observable angular size in B array at 4.8 GHz is 36". The same 12 sources as before, 
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Source z Pi5imhz Angular size (T4.8GHz (Jiaguz 

(W/Hz) (arcsec) {/j, Jy) (/i Jy) 



6C nQ43+3q 


1 


04 


1. 


n 


X 


10"^ 

J- \J 


10 


20 


55 


an im 1-1-36 


1 


04 


1 


1 


X 


10^7 


49 


15 


65 


an 1018-1-37 





81 


8 


8 


X 


1Q26 


64 


24 


50 


6C 1129+37 


1 


.06 


1, 


,1 


X 




15 


21 


60 


6C 1256+36 


1 


.07 


1, 


,3 


X 




14 


18 


70 


6C 1257+36 


1 


.00 


1, 


,1 


X 


10^^ 


38 


33 


100 


7C 1745+642 


1, 


.23 


9. 


,5 


X 


1026 


16 


30 


60 


7C 1801+690 


1 


.27 


1. 


,0 


X 


10^7 


21 


24 


75 


7C 1813+684 


1 


.03 


7. 


,1 


X 


1026 


52 


23 


47 


3C 65 


1 


1 8 


1 


n 


X 


iq28 


17 


22 


61 


68 1 

L/C . J. 


1 


.24 


1 


1 


X 


1Q28 


52 


41 


100 

X \J\J 


3C 252 


1 


.11 


7, 


,2 


X 


10^7 


60 


23 


100 


3C 265 





.81 


6. 


,5 


X 


10^7 


78 


32 


70 


3C 268 1 





.97 


1. 


,0 


X 


1028 


46 


32 


68 


3C 267 


1 


.14 


1. 


,0 


X 


1028 


38 


25 


61 


3C 280 


1 


.00 


1. 


,2 


X 


1028 


15 


33 


75 


3C 324 


1 


.21 


1. 


,3 


X 


1028 


10 


26 


70 


4C I6.49 


1 


.29 


9. 


,9 


X 


10^7 


16 


25 


81 


3C 16 





.41 


7, 


,4 


X 


1026 


63 


36 


80 


3C 42 





.40 


7. 


,5 


X 


1026 


28 


29 


70 


3C 46 





.44 


7. 


,9 


X 


1027 


168 


20 


53 


3C 299 





.37 


6. 


,9 


X 


1026 


11 


39 


11 


3C 341 





.45 


8. 


,8 


X 


1026 


70 


26 


88 


3C 351 





.37 


7, 


,6 


X 


1026 


65 


31 


100 


3C 457 





.43 


9. 


,6 


X 


1027 


190 


18 


80 


4C 14.27 





.39 


8. 


,8 


X 


1026 


30 


30 


90 



Table 2.1: Details of the sources in sample A, B & C. Sources in italics are quasars, 
gives the Icr noise level in the final total flux maps at frequency v. 

were then observed at 4.8 GHz with C array. This ensured that both 1.4 GHz and 
4.8 GHz observations were equally matched in sensitivity and resolution. Details of 
the observations are given in Tables 2.2 to 2.4. 

Sources in the 6C and 7C samples have a typical bridge surface brightness of 
~ 70// Jy beam~^ in 4.8 GHz A-array observations (Best et al., 1999). In order 

to detect 10% polarisation at 3(7 in B-array observations I required an rms noise 
level of 20/i Jy bcam^^, corresponding to 70 mins of integration time. At 1.4 GHz, 
assuming a = —1.3, bridges will be a factor of 4 more luminous. At this frequency, 
the integration time is set by the requirement to have an adequate amount of uv- 
coverage to map the bridge structures. 20 minute observations were split into 4x5 



2. Observations 



-18- 



minute intervals. This observation splitting to improve uw-coverage was also done 
for the 4.8 GHz data. 

As Table 2.2 demonstrates, for many of my sources the integration times at 
4.8 GHz are considerably less than the 70 min requirement, due to telescope time 
constraints. At 1.4 GHz the integration time on all the sources is above the minimum 
required for good signal-to-noise. This is not the case at 4.8 GHz. Many of the 
observed properties that depend on polarisation observations (e.g depolarisation and 
rotation measure), are therefore poorly measured in the fainter components at 4.8 
GHz. The values obtained are then only representative of the small region detected 
and not the entire component. Spectral index is independent of the polarisation 
measurements and so it is relatively unaffected by the short integration times. 

The 3GRR sources are more luminous but much of this is due to the increase in 
the luminosity of their hotspots; their bridge structures are only a few times brighter 
than those of the 6C/7C HI sources. To reach a Sa detection of 7% polarisation on 
the bridge structures, a total integration time of 30 mins was required at 4.8 GHz 
and 20 mins at 1.4 GHz, split into 3-4 minute intervals to improve ttv-coverage. The 
vast majority of the sources in sample B and C had at least this minimum amount 
of time on source at each frequency (see Tables 2.3 and 2.4). 



2.1.3 Observational Programs 

AD429 Most observations at 1.4 GHz in A array were obtained on 31/07/99 with 
the observing program AD429. The data from this day is strongly affected by a 
thunderstorm at the telescope site during most of the observations. Even after 
removal of bad baselines and antennas the noise level in this data remained at least 
twice the theoretical value. However, careful calibration and CLEANing reduced 
this effect to a minimum. Sample A was most affected by the thunderstorm and the 
lack of observing time at all frequencies. However, I find that the results obtained 
by Best et al. (1999) for some of the sources in sample A are in good agreement with 
my results. I am confident that my data is reliable for fluxes above the Sa^^ise level. 
The polarisation calibration of AD429 was compared to the B-array data at 1.4 
GHz (for the 12 sources that had B-array data). This confirmed that the position 
angle (PA) of the polarisation vector in both data sets agreed to within 15 degrees in 
all sources. This additional check was used to ensure that the 1.4 GHz polarisation 
angle calibration was accurate. AD429 also contains observations of sources in C 
array at 4.8 GHz. As these observations were taken on 12/06/00, these observations 
are obviously not affected by the thunderstorm that affects the 1.4 GHz data. 



AD444 Tables 2.2 to 2.4 show that AD444 was used to observe some of the sources 
in B array at 4.8 GHz and also to complement the AD429 1.4 GHz observations in 
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Source 


Array 


Frequency 


Bandwidth 


Observing program 


Int. 




configuration 


(MHz) 


(MHz) 


(dd/mm/yy) 


(min) 


6C 0943+39 


A 


1465,1665 


25 


31/07/99 (AD429) 


16 




B 


4885,4535 


50 


20/05/01 (AD444) 


31 


6C 1011+36 


A 


1465,1665 


25 


31/07/99 (AD429) 


16 




B 


1452,1652 


25 


20/05/01 (AD444) 


17 




B 


4885,4535 


50 


25/02/97 (AL397) 


21 






4885,4535 


50 


20/05/01 (AD444) 


30 




C 


4885,4535 


50 


12/06/00 (AD429) 


20 


6C 1018+37 


A 


1465,1665 


25 


31/07/99 (AD429) 


16 




B 


1452,1652 


25 


20/05/01 (AD444) 


17 




B 


4885,4535 


50 


20/05/01 (AD444) 


31 




C 


4885,4535 


50 


12/06/00(AD429) 


20 


6C 1129+37 


A 


1465,1665 


25 


31/07/99 (AD429) 


16 




B 


4885,4535 


50 


20/05/01 (AD444) 


17 






4885,4535 


50 


25/02/97 (AL397) 


21 


6C 1256+36 


A 


1465,1665 


25 


31/07/99 (AD429) 


16 




B 


4885,4535 


50 


27/02/93 (AR287) 


15 


6C 1257+36 


A 


1465,1665 


25 


31/07/99 (AD429) 


16 




B 


4885,4535 


50 


20/05/01 (AD444) 


16 






4885,4535 


50 


25/02/97 (AL397) 


22 


7C 1745+642 


A 


1465,1665 


25 


31/07/99 (AD429) 


16 




B 


4885,4535 


50 


20/05/01 (AD444) 


11 






4885,4535 


50 


23/11/97 (AL401) 


31 


7C 1801+690 


A 


1465,1665 


25 


31/07/99 (AD429) 


16 




B 


4885,4535 


50 


26/03/96 (AB978) 


29 






4885,4535 


50 


23/11/97 (AL401) 


17 


7C 1813+684 


A 


1465,1665 


25 


31/07/99 (AD429) 


16 




B 


1452,1652 


25 


20/05/01 (AD444) 


16 




B 


4885,4535 


50 


20/05/01 (AD444) 


39 






4885,4535 


50 


23/11/97 (AL401) 


19 




C 


4885,4535 


50 


12/06/00 (AD429) 


20 



Table 2.2: Details of the VLA observations for sample A with the integration times 
included. See Tables 2.3 and 2.4 for samples B and C respectively. 
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Source 


Array 


Frequency 


Bandwidth 


Observing 


program 


Int. 




configuration 


(MHz) 


(MHz) 


(dd/mm/yy) 


(min 


3C 65 


A 


1465,1665 


25 


31/07/99 


(AD429) 


16 




B 


4885,4535 


50 


20/05/01 


(AD444) 


20 


3C 68.1 


A 


1417,1652 


25 


31/07/99 


(AD429) 


16 




B 


1417,1652 


25 


13/07/86 


(AL113) 


20 




B 


4885,4535 


50 


19/07/86 


(AB369) 


300 




C 


4885,4535 


50 


12/06/00 


(AD429) 


20 


3C 252 


A 


1465,1665 


25 


31/07/99 


(AD429) 


16 




B 


1465,1665 


25 


20/05/01 


(AD444) 


27 




B 


4885,4535 


50 


19/07/86 


(AB369) 


97 




C 


4885,4535 


50 


12/06/00 


(AD429) 


20 


3C 265 


A 


1417,1652 


25 


31/07/99 


(AD429) 


16 




B 


1417,1652 


25 


13/07/86 


(AL113) 


30 




B 


4873,4823 


50 


17/12/83 


(AM224) 


238 




C 


4873,4823 


50 


12/06/00 


(AD429) 


20 


3C 267 


A 


1465,1665 


25 


31/07/99 


(AD429) 


16 




B 


4873,4823 


50 


17/12/83 


(AM224) 


56 


3C 268.1 


A 


1417,1652 


25 


31/07/99 


(AD429) 


16 




B 


1417,1652 


25 


13/07/86 


(AL113) 


30 




B 


4885,4835 


50 


15/08/88 


(AR166) 


20 






4885,4835 


50 


01/06/85 


(AR123) 


21 




C 


4885,4835 


50 


06/11/86 


(AL124) 


102 


3C 280 


A 


1465,1665 


25 


31/07/99 


(AD429) 


16 




B 


4873,4823 


50 


17/12/83 


(AM224) 


46 


3C 324 


A 


1465,1665 


25 


31/07/99 


(AD429) 


16 




B 


4873,4823 


50 


17/12/83 


(AM224) 


51 


4C 16.49 


A 


1465,1652 


25 


31/07/99 


(AD429) 


16 




B 


4885,4535 


50 


04/03/97 


(AB796) 


30 



Table 2.3: Details of the VLA observations for sample B 



2.1. Observational program 



-21- 



Source 


Array 


Frequency 


Bandwidth 


Observing prog. 


Int. 




coiifiguratioii 


(MHz) 


(MHz) 


((l(l/iriiri/A'A') 


(iiiin) 


3C 16 


A 


1452,1502 


25 


14/09/87 (AL146) 


59 




B 


1452,1502 


25 


25/11/87 (AL146) 


39 




B 


4885,4535 


50 


20/05/01 (AD444) 


20 






4885,4535 


50 


17/11/87 (AH271) 


10 




c 


4885,4535 


50 


12/06/00 (AD429) 


20 


3C 42 


A 


1452,1502 


25 


14/09/87 (AL146) 


40 




B 


4885,4535 


50 


23/12/91 (AF213) 


67 


3C 46 


A 


1452,1502 


25 


31/07/99 (AD429) 


16 




B 


1452,1502 


25 


25/11/87 (AL146) 


35 




B 


4885,4535 


50 


20/05/01 (AD444) 


20 




c 


4885,4535 


50 


12/06/00 (AD429) 


20 


3C 299 


A 


1452,1502 


25 


31/07/99 (AD429) 


16 




B 


4835,4535 


50 


20/05/01 (AD444) 


15 






4885,4835 


50 


28/01/98 (AP331) 


15 


3C 341 


A 


1452,1502 


25 


14/09/87 (AL146) 


38 




B 


1452,1502 


25 


25/11/87 (AL146) 


47 




B 


4885,4535 


50 


20/05/01 (AD444) 


11 






4935,4535 


50 


26/10/92 (AA133) 


25 




C 


4885,4535 


50 


12/06/00 (AD429) 


20 


3C 351 


A 


1452,1502 


25 


31/07/99 (AD429) 


16 




B 


1452,1502 


25 


25/11/87 (AL146) 


56 




B 


4885,4835 


50 


19/07/86 (AB369) 


12 






4885,4535 


50 


20/5/01 (AD444) 


16 




C 


4885,4835 


50 


09/10/87 (AA64) 


22 


3C 457 


A 


1452,1502 


25 


31/07/99 (AD429) 


16 




B 


1452,1502 


25 


25/11/87 (AL146) 


30 




B 


4885,4535 


50 


20/05/01 (AD444) 


50 




C 


4885,4535 


50 


12/06/00 (AD429) 


20 


4C 14.27 


A 


1452,1502 


25 


14/09/87 (AL146) 


28 




B 


4885,4535 


50 


20/05/01 (AD444) 


17 



Table 2.4: Details of the VLA observations for sample C 
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B array. The observations were all taken on the 20/05/01 and were not affected by 
any serious problems. 



Archival Data For most sources some archival data has been used to supplement 
AD429 and AD444. Tables 2.2 to 2.4 detail which archival programs were used and 
the dates on which the sources were observed. 



2.2 Data Reduction 

The data reduction of the VLA data was done using the National Radio Astronomy 
Observatory's Astronomical Imaging Processing System or AIPS. All data from the 
archives was re-reduced to maintain consistency of data reduction methods when 
analysing the results. 



2.2.1 Flux calibration. 

The first step in the data reduction process involves removing the first 10s (usually) 
of each observation in the data set. This process is used to remove any residual er- 
rors from antenna tracking problems after moving to the target source. Once these 
initial scans are removed then the data is 'flagged' to remove any further errors in 
the data. Data is usually flagged for many reasons, such as interference from ter- 
restrial sources (e.g mobile phones), errors with the antenna and receivers and bad 
weather. If a thunderstorm is present during an observing run, then this will cause 
large errors in the data due to the highly ionised clouds. This was unfortunately the 
case for most of the data set AD429 at 1.4 GHz and results in a higher noise level 
for the sources, see section 2.1.3. All scans flagged are entered into a 'FG' table 
within AIPS and are no longer used in the rest of the data reduction process. 

A good example of where flagging is essential can be seen in the source 7C 
1745-1-642 as observed as part of AD429. Figure 2.3(a) shows the effect of the 
thunderstorm on 7C 1745-1-642. The usual flux of this source is ~ 80 mJy, but in 
several timesteps the flux jumps to more than 700 Jy, completely swamping the 
weak source underneath. By careful flagging of the bad timesteps 1 was able to 
retrieve the flux from the source which can be seen in Figure 2.3(b). Although this 
is an extreme case, it does demonstrate the need for flagging. After all the corrupted 
data has been removed the observations are ready to be calibrated. 
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(a) 



(bi 



20 40 60 80 100 120 140 160 180 
Kilo Wavlngth 



60 SO 100 120 140 
Kilo Wavlngth 



Figure 2.3: (a) A plot of the flux against u-v for the source 7C 1745+642. The radio 
flux of this source hes in the mJy range indicating that the burst exceeding 700 Jy 
is from a hghtning bolt near the array. The lightning bolt affected 20 seconds of 
the data, (b) Removal of the affected 20 seconds results in a realistic UV spectrum. 
Only every 10th visibility is plotted. 



2.2.1.1 Calibrating the primary calibrator. 



The AlPS task SETJY uses a primary flux calibrator such as 3C 286 or 3C 48 to 
determine the absolute flux density scale on which all other calibrations are based. 
3C 286 and 3C 48 are observed every few years with the VLA in D-configuration. 
Their flux density are recalculated accordingly using the modifled Baars scale (Baars 
et al., 1977), 

logS ^ A + Blogiy + C(log + D(log uf (2.1) 

where S is in Jansky, u is in MHz and Table 2.5 gives the values of the constants 
A, B, C and D for 3C 286 and 3C 48 which I used here. SETJY uses the calculated 
flux densities of these primary calibrators to set the absolute flux calibration. 

Source A B C D 

3C 48 1.31572 -0.74090 -0.16708 0.01525 
3C 286 1.23734 -0.43276 -0.14223 0.00345 

Table 2.5: The constants for equation 2.1 with the 1999.2 values. 



The task VLACALIB was then used to determine the antenna based calibration 
for the primary calibrator. Ideally the primary calibrators should be point sources 
that do not have any time variability. However, 3C 286 is partially resolved with 
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most array combinations of the VLA and most wave bands. It is resolved on two 
scale lengths, it has a small secondary core located 2.5" from the central core and 
the core itself is partially resolved at longer baselines. To over come this it is possible 
to restrict the uv-range and (if necessary) the number of inner antennas used in the 
cahbration fit. 3C 48 is also resolved with some array configurations and thus it 
also has some restrictions on its use. It also shows some time variability on the time 
scale of years. However, both calibrators are still the preferred primary calibrators 
and the restrictions do not adversely affect the flux calibration. For full details see 
the AIPS cahbrator manual. 



2.2.1.2 Calibrating the secondary 

The primary calibrator sets the absolute flux density of the observations. Several 
secondary calibrators that lie as close to the target sources as possible are observed 
as frequently as the target sources. These secondary calibrators are essential in 
the phase calibration of the sources. The secondary calibrators are calibrated using 
the task GETJY which calculates the flux density of the secondaries with respect 
to the primary, thus 'bootstrapping' the initial calibration on to the secondaries. 
VLACLCAL then interpolates between individual observations of the secondary cal- 
ibrator, thus applying the calibration to the entire data set. The last step in the 
flux calibration is to self-calibrate all the calibrators, still using VLACLCAL. This 
ensures that the flux calibration is complete. A flnal check on the flux density cali- 
bration is done using the task LISTR with OPTYPE='MATRIX' and DPARM = 31 
0. This displays all the rms noise values for the calibrators. A high value of the rms 
noise indicates a bad calibration. In this case the calibration process was repeated 
until a satisfactory rms value was achieved. 



2.2.2 Polarisation calibration. 

Polarisation calibration is more straight forward than the flux density calibration. 
However, the range of polarisation angles (PA) must be known for the calibrators as 
a successful calibration requires a polarisation angle change of more than 90 degrees 
over the course of individual observing programmes. PCAL reads in the uv-da.ta. 
and determines the effective feed parameters for each antenna and also for each IF 
(intermediate frequency)*. All corrections are then applied to the 'AN' (antenna) 
table. PCAL removes the parallactic angle from the phase of the visibility data and 
applies the corrections to the instrumental corrections that are placed in 'AN'. 

*Thc IF is the frequency that comes from mixing the local oscillator (LO) frequency with the 
radio frequency (RF). The RF is the incoming frequency of the observations and the LO is a 
sinusoidal signal that converts the RF to an IF. 
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2.2.3 Calibration of source files 



The flux calibration and the polarisation calibration is now complete for the ob- 
servations as a whole, but the sources must be individually 'cleaned'. There are 
several calibration methods that can be used, the main two used are CLEAN and 
the Maximum Entropy Method (MEM). The CLEAN method by Hogbom (1974) 
represents a radio source as a collection of point sources in an empty field. An iter- 
ative process is used to find the flux and position of these point sources, removing 
10% from the fiux of the point sources, with the maximum fiux, in each iteration. It 
then selects the most probable image from a set of images. MEA-I is not a procedural 
method, it selects the image with the largest entropy that fits the data and lies 
within the noise level. For full details of both methods see Cornwell et al. (1999). I 
use CLEAN to calibrate all my sources, but in order to check the quality of my maps 
I also use MEM, implemented in the AlPS task VTESS, instead of the CLEAN al- 
gorithm. The resulting maps are not significantly different from those produced by 
the CLEAN algorithm. In fact, despite common opinion VTESS is not necessarily 
superior to CLEAN in producing accurate maps of extended low surface brightness 
regions (Rupen, 1997). 



2.2.3.1 Clean 



Before the maps of individual sources can be cleaned they are first spht off from the 
entire data set using the task SPLIT. The source is then loaded into IMAGR and 
the number of CLEAN iterations set via NITER. NITER determines the number of 
times IMAGR will subtract 10% from the current peak flux. By choosing too small a 
value of NITER the task IMAGR will not reach the noise level. By choosing too large 
an NITER the CLEAN algorithm will clean past the noise level, i.e. it will assume 
that the noise is part of the source structure and will try to incorporate it into a 
model of the source. This 'over-cleaning' can lead to errors in the final source map. 
Thus it is a sensitive balance and a conservative value for NITER is usually the best 
option for the first round of cleaning. When IMAGR is finished it produces a set 
of clean components (CC's). The CC's are essentially a representation of a map of 
the source including the noise. Selecting a set of CC's excluding the noise provides 
a template map for the observed source. This subset corresponds to a specific uv- 
range of the source and this range must be used in CALIB for self-calibration. The 
sub-set of the CC's and their corresponding uv-range are used in conjunction with 
SOLMODE='p' which means that the initial calibration is only concerned with the 
phase of the source. The template map of the source is used in self-calibration to 
reduce the residual noise, in other words, the template is used to modify the uv data 
to more closely fit the template. The new uv data set is then run through IMAGR 
again and the whole process is repeated. In general the self-calibration is used with 
SOLMODE='p' for two iterations and then SOLMODE='a&p' (i.e. both amplitude 
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and phase) for all remaining iterations. After the final self-cleaning step the uv data 
set is ran through IMAGR with NITER set to a much higher value to remove the last 
remnants of any remaining noise from the stokes T image, i.e. the total intensity 
image. IMAGR is then run on a reduced number of NITER with stokes 'Q' and then 
'U^'^ to produce the raw maps that can be combined to produce polarisation maps 
described below. 



2.3 Map production 

Total intensity maps were made from the Stokes I parameters at each frequency. 
Polarisation maps were also made at all frequencies by combining the Stokes Q and 
U polarisation parameters. A map was then produced that contained the polarised 
flux, P = (Q^ + f/^)^/^, and the electric field position angle, PA = 0.5 tan" -'^(^), 
at a given frequency. The AIPS task POLCO was used to correct for Ricean bias, 
which arises when the Stokes Q and U maps are combined without removing noise- 
dominated pixels. By careful setting of the PCUT parameter this bias is removed. 
All maps only contain pixels where the polarised flux and the total intensity flux 
are above Sunoise at 4.8 GHz and Sunoisc at 1.4 GHz. The lower threshold at 1.4 GHz 
was necessary because the 1.4 GHz data had a higher noise level, so blanking flux 
below 5a resulted in large regions of polarised flux being lost. 

At all frequencies the individual maps were made such that the beam size, the 
cell size of the image and the coordinates of the observations were exactly the same. 
If any of these properties of the map differed between frequencies, then the resultant 
multi-frequency map would contain false structures that would be directly related 
to the mis-alignment of the maps. To make sure that the coordinates (and cell 
size) were always within acceptable tolerances the AIPS task HGEOM was used to 
realign maps at one frequency to maps at another frequency. In sources where an 
identiflable core exists at both frequencies, the core positions were used as a check 
on the alignment from HGEOM. In general HGEOM is adequate in aligning the 
multi-frequency data. Sources with a distinct core at all frequencies were aligned 
within 0.03", where no core existed the hotspots were aligned within 0.045". In 4 
sources this was not sufficient. 3C 68.1 had to be shifted 0.05" east and 0.07" north, 
3C 265 had to be shifted 0.04" west and 0.02" north, 3C 299 had to be shifted 0.1" 
east and 0.03" north and flnally 3C 16 had to be shifted 0.1" east and 0.1" north. 
All shifts were applied to the 4.8 GHz observations. 

Spectral index maps were made between 4.8 GHz and 1.4 GHz, where the spectral 
index, a, is deflned by oc z/". Depolarisation maps were made by dividing the 
map of the percentage polarised flux at 4.8 GHz by the corresponding map at 1.4 

= {RL + LR)/2 and U = \{RL — LR)/2, where R is the right circular polarisation and L is 
the left circular polarisation (Fomalont & Perley, 1999). 
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GHz. Rotation measure maps were made using the polarisation angle maps at three 
frequencies, 

PA{X) = PA, + RMX^, (2.2) 

where PA is the observed polarisation position angle of a source, PA, is the initial 
polarisation angle before any Faraday rotation, RM is the rotation measure and 
A is the wavelength of the observations. Clearly, polarisation angle measurements 
are ambiguous by ±n7r and this can introduce ambiguities in the rotation measure 
maps. A change of tt between 1.4 GHz and 4.8 GHz introduced by the used fitting 
algorithm will cause a change of 80 rad m~^ in rotation measure. To determine 
if any strong rotation measure change is real within a rotation measure map a plot 
of the polarisation angle (measured on both sides of the observed jump) against A^, 
including nvr ambiguities can be produced. The best fit from AIPS is then overlayed. 
Any true feature will not show any ambiguities in titt. This has been done for two 
sources: 4C 16.49 and 6C 1256+36. The resulting fits are presented in Figures 2.4 to 
2.5 and their corresponding values for the fits arc discussed in the relevant notes 
on these sources below. Another test that a feature is real is that a true jump in 
rotation measure causes depolarisation near the jump, but the magnetic field map 
shows no corresponding jump in the region. 

The rotation measure maps of 7C 1813+684, 3C 65 and 3C 268.1 contained 
obvious jumps in position angle which I was not able to remove. Plots analogous 
to Figures 2.4 to 2.5 indicated that there were regions that obviously contained 
errors caused by utt ambiguities. As previously noted the A array AD429 data was 
problematic and this was found to be the cause of the jumps. To overcome this 
problem I artificially shifted the position angles at 1.4 GHz data down by 10 to 15 
degrees before the PA maps were produced. This shift is within the position angle 
error and resolved any ambiguities. 

Table 2.4 shows that all sources in sample G were observed with IFs separated 
by only 50 MHz or less at around 1.4 GHz. This means that they were not well 
enough separated at 1.4 GHz to overcome the mr ambiguities. To compensate for 
this lack of separation the 4.8 GHz observations were split into their two component 
frequencies, 4885 MHz and 4535 MHz. I then used 4 frequencies for the fit instead 
of 3, but the sources in sample G are still only marginally sensitive to m: jumps. 
The resulting rotation measure maps cover the same frequency range as samples A 
and B, but use different frequencies for the fit. This was not possible in the case of 
3G 351 and 3C 299, resulting in larger uncertainties in the rotation measurements 
for these sources. 

In the case of sample C any source that has a large range of rotation measures 
(> 80 rad m^^), the AIPS task RM will force the rotation measure into a range ±40 
rad m~^ around the mean rotation measure. This is due to the lack of frequency 
separation at 1.4 GHz and it can cause jumps. In the case of 3C 457 these jumps 
were severe and I was unable to resolve them. The rotation measure and magnetic 
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field maps for this source arc not included in the analysis. The rotation measure 
varies smoothly over all other sources in this sample. The error affects the absolute 
value of the rotation measure for each source and therefore it does not affect the 
difference in the rotation measure between the two lobes of a given source, dRM or 
the rms variation in the rotation measure. 

Depolarisation and rotation measure properties are discussed in more detail in 
chapter 3. 



2.3.1 Notes on individual sources 

In appendices A to C (Figures ?? to ??) I present the maps of the radio properties 
discussed above for all the sources from the 3 samples. Each figure shows the 
depolarisation map (where the polarisation was detected at all frequencies), the 
spectral index map, the rotation measiu^e map and the magnetic field direction map 
(when a rotation measure is detected). Table 2.6 contains a listing of previously 
published maps for all sources, for completeness I have included all my reduced 
maps in the appendices but these maps are not published in Goodlet et al.(2004). 
In all cases only regions from the top end of the grey-scales saturate, as I have always 
kept the lowest values well inside the grey-scales, to ensure that no information has 
been lost. 



2.3.1.1 Sample A: 

6C 0943+39: (Figure A.l) No core is detected in my observations. Best et al. 
(1999) detected a core at 8.2 GHz and minimally at 4.8 GHz. My non-detection 
is probably due to the different resolution of the two data sets. The value of the 
rotation measure in the Eastern lobe must be considered with caution as it is based 
on only a few pixels. 



6C 1011+36: (Figure A. 2) This is a classic double-lobed structure, showing a 
strong core at both 4.7 GHz and 1.4 GHz with an inverted spectrum. 



6C 1018+37: (Figure A. 3) The maps were made with the smaller arrays only at 
each frequency. In the 1.4 GHz A-array data set the lower lobe was partially resolved 
out, but this was compounded by the high noise. So no feasible combination of the 
A and B array was possible. To maintain consistency the B-array 4.7 GHz data was 
also excluded. 
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Source 


Map 


Frequency 
(GHz) 


Reference 


6C 0943+39 


P 


4.8 


Best et al. (1999) 


6C 1011+36 


P 


4.8 


Best ct al. (1999) 




TI 


1.4 


Law-Green et al. (1995) 


6C 1129+37 


P 


4.8 


Best et al. (1999) 




TI 


1.4 


Law-Green et al. (1995) 


6C 1256+36 


P 


4.8 


Best et al. (1999) 




TI 


1.4 


Law- Green ct al. (1995) 


6C 1257+36 


P 


4.8 


Best et al. (1999) 




TI 


1.4 


Law-Green et al. (1995) 


3C 65 


TI 


1.4, 4.8 


Polatidis et al. (1995) 


3C 68.1 


P 


4.8 


Bridle et al. (1994) 




TI 


1.4 


Leahy et al. (1989) 


3C 252 


P 


4.8 


Fernini et al. (1993) 


3C 265 


P 


4.8 


Fernini et al. (1993) 


3C 267 


TI 


4.8 


Best et al. (1997) 






1.4 


Leahy et al. (1989) 


3C 268.1 


P 


4.8 


Laing (1981) 




TI 


1.4 


Leahy et al. (1989) 


3C 280 


P 


1.4, 4.8 


Liu & Pooley (1991b) 




S, D 


1.4, 4.8 


Liu & Pooley (1991b) 


3C 324 


TI 


4.8 


Best et al. (1998) 




P 


1.4 


Fernini et al. (1993) 


4C 16.49 


P 


4.8 


Lonsdale ct al. (1993) 


3C 16 


TI 


4.8 


Giovannini et al. (1988) 




P 


1.4 


Leahy & Perley (1991) 


3C 42 


P 


4.8 


Fernini et al. (1997) 


3C 46 


TI 


4.8 


Giovannini et al. (1988) 






1.4 


Grcgorini ct al. (1988) 


3C 299 


P 


1.4, 4.8 


Liu k Pooley (1991b) 




S, D 


1.4, 4.8 


Liu & Pooley (1991a) 


3C 341 


P 


1.4 


Leahy & Perley (1991) 


3C 351 


TI 


4.8 


Bridle et al. (1994) 




P 


1.4 


Leahy & Perley (1991) 


3C 457 


P 


1.4 


Leahy & Perley (1991) 


4C 14.27 


P 


1.4 


Leahy & Perley (1991) 



Table 2.6: Details of previously published maps. TI = total intensity, P = polari- 
sation, S = spectral index & D = depolarisation 
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6C 1129+37: (Figure A. 4) The SE lobe contains two distinct hotspots, Best et al. 
(1999) found 3 hotspots. This discrepancy is probably due to the different reso- 
lutions of the two observations. The source shows distinct regions of very strong 
depolarisation, however these regions are slightly smaller than the beam size. 

6C 1256+36: (Figure A. 5) The rotation measure map shows distinct changes in 
the values of the rotation measure. Figure 2.4 shows that although the jump in RM, 
in the southern lobe, does not correspond to a jump in depolarisation, it is not due 
to any error in the fitting program. The corresponding reduced values for the 
fits are given in Table 2.7. 




2 / \ «.S 06£ 06.4 06:2 06.0 05.S 05.6 05.4 

RIGHT ASCENSION (J2000) 



Figure 2.4: (a) A plot of the polarisation angle against A^, allowing for n7r variations, 
for the southern lobe of 6C 1256+36 with the dashed hnes showing the best fit 
models. All reasonable utt solutions for the 4.7 GHz data are considered and plotted. 
Both sides of the jump are plotted with 'o' indicating one side of the jump and 'x' 
the other, (b) The rotation measure map for 6C 1256+36 (rad m~^) between 4710 
MHz, 1665 MHz and 1465 MHz. All contours are at 5(7 at 4710 MHz (0.25 mJy 
beam-^)x(-l, 1, 2, 4.. .,1024) with a beam size of 2.5" x 1.4". 



Source 


Lobe 


n = -1 


n—Q n—1 


Symbol 


6C 1256+36 


S 


382.0 


51.0 1.9 


X 






981.6 


1.2 976.5 


o 



Table 2.7: Reduced values for the rotation measure fits for 6C 1256+36. 

6C 1257+37: (Figure A. 6) A core was detected at 4.7 GHz but was absent from 
the 1.4 GHz data. The high noise level and short observation time meant that the 
S lobe had very little polarised flux above the noise level. Reliable values for the 
rotation measure were found in only a few pixels around the hotspots. 
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7C 1745+642: (Figure A. 7) This is a highly core dominated source, with the north- 
ern lobe appearing faintly. There is an indication of a jet-like structure leading down 
from the core into the southern, highly extended, off-axis, lobe. The source is a weak 
core dominated quasar (Barkhouse & Hall, 2001). 



7C 1813+ 684: (Figure A. 9) This is the faintest of the sources in sample A and 

is also a quasar (Barkhouse & Hall, 2001). The source shows a compact core that 
is present at all observing frequencies, but it is too faint to detect any reliable 
polarisation properties. 



2.3.1.2 Sample B: 

3C 65: (Figure B.l) The W lobe shows a strong depolarisation shadow that is 
smaller than the beam size. Best (2000) found the source to lie in a cluster which 
might account for the presence of the depolarisation shadow and the large depolar- 
isation overall. 



3C 68.1: (Figure B.2) The source is a quasar (Bridle et al., 1994). A core has been 
detected by Bridle et al. (1994) in deeper observations. 



3C 252: (Figure B.3) The SE lobe shows a sharp drop in the polarisation between 
the 4.7 GHz and 1.4 GHz observations. 



3C 265: (Figure B.4) The NW lobe shows evidence of a compact, bright region with 
a highly ordered magnetic field which at higher resolutions Fernini et al. (1993) show 
is the primary hotspot. 



3C 267: (Figure B.5) The E lobe is highly extended, reaching to the core position, 
which can be seen in the 1.4 GHz image. The large depolarisation region in the W 
lobe coincides with a region with no observed rotation measure. The core is strongly 
inverted with a = 0.48. 



3C268.1: (Figure B.6) This is a classic double-lobcd source but there is no core 
detected at either 1.4 GHz or 4.8 GHz. 
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3C 280: (Figure B.7) The value of the rotation measure and the magnetic field 
direction in the E lobe must be treated with caution as it is based on only a small 
region of the entire lobe. The sharp changes in the rotation measure map are not 
seen in the magnetic field map and the depolarisation map shows a similar structure 
suggesting that is not due a fitting error. 



3C 324: (Figure B.8) The NE lobe shows evidence of a depolarisation shadow. 
Best (2000) found the source to lie in a cluster which may explain the faint shadow. 



4C 16.49: (Figure B.9) The source is a quasar (Barkhouse & Hall, 2001), that 
shows a strong radio core, jet structure and possibly a small counter-jet. The source 
is highly asymmetric with the southern lobe almost appearing to connect to the core. 
It has a very steep spectral index, a < —1.0 making it an atypical source. Figures 
2.5(b) and 2.5(c) demonstrates that the sharp changes in the rotation measure map, 
see Figure 2.5(a), are not due to any fitting errors. The corresponding reduced 
values for the fits are given in Table 2.8. 



Source 


Lobe 


n = —1 


n = 


n = 1 


Symbol 


4C 16.49 


N(b) 


1.6 


669.0 


2917.0 


o 






2772.0 


90.1 


1.7 


X 


4C 16.49 


S(c) 


151.9 


0.9 


62.0 


o 






45.4 


0.4 


34.8 


X 



Table 2.8: Reduced values for the rotation measure fits for Figure 2.5. 



2.3.1.3 Sample C: 

3C 16: (Figure C.l) The source shows a strong SW lobe, with a relaxed NE lobe. 
The SW lobe shows a strong depolarisation feature that is narrower than the beam 
size. The strong rotation measure feature is evident in the depolarisation map, 
but not the magnetic field map, indicating that it is not an error in the fitting 
program. No value for the rotation measure was obtained for the NE lobe because 
the polarisation observed was too weak. 



3C 42- (Figure C.2) The core was detected at 4.7 GHz, but was absent at the 
lower frequencies. The source has been observed to lie in a small cluster by de Vries 
et al. (2000). Fernini et al. (1997) observed that the N hotspot was double, but this 
is not evident in my observations which can be attributed to the differences in the 
resolutions of the two observations. 
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Figure 2.5: (a) Map of the rotation measure of the radio source 4C 16.49 (rad m ^) 

between 4710 MHz, 1665 MHz and 1465 MHz. All contours are at ha at 4710 MHz 
(0.8 mJy beam-i)x(-l, 1, 2, 4.. .,1024) with a beam size of 2.2" x 1.8". (b) A plot 
of the polarisation angle against A^, allowing for utt, for the northern lobe of 4C 
16.49 with the dashed lines showing the best fit models. All reasonable utt solutions 
for the 4.7 GHz data are considered and plotted. Data for two small regions, one 
on each side of the jump are plotted with 'o' indicating one side of the jump and 
'x' the other. The jump plotted lies SW of the central intensity contour, (c) Same 
as (b) but for the southern lobe, the region plotted lie either side of the jump south 
of the peak intensity contour. 



3C 46: (Figure C.3) The source has a prominent core at 4710 MHz, but it is 
indistinguishable from the extended lobe at 1452 MHz. 



3C 341 : (Figure C.4) The source is a classic double with a resolved jet-like structure 
running into the SW lobe. The jet is more prominent in the higher frequency 
observations than at the lower frequencies. 
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3C 351: (Figure C.5) The source is an extended and distorted quasar (Bridle et al., 
1994). Both lobes expand out to envelope the core. The NE lobe is highly extended, 
off-axis and shows two very distinct hotspots. The depolarisation increases towards 
the more compact SW lobe supporting the idea that the environment around the 
SW lobe is denser, stopping the expansion seen in the NE lobe. There is evidence 
of a rotation measure ridge in the NE hotspots which corresponds to a narrow ridge 
of depolarisation, but there is no corresponding shift in the magnetic field map. 

3C 457: (Figure C.6) The SW lobe shows a prominent double hotspot. The small 
compact object just south of the SW hotspots is most likely an unrelated background 
object. The inverted core was observed to be present at all frequencies. This source 
has no rotation measure map or magnetic field measure map as I was unable to 
remove all nvr ambiguities from this source. This was due to the small separation of 
observing frequencies around 1.4 GHz and 4.8 GHz, see section 2.3. 

3C 299: (Figure C.7) The source is the least luminous of all the sample members. 
The source shows a large change in rotation measure between the lobes but the 
difference is probably due to the small number of pixels with rotation measure 
information in the NE lobe. 



4C 14-27: (Figure C.8) There is no core detected at any frequency even in a better 
quality map by Leahy & Per ley (1991). 



Chapter 3 

Source properties and trends 



The average spectral index, average depolarisation and average rotation measure are 
given by a, DM and RM respectively. Each of these parameters are the average over 
each individual lobe of a source, a is calculated using oc v~'^ between 1.4 GHz 
and 4.8 GHz, where the 5*1^ are the Total Intensity Flux values given in Tables 3.1 to 
3.3. The average DM is calculated using the percentage polarised flux measurements 
given in Tables 3.1 to 3.3. The average RM is measured directly from the rotation 
measure maps (see appendix A to G) using the AIPS task IMEAN. 

The differential spectral index, differential depolarisation and differential rotation 
measure over the lobe of an individual source are given by da, dDM and dRM 
respectively. All differential properties are taken to be the difference of the averages 
over individual lobes of a source e.g. dDM = DMjofeei - DM;o6e2- The rms variation 
of the rotation measure is defined by orm- 

To test if any property is varying between samples 1 compare the average of the 
property over a sample with the average taken from the other two samples. This 
allows a simple statistical test to be applied to the data. I also average properties 
over samples with equal distributions in redshift (A+B) and equal distributions in 
radio- luminosity (A+G). For example, if depolarisation increases with redshift but 
not with radio-luminosity 1 would expect samples A and B to show similar averages 
(within errors) but to be statistically greater than sample G. In this case I would 
also expect to see the average of sample A-l-G (low radio luminosity sources) to be 
statistically similar to the average of sample B (high radio luminosity sources). A 
more thorough statistical approach using Spearman Rank, Partial Spearman Rank 
and Principal Gomponent is presented in chapter 4. 
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3.1 Polarised Flux and Spectral Index. 

The flux observed from both lobes of aU sources was found to be polarised at levels 
greater than 1%, the only exceptions being 6C 0943+39 and 3C 299, see Tables 3.1 
to 3.3. At lower rcdshifts the polarisation exhibits the largest range from 0.8% in 3C 
299 to 28.0% in 3C 341. This range is not evident in either of the other two samples. 
Statistically the sources at low redshift (sample C) are slightly more polarised than 
those at high redshift, both at 1.4 GHz (sample A: 2.1a, sample B: 2.3a) and at 
4.8 GHz (sample A: l.Scr, sample B: 1.8a). There is no difference between the low 
radio luminosity sources (sample A) and the high luminosity objects (sample B). 
This suggests that percentage polarisation decreases for increasing redshift but is 
less dependent on radio luminosity. 

Table 3.4 shows that there are no significant correlations between spectral index, 
a, and redshift or radio luminosity as suggested by Onuora (1989), Veron et al. 
(1972) and Athreya & Kapahi (1999), which use much larger samples. This suggests 
the trend may be present, but may be masked by my small sample size. However, 
there are trends found with the difference in spectral index between the two lobes. 
Sample B shows a larger average difference in the spectral indices, da between the 
two lobes of a given source than the other samples (sample A: 2.5a, sample C: 3.3a). 
Sample B contains the sources with the highest radio luminosity and so 1 find that 
in my samples the difference in spectral index, da increases with radio luminosity 
rather than with redshift. The trend of the difference in the spectral index between 
the two lobes may be related to the extra luminosity of the 3CRR hotspots compared 
to the 6C/7C hotspots (see section 2.1.2). On average 1 find that the hotspots of 
my sources have shallower spectral indices. The average spectral index integrated 
over the entire sources will therefore depend on the fraction of emission from the 
hotspots compared to the extended lobe, thus creating the observed trend. 



3.2 Rotation measure 

Faraday rotation occurs when polarised light passes through a plasma. AGN's emit 
synchrotron radiation which is almost completely linearly polarised, this can be 
separated into right-handed (RH) and left-handed (LH) circular polarisation which 
have different phase velocities. At a distance dl from a source the two circular 
polarisations will be out of phase by A^, which is given by: 



(l>RH - (t>LH 

2 



(yi^RH) - \f{eLH)^ 



oodl 



2c 



(3.1) 



where 




(3.2) 
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Table 3.1: Properties of the sample A radio source components. Errors are 5% 
or less unless stated otherwise. The spectral indices are the mean values for each 
component, calculated between approximately 4800 MHz and 1465 MHz. 
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Table 3.2: Properties of the sample B radio source components. Errors are 5% 
or less unless stated otherwise. The spectral indices are the mean values for each 
component, calculated between approximately 4800 MHz and 1465 MHz. 
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Table 3.3: Properties of the sample C radio source components. Errors are 5% 
or less unless stated otherwise. The spectral indices are the mean values for each 
component, calculated between approximately 4800 MHz and 1452 MHz. 
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1.11 ± 0.05 
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Table 3.4: Mean properties of the sources averaged over each sample with the as- 
sociated error. Differential spectral index (da) is derived by taking the difference 
in the spectral index between the two lobes of each source and then averaging this 
difference over each sample. Source size is in kpc and radio-luminosity is in W/Hz. 



where oj^ is the phase velocity of the polarisation, ojg is the non relativistic gyro- 
frequency of the electrons in a magnetic field, uj = 2t:v where v is the frequency of 
the observations and c is the speed of light. The total phase difference is then 



= o 2 2 / "^eBwdl (3.3) 

where is the line-of-sight component of the magnetic field, Ue, the column density 
around the source, dl is the path length to the source, e is charge on the electron, 
rrte is the mass of an electron and eo is the electric constant. From section 2.3 
the observed rotation measure, RM, of a source depends on the frequency of the 
observations, RM = A9/X^, hence equation 3.3 becomes: 

'''' = s^J^ f^-^) 

where RM is in rad m~^. 



Observationally it is not possible to measure the electron density or the magnetic 
field strength without some prior assumptions, thus the rotation measure must be 
calculated from the observed polarisation angle. The rotation measure a source 
displays is related to the degree of rotation of the polarisation position angle over 
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a set frequency range, see equation 2.2. To gain an accurate value of the rotation 
measure three frequencies (1.4 GHz, 1.6 GHz and 4.8 GHz) are used to overcome 
the nn ambiguities discussed in chapter 2.3 that arise when fitting to the observed 
polarisation angles (Simard-Normandin et al., 1981; Rudnick et al., 1983). 

The observed RM and the degree of polarisation in a source may be caused by the 
presence of plasma either inside the radio source itself (internal depolarisation) or by 
a Faraday screen in between the source and the observer (external depolarisation). 
In the latter case the screen may be local to the radio source or within the Galaxy, 
or both. Only in the case of an external Faraday screen local to the radio source do 
the measurements contain any information on the source environment. 

The average RM observed in my sources is consistent with a Galactic origin 
(Leahy, 1987). This is also consistent with the absence of any significant differences 
of RM between the samples (see Table 3.5). However, large variations of RM are 
observed within individual lobes on small angular scales, e.g. dRM and (Trm, are 
probably caused by a Faraday screen local to the source (Leahy, 1987) and therefore 
must be corrected for the source redshift by multiplying by a factor {1 + z)"^. This 
allows a valid comparison between sources. The variation of RM on angular scales 
of order 10s of arcseconds, i.e. between the two lobes of a source, dRM, may still 
be somewhat influenced by the Galactic Faraday screen. Nevertheless, the large 
variations of RM found on arcsecond scales measured by (Trm suggest an origin 
local to the source. There is no statistically significant trend of dRM with respect 
to redshift or radio luminosity as is seen in Table 3.5. This is again consistent with 
a Galactic origin of the RM properties of the sources on large scales. On small 
angular scales the variation of RM as measured by its rms variation, a^M, shows a 
significant trend between the low redshift sample (C) and the high redshift samples 
(A: 2.6a, B: 3.3a). The exclusion of 3C 457 due to n7r-ambiguities in the rotation 
measure (see Section 2.3) from the rotation measure averages could bias the sample 
C results towards a small value of auM- In principle, similar ambiguities could also 
affect other sources in sample C. However, as stated above the rotation measure is 
found to vary smoothly across the other sources in this sample. The exclusion of 
one source will not remove the trend found here. There is no significant difference 
between the low and high radio luminosity samples A and B, suggesting that at 
least the range of rotation measures on small angular scales in a source does depend 
on the source redshift and not the source radio luminosity. 



3.3 Depolarisation 



The observed depolarisation depends on the amount of polarised flux at two frequen- 
cies but also on the total intensity flux at the same two frequencies. It is defined 
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Table 3.5: Mean properties of the sources averaged over each sample with the asso- 
ciated error. Differential rotation measure (dRM) is derived by taking the difference 

of the rotation measure between the two lobes of each source and then averaging 
this difference over each sample. a^M is the variation of the RM over the source. 
Properties in italics are in the source frame of reference. 



as: 

nj\^4.8GHz _ -P-^4.8GHz _ "^6cm /„ 

^^^hAGHz - - — > ['^■^) 

-r-ri.4GHz "7.20 cm 

where the PF^ is the fractional polarisation, (polarised fiux)/(total intensity flux), 
at a given frequency, and mx is the percentage polarisation at an observed wave- 
length A. 

The observed dcpolarisation also depends on the distribution of the Faraday 
depths covered by the projected area of the telescope beam, i.e. the variation of 
RM on the smallest angular scales. If the depolarisation is caused by a local but 
external Faraday screen and the distribution of Faraday depths in this screen is 
Gaussian with standard deviation A, then (e.g. Burn, 1966) 



mx = moexp{-2A^[X/{l + z)f}, (3.6) 



where, mo is the initial percentage polarisation before any depolarisation and z is 
the redshift of the source. Since I measure mx at two observing frequencies (1.4 GHz 
and 4.8 GHz), equation (3.6) can be solved for A as a function of the depolarisation 
measure. 
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A 



2{Xi 



4 GHz 



X4 

^4.8 GHz 



radm 



(3.7) 



I now drop the GHz subscript. 



For each source, it is possible to compare the value of A as derived from the 
measured depolarisation with the observed rms of the rotation measure, aRM- If 
the Faraday dispersion is less than ajiM, then the observations are consistent with 
an external Faraday screen (Garrington & Conway, 1991). Figure 3.1 displays the 




10 20 30 40 50 

A 



Figure 3.1: Plot of the Faraday dispersion. A, against the rms of the rotation 
measure for each source. 

Faraday dispersion. A, for each lobe of the sources against the rms of the rotation 
measures observed. It is evident from the plot that the value of aRM > A for most 
components, see Tables 3.6 to 3.8 but there are a few sources where this is not 
the case. However, these components belong to sources where the depolarisation or 
rotation measure is only determined reliably for a few pixels so an accurate value 
is not obtainable for a^M or A. There is little correlation between A and a^M- 
This is a strong indicator that the Faraday medium responsible for variations of 
RM on small angular scales, and thus for the polarisation properties of the sources, 
is consistent with being external but local to the sources'. 



Comparing the average depolarisation of individual samples with each other I 
find only very weak trends with redshift or radio luminosity. When the samples are 
averaged together I find a trend in depolarisation with redshift but none with radio 
luminosity. Samples A+C (low radio luminosity) are statistically identical to sample 
B (high radio luminosity), suggesting that there is no trend with radio luminosity 
in my sources. By considering the averaged depolarisation of samples A+B (high 
redshift) compared with that of sample C (low redshift) there is a weak trend, 1.7a, 
with redshift. This trend is echoed in the dDM values {1.9a). This may confirm 
the results of Kronberg et al. (1972): redshift is the dominant factor compared to 
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Source Component RM DMfl Average grm 

A 









(rad m ^) 




(rad m ^) 


(rad m ^) 


6C 


0943+37 


W 


1.6 ± 2 


2.19 


15.01 


16.1 






E 


-19.1 ± 6 


3.86 


19.7 


101.4 


6C 


1011+36 


N 


30.8 ± 5 


1.53 


11.1 


12.8 






S 


12.6 ± 4 


1.06 


4.1 


10.1 


6C 


1018+37 


NE 


0.85+ 3 


1.23 


7.71 


6.1 






SW 


11.2 ± 2 


2.41 


15.9 


2.9 


6C 


1129+37 


NW 


-19.3 ± 4 


2.51 


16.3 


54.0 






SE 


0.08+ 3 


5.09 


21.61 


22.0 


6C 


1256+36 


NE 


5.9 ± 3 


1.35 


9.3 


21.8 






SW 


15.4 ± 3 


1.01 


1.7 


14.0 


6C 


1257+36 


NW 


-115.3 ± 9 


1.27 


8.3 


10.0 






SE 


-115.6 ±10 


1.88 


13.5 


16.0 


7C 


1745+642 


N 




1.85 


13.3 








Core 


65.4 ± 4 


1.69 


12.3 


20.9 






S 


12.8 ± 3 


0.91 




3.1 


7C 


1801+690 


N 


44.8 ± 3 


1.11 


5.5 


16.0 






Core 


30.3 ± 3 


0.45 




1.80 






S 


20.8 ± 2 


1.39 


9.7 


10.0 


7C 


1813+684 


NE 


13.9 ± 3 


0.81 




85.0 






SW 


-68.4 ± 7 


1.14 


6.1 


41.7 



Table 3.6: Properties of the sample A radio source components. Errors are 5% 
or less unless stated otherwise. The depolarisation measures are the mean values 
of the ratio of the fractional polarisation between approximately 4800 MHz and 
1465 MHz for each component. The rotation measures are the mean values between 
approximately 4800 MHz, 1665 MHz and 1465 MHz and are quoted in the observer's 
frame of reference. The Faraday dispersion, A, is given in equation 3.7. (Trm is the 
rms in the rotation measure. All mean values take into account pixels above Barms 
at 4.8 GHz and above 3(7^^, at 1.4 GHz. 
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Source Component RM DMf| Average o^rm 

A 







(rad m 






(rad m ^) 


(rad m ^) 


3C 65 


W 


-82.6± 


7 


3.57 


19.1 


21.4 




E 


-86.1± 


6 


1.28 


8.4 


9.7 


3C 252 


NW 


15. 7± 


6 


1.08 


4.7 


20.2 




SE 


58.5± 


6 


2.07 


14.5 


45.4 


3C 267 


E 


-9.6± 


3 


1.31 


8.8 


24.4 




W 


-21.5± 


3 


1.06 


4.1 


23.4 


3C 280 


E 


-37.7± 


8 


1.82 


13.1 


20.1 




W 


-7.5± 


4 


1.54 


11.1 


1.5 


3C 324 


NE 


22. 1± 


4 


1.66 


12.1 


23.0 




SW 


43.0± 


5 


1.95 


13.9 


16.8 


4C 16.49 


N 


-4.3± 


5 


0.64 




56.4 




Jet 


30.1± 


4 


2.01 


14.2 


57.2 




Core 






0.77 








S 


0.9± 


3 


1.09 


5.0 


56.4 


3C 68.1 


N 


-26.6± 


5 


1.23 


7.7 


27.7 




S 


57.9± 


2 


1.73 


12.5 


59.8 


3C 265 


NW 


42.2± 


6 


1.72 


12.5 


21.9 




SE 


32.8± 


3 


1.48 


10.6 


22.8 


3C 268.1 


E 


21. 7± 


5 


1.47 


10.5 


10.8 




W 


26.8± 


6 


1.74 


12.6 


13.3 



Table 3.7: As Table 3.6 but for sample B. 
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Source 


Component 


RM 
(rad m~^) 


DMf;| 


Average 
A 

(rad m"^) 


<yRM 

(rad m~^) 


3C 42 


NW 


-2.4 ±5 


1.02 


2.4 


14.1 




SE 


5.0 ±5 


1.01 


1.7 


14.2 


4C 14.27 


NW 


-13.0 ±4 


1.41 


9.9 


5.4 




SE 


-17.3 ±5 


1.03 


2.9 


4.8 


3C 46 


NE 


-4.8 ±5 


1.27 


16.9 


3.9 




SW 


-2.9 ±1 


1.14 


12.5 


4. 


3C 457 


NE 




1.01 


2.1 






SW 




1.19 


5.6 




3C 351 


NE 


1.00±5 


0.84 




13. 




Diffuse 


-8.7 ±1 


1.93 


13.7 


12. 




Core 


-0.53±2 


2.27 


15.3 


5. 




SW 


4.4 ±3 


2.56 


16.4 


10.8 


3C 341 


NE 


20.3 ±5 


1.43 


10.1 


10.0 




SW 


18.2 ±5 


1.53 


11.1 


12.0 


3C 299 


NE 


-126.3 ±8 


1.84 


13.2 


75.9 




SW 


16.0 ±5 


1.19 


7.10 


6.9 


3C 16 


NE 




2.20 


15.0 






SW 


-4.3 ±3 


1.79 


12.9 


17.4 



Table 3.8: As Table 3.6 but for sample C. 
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radio luminosity in determining the depolarisation properties of a source but since 
the significance levels are low it is not possible to say with any confidence. 



Average 




Sample 




Prop. 


A 


B 


C 


z 


1.06 ± 0.04 


1.11 ± 0.05 


0.41 ± 0.01 


Pl51 


(1.02 ± 0.06) X 10^7 


(1.01 ± 0.07)xl028 


(8.02 ± 0.40) X 10^6 


D 


220.34 ± 46.40 


262.69 ± 55.41 


399.57 ± 122.50 


DM 


1.82 ± 0.32 


1.61 ± 0.13 


1.42 ± 0.12 


(IBM 


0.93 ± 0.26 


0.61 ± 0.22 


0.13 ± 0.18 




Average 


Sample 




Prop. 


A+B 


A+C 




z 


1.08 ± 0.03 


0.75 ± 0.08 




Pl51 


(5.57 ± 1.16) X 10^7 


(8.97 ± 0.77) X 10^6 




D 


241.51 ± 35.43 


304.69 ± 64.50 




DM 


1.71 ± 0.17 


1.63 ± 0.18 




dDM 


0.77 ± 0.17 


0.69 ± 0.17 



Table 3.9: Mean properties of the sources averaged over each sample with the asso- 
ciated error. Difi^ercntial depolarisation (dDM) is derived by taking the difference 
of the depolarisation between the two lobes of each source and then averaging this 
difference over each sample. 



Strom (1973); Strom & Jagers (1988); Pedelty et al. (1989); Ishwara-Chandra 
et al. (1998) find an anti-correlation of depolarisation with linear size. Figure 3.2 
shows that there is a weak anti-correlation between physical source size and depo- 
larisation. However, this trend is due to the two largest sources in sample C, 3C 
46 and 3C 457. Removing these sources yields an average depolarisation measure 
of 1.50 ± 0.15 for sample C, which is not significantly different from the average 
with these two sources included. So I can rule out the possibility that the larger 
depolarisation at high redshift is caused by selecting preferentially smaller sources 
at high redshift. 



3.3.1 Laing— Garrington effect 

Assuming that a radio source does not lie flat to the observers line-of-sight and 
that both lobes are embedded in the same Faraday medium then one lobe will lie 
further away. The polarised flux from this lobe will have to travel through a larger 
path length and will become more depolarised than the nearer lobe. This creates 
an asymmetry in the depolarisation that is known as the Laing-Garrington effect. 
8 sources show dDM > 1 over the source and a further 6 show 0.5 < dDM < 1. 
This implies that 14 sources out of 26 sources show a significant asymmetry in 
the depolarisation of their lobes. I would expect a proportion of the sources to be 
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Figure 3.2: Average depolarisation against projected linear size (kpc) for all 3 sam- 
ples. Symbols as in Figure 2.1. Figure (a) assumes Ho = 75 kms~^Mpc~^, and 



0.5, Q\ = 0. Figure (b) assumes Ho = 75 kms ^Mpc ^, and Qr, 



0.35, 



0.65. 



observed at angles considerably smaller than 90° to my line-of-siglit. Therefore it 
is not surprising that so many sources are found to be candidates for the Laing- 
Garrington effect. 



3.3.2 Depolarisation shadows 



Depolarisation shadows (regions where the depolarisation is appreciably greater than 
in the surrounding area) are seen in a few sources, e.g. Figure 3.3, and may be real 
features. These were first found by Fomalont et al. (1989) in a study of Fornax A. 
Depolarisation shadows can be caused by the parent galaxy as in the case of 3C 
324 or by an external galaxy in the foreground of the source, causing a depolarising 
silhouette (Best et al., 1997). Several sources show signs of depolarisation shadows 
in at least one of the lobes. Of these, 6C 1256+36 was observed to lie in a cluster 
by Roche et al. (1998), as were 3C 65 and 3C 324 (Best, 2000). 3C 324 has been 
observed by Best et al. (1998) in the sub-mm (850/im) and found a large dust mass 
centred around 3C 324. The host galaxy was shown to be the cause of the very 
strong depolarisation, (Best et al., 1998). 
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Figure 3.3: Depolarisation map of the radio source 3C 324 between 4848 MHz 
and 1465 MHz. All contours are at 5a at 4848 MHz (1.5 mJy bearn"^) x(-l, 1, 2, 
4.. .,1024) with a beam size of 2.2" x 1.6". 

3.4 Using Burn's law 
3.4.1 Polarisation 

So far I have only considered the percentage polarisations of a source measured in 
the observing frame. Variations of the RM on small angular scales which determine 
the degree of polarisation are caused in Faraday screens local to the sources. The 
trend with redshift may therefore simply reflect the different shifts of the observ- 
ing frequency in the source rest-frame for sources at low and high redshift. Using 
Burn's law in the form of equations (3.6) and (3.7) it is possible to determine, for 
each source, the percentage polarisation expected to be observable at a frequency 
corresponding to a wavelength of 5 cm in its rest-frame. This wavelength was chosen 
as it lies close to the wavelengths already observed in all the sources. The results 
for individual sources are presented in Table 3. 10 and the sample averages are sum- 
marised in Table 3.11. Again sources at low redshift (sample C) are slightly more 
polarised than sources at high redshift (sample A: 1.9(7, sample B: 2.1a). As before, 
there is no trend with radio luminosity, indicating that the trend with redshift is 
dominant. This is not caused by pure Doppler shifts of the observing frequencies. 



3.4.2 Depolarisation 

Cosmological Doppler shifts of the observing frequencies influence the trend of DM 
with redshift, and in fact the true trend is stronger than that naively observed. To 
demonstrate this equation (3.7) can be used to derive the standard deviation of 
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Source 


z 




DM^=i 




Polarisation 

/^rest — ^ ^^^^ 


6C 


0943+37 


1.04 


74.12 


3.31 


1.98 


3.10 


6C 


1018+37 


0.81 


42.95 


1.49 


0.52 


8.32 


6C 


1011+36 


1.04 


36.11 


1.33 


0.66 


5.63 


6C 


1129+37 


1.06 


83.07 


4.49 


3.43 


3.09 


6C 


1256+36 


1.07 


29.53 


1.21 


0.40 


8.26 


6C 


1257+36 


1.00 


45.51 


1.58 


0.61 


9.39 


7C 


1801+690 


1.27 


41.24 


1.45 


0.54 


3.33 


7C 


1813+684 


1.03 


6.96 


1.01 


0.25 


8.47 


7C 


1745+642 


1.23 


47.81 


1.65 


1.36 


7.32 


3C 


65 


1.18 


75.33 


3.50 


4.61 


6.26 


3C 


252 


1.11 


50.66 


1.76 


1.36 


6.36 


3C 


267 


1.14 


32.34 


1.26 


0.35 


5.04 


3C 


280 


1.00 


48.79 


1.68 


0.28 


5.48 


3C 


324 


1.21 


54.34 


2.42 


0.58 


4.74 


3C 


265 


0.81 


35.56 


1.37 


0.14 


5.05 


3C 


268.1 


0.97 


42.13 


1.57 


0.25 


3.07 


3C 


68.1 


1.24 


45.11 


1.85 


0.98 


10.20 


4C 


16.49 


1.29 


31.04 


1.23 


0.47 


6.76 


3C 


42 


0.40 


4.59 


1.00 


0.01 


9.85 


4C 


14.27 


0.39 


14.34 


1.05 


0.08 


9.08 


3C 


46 


0.44 


15.06 


1.05 


0.03 


12.45 


3C 


457 


0.43 


10.50 


1.03 


0.04 


12.88 


3C 


351 


0.37 


23.82 


1.14 


0.15 


8.59 


3C 


341 


0.45 


21.91 


1.11 


0.02 


18.57 


3C 


16 


0.41 


27.54 


1.19 


0.06 


6.64 


3C 


299 


0.37 


20.20 


1.10 


0.10 


1.52 



Table 3.10: Recalculated average depolarisation and dDM for each source if it was 
located at z = 1 and average percentage polarisation of all sources if it was emitted 
at 5 cm in the rest frame. 

Faraday depths, A, for each source. By setting z = 1 it is possible to then rescale 
all the depolarisations, DM^ ^, to the same rcdshift. This allows all 3 samples to 
be compared without any bias due to pure redshift effects (see Tables 3.10 and 
3.11). If there was no intrinsic difference between the high- redshift and low-redshift 
samples then I would expect these values to be consistent with each other. This 
is evidently not the case. The high redshift samples are, on average, significantly 
more depolarised (sample A: 2.2a, sample B: 3.2ct) than their low-redshift counter- 
parts (sample C). Comparing sample A with sample B 1 find no trend with radio 
luminosity. However, a note of caution must be issued as the corrections applied 
use Burn's law and may actually be too large. Considering the prccorrccted and the 
corrected values together it is obvious that there is a connection between redshift 
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and depolarisation but there is no significant trend of depolarisation witli radio lu- 
minosity. There is also a connection between the difference in the depolarisation, 
dDM, and the redshift of the source, but no significant trend of the difference in 
the depolarisation with the radio luminosity of the source. As noted in section 2.3 
regions with signal-to-noise < 3a were blanked in the map production. In individual 
sources blanking of low S/N regions in the polarisation maps will cause the measured 
depolarisation to be underestimated. Sources in sample A are more affected by this 
problem than objects in the other samples. Therefore I probably underestimate the 
average depolarisation in sample A implying that the trend with redshift could be 
even stronger than my findings suggest. 



Average 






Sample 






Prop. 


A 


B 


C 


A+B 


A+B 


DM,=i 


1.95 ± 0.39 


1.85 ± 0.24 


1.08 ± 0.02 


1.90 ± 0.22 


1.54 ± 0.23 


dDM^=i 


1.08 ± 0.34 


1.00 ± 0.47 


0.06 ± 0.02 


1.04 ± 0.28 


0.60 ± 0.22 


^^\rest=5 cm 


6.32 ± 0.86 


5.88 ± 0.65 


9.82 ± 1.85 


6.10 ± 0.53 


7.97 ± 1.04 



Table 3.11: Average DM and dDM now shifted so that all the measurements are 
taken at z = 1. Average polarisation of all the sources, shifted to a common rest 
frame wavelength of 5 cm using equation 3.7. 



The trends of percentage polarisation and of depolarisation with redshift are 
probably related in the sense that a lower degree of depolarisation at low redshift 
also leads to a higher observed degree of polarisation. Clearly a variation of the 
initial polarisation, mo, with redshift would lead to variations of the observed mx 
independent of the properties of any external Faraday screen. Therefore both trends 
could also be caused by a significantly higher level of mo for sources at low redshift 
(sample C). Using equation (3.6) 1 find mo = 9.3 ±1.1 for sample A, mo = 9.1 ± 1.0 
for sample B and mo = 13.5 ± 2.5 for sample C. The uncertainties associated with 
the use of Burn's law in extrapolating from my observations to A = are large. 
There is no difference between the average initial polarisation of the sources in 
samples A and B. The difference found for mo at low redshift (sample C) compared 
to high redshift (samples A and B) is small compared to the difference found for the 
depolarisation comparing the same samples. This suggests that the differences in 
percentage polarisation and depolarisation are due to variation with redshift of the 
Faraday screens local to the sources rather than to differences in the initial degree 
of polarisation. However note, that the variation of mo is not significantly smaller 
than the trend of percentage polarisation with redshift. 

Burn's law predicts a steep decrease of percentage polarisation with increasing 
observing frequency. Although the decrease may well be 'softened' by geometrical 
and other effects in more realistic source models (Laing, 1984), the value of DM may 
be small for sources in which the observing frequencies are lower than the frequency 
at which strong depolarisation occurs. For such sources I would expect to measure a 
low value of DM associated with a low percentage polarisation. Akujor & Garrington 
(1995) show that one of my sources, 3C 299 (sample C), is strongly depolarised 
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between 1.6 GHz and 8.4 GHz with almost all of the depolarisation taking place 
between 4.8 GHz and 8.4 GHz. I measure only a very low percentage polarisation 
for this source and the average value of DM^;4 = 1.5 is also lower than DM^;g ~ 4 
as measured by Akujor & Garrington (1995). If many of the sources at low redshift 
are affected by strong depolarisation at frequencies higher than 4.8 GHz, then this 
may cause the trend of DM with redshift noted above. The absence of any other 
sources with very low degrees of polarisation combined with a low value for DM, 
at least in the low redshift sample C, argues against this bias. In fact, Tabara & 
Inoue (1980) show that the sources 3C 42, 3C 46, 3C 68.1, 3C 265, 3C 267, 3C 
324 and 3C 341 depolarise strongly only at frequencies lower than my observing 
frequencies. 3C 16, 3C 65, 3G 252, 3G 268.1 and 3C 280 do depolarise strongly 
between 1.4 GHz and 4.8 GHz. In all the sources mentioned in Tabara & Inoue 
(1980) none have strong depolarisation at frequencies higher than 4.8 GHz. There 
is no information for 4C 14.27, 4C 16.49, 3C 351 and 3C 457. Those sources which 
do depolarise strongly between my observing frequencies, 1.4 and 4.8 GHz, could in 
principle have inaccurate values of the rotation measure because of this. However, 
the pixels containing most depolarised regions of the source are likely to have been 
blanked because of insufficient signal-to-noise in their polarisation at 1.4 GHz; the 
rotation measure is determined from the unblanked (less depolarised) regions of the 
source, and will therefore be reliable. 



3.5 Discussion 



In this chapter I present the complete data set of the three samples of radio galaxies 
and radio-loud quasars. The three samples were defined such that two of them 
overlap in redshift and two have similar radio luminosities. Allowing the effects of 
redshift and radio luminosity on various source properties to be analysed. 

There is little correlation between A and a^M, suggesting that the Faraday 
medium responsible for variations of RM on small angular scales, is consistent with 
being external but local to the sources. This implies that any trends with radio lumi- 
nosity and/or redshift reflect changes of the source environment depending on these 
quantities. There is also little correlation between rotation measure and redshift or 
radio luminosity which is consistent with a Galactic origin of the RM properties of 
the sources on large scales. However, I find that the rms fluctuations of the rotation 
measure correlate with redshift but not radio luminosity to a confldence level of 
> 99.9%, determined in the sources' frame of reference. 

I flnd that the polarisation of a source anti-correlates with its redshift but is in- 
dependent of its radio luminosity, resulting in the low redshift sample having much 
higher degrees of polarisation, in general. I also detect higher degrees of depolari- 
sation in the high redshift samples (A and B) compared to sources at lower redshift 
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(sample C). This suggests that depolarisation is correlated with redshift. These 
two resuhs are probably related in that lower depolarisation at low redshift leads 
to both lower depolarisation measurements and also higher degrees of observed po- 
larisation. According to Burn's law (Burn, 1966), this implies an increase in the 
source environments of either the plasma density or the magnetic field strength or 
both with redshift. Such an interpretation is also supported by the increased depo- 
larisation asymmetry of sources at high redshift compared with their low redshift 
counterparts. 

The findings on the rotation measurements and polarisation properties of my 
sources are indicative of an increase of the density and/or the strength of the mag- 
netic field in the source environments with increasing redshift. 

I do not find a strong correlation between the projected sizes of the sources and 
their depolarisation measure as found by Strom (1973); Strom & Jagers (1988); 
Pedelty et al. (1989); Ishwara-Chandra et al. (1998). 

I find no correlation between the spectral index and redshift or radio luminosity. 
However, I do find the difference in the spectral index, across individual sources, in- 
creases for increasing radio luminosity of the source and also with increasing redshift 
of a source. 

The correlations found are only general trends within the samples. To completely 
break the z-P degeneracy a much more rigorous and detailed approach to the sta- 
tistical analysis is needed. In the next chapter I will investigate the observational 
properties using Principal Component Analysis and Spearman Rank and determine 
the significance of the results presented here. 



Chapter 4 

Analysing the observations 



In this chapter I investigate the trends and correlations of a number of observational 
parameters (spectral index, rotation measure and depolarisation measure) with the 
'fundamental' parameters, low frequency radio luminosity, redshift and the physical 
size of the radio structure. The observational parameters mainly constrain the 
properties of the gas in the vicinity of the radio structure which acts as a Faraday 
screen, see sections 3.2 and 3.3. 

In my statistical analysis I use Spearman Rank as well as Partial Spearman Rank 
and Principal Component Analysis techniques. Thus I am able to assess whether any 
observed correlations of parameters with redshift and radio luminosity are significant 
in their own right or whether they simply arise from a Malmquist bias. 

The 3 samples can also be used to determine if there are any asymmetries present 
between the lobes of sources and if so, what causes these asymmetries. Traditionally 
this has been done with sources in which jets arc detected so that the asymmetries 
between the jet side and counter-jet side can be studied, (e.g. Garrington & Conway, 
1991; Dennett-Thorpe et al., 1997; Leahy, 1987). This type of study allows a direct 
insight into the orientation of a source. In my sample only a few sources have a well 
defined jet at cither 4.8 GHz or 1.4 GHz, so this type of analysis is not possible. 
However, I compare the brighter lobe to the fainter lobe to determine if there are any 
asymmetries present in the samples and if possible, explain the underlying physical 
cause(s) of the asymmetries. 
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4.1 Statistics 



4.1.1 Spearman Rank 



The Spearman Rank test is a non-parametric correlation test assigning a rank to 
given source properties, X and Y and then performing a correlation on the rank. 
The statistic is given by 

r^^ = ^"^yj (4.1) 



•n? — n 



where d is the difference between the rank in the X direction and the rank in the Y 
direction and txy has a Student-t distribution: 



D^. = ril^ (4.2) 



where n is the number of sources. 



The null hypothesis states that no correlation is present between two properties 
X and Y, corresponding to rxY=0 and a (anti-) correlation corresponds to (-)l. A 
correlation with a confidence level above 95% will show a Student-t value greater 
than 2 for my sample size (Riley et al., 1997). 

Partial Spearman Rank is used to determine if a correlation between two variables 
(X, Y) depends on the presence of a third, A (Macklin, 1982). The null hypotheses 
are: 



• the A-X correlation is entirely due to the X-Y and the Y-A independent cor- 
relations, and 

• the Y-A correlation is entirely due to the X-Y and the X-A independent cor- 
relations 



The Partial Spearman Rank correlation coefficient is given by: 

V[(1-4y)(1-'^ay)] 

where tax is the normal Spearman Rank correlation coefficient described above. 
The significance level for the A-X correlation, independent of Y is 

OAx.Y = i("-4)4ta(i±^). (4.4) 

^ V -'- ~ ^AX,Y / 
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To fully determine the relationships between X,Y and A, the triplet tax.y, ''"ay.x 
and rxY,A must be calculated and their corresponding significance levels found. By 
comparing this triplet it is possible to determine which relationship, if any, between 
X, Y and A dominates. 



4.1.2 Principal Component Analysis, PCA 



Principal Component Analysis, PCA, is a multi-variate statistical test which is de- 
scribed in detail by Deeming (1964) and Efstathiou & Fall (1984). PCA is a linear 
self-orthogonal transformation from an original set of n objects and m attributes 
forming a n x m matrix, with zero mean and unit variance, to a new set of parame- 
ters, known as principal components. The principal components of the new dataset 
are all independent of each other and hence orthogonal. The first component de- 
scribes the largest variation in the data and has, by definition, the largest eigenvalue 
of the n X m matrix. In physical terms the magnitude of the eigenvalue determines 
what fraction of the variance in the data any correlation describes. Thus a strong 
correlation will be present in the first eigenvector and will only strongly reverse 
in the last. PCA in its simplest terms is an eigenvector-eigenvalue problem on a 
transformed, diagonalised and standardised set of variables, n objects will create 
an eigenvector-eigenvalue problem in n — 1 dimensions. To avoid problems with the 
interpretation of my results 1 do not use more than 4 parameters in any part of the 
analysis. 



4.2 Fundamental Parameters 



In the following I investigate the relationship of the observational properties of my 
sources with the 'fundamental' properties redshift, z, radio luminosity at 151 MHz, 
Pi5i, and physical size, Dgource- For this I first need to understand the relations 
between these three fundamental parameters. Table 4.1 contains the associated 
Spearman Rank results. Redshift and radio luminosity are highly correlated, to a 
significance of ~ 99.9% but this is simply due to the way the samples were selected. 
Using the statistical techniques described above this correlation can be isolated from 
other correlations. In fact, it provides a 'bench mark' for other correlations. 



Parameters 


rs value 


t value 


Z Pl51 


0.69846 


4.78 


Z D source 


-0.22940 


-1.15 


Pl51 ^source 


-0.13026 


-0.64 



Table 4.1: Spearman rank values for z, P151 and D 
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Figure 4.1: (a - left) Redsliift against average source sizes in kpc and (b - right) 
with individual lobe sizes in kpc. Sample A is represented by 'x', sample B by 'o' 
and sample C by '*'. 



There is also a weak z - D source anti-correlation, but as Table 4.1 demonstrates 
this anti-correlation is not significant. Figure 4.1(a) shows how the source size 
is distributed over the redshift range of my samples. The presence of the two low 
redshift giants, 3C 46 and 3C 457 are the cause of the weak correlation. By removing 
these two giants the anti-correlation completely disappears. However, as the anti- 
correlation is below 95% significant, even with the two giants included in the sample, 
any correlations found with size will not be biased towards lower redshift sources. 
Figure 4.1(b) demonstrates that the sizes of individual lobes are evenly distributed 
over the redshift range of the three samples and is also not unduly affected by the 
presence of the two giants. 

Partial Spearman Rank was used to determine the dependencies between z, P151 

and smi.rr.p.' 



0.04, 

rzPl51,D source = 0-69, 
fPlilDsourccZ = -0.20, 



D = 0.20 
D = 4.00 
D = 0.93 



^source auti-correlatiou. 



Using Partial Spearman Rank I find no independent z 
proving that my sample is not biased towards larger sources at lower redshifts, 
allowing a fair comparison with higher redshift sources. Interestingly Table 4.1 in- 
dicates that there is no relationship between a sources' radio-luminosity and its size 
and this result is confirmed by the Partial Spearman Rank results at any fixed red- 
shift. In a study by Kapahi (1989) it was found that Dsource oc (1 + ^)-3±o.5 ^gj^g 
= 1, but at constant redshift it was found that Dsource oc p0-3±0-i. Singal (1993) 
also found that for his sample of 789 sources, size predominantly correlated with 
radio- luminosity and was only a very weak function of redshift. However it is inter- 
esting to note that Barthel & Miley (1988) found that D oc (1 + ;2)-i-5±i.4p-o.o3±o.3 
using a sample of steep spectrum radio quasars. My results are consistent with 
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the results from Barthel & Milcy (1988). These results demonstrate that it is very 
complicated to disentangle the redshift, radio-luminosity and size correlation and 
my sample is simply too small to add anything significant to this well studied field. 

In Tables 4.2 to 4.7 the relationships between the observational parameters with 
the fundamental parameters are given. 



4.3 Observational Parameters 



In order to compare the differential observational parameters e.g. dDM^.* with 
averaged properties like DM, all averaged properties are taken over the entire source 
as opposed to individual lobes. This is the only method possible as for each source 
I only have one measurement of the differential properties and two measurements 
of the average properties (taken from each lobe). In the case of the RM properties 
(RM, dRM^ & (JRM^) I only use 23 of the 26 sources in the analysis. This is simply 
due to the fact that 7C 1745+642, 3C 16 and 3C 457 do not have rehable RM 
information in at least one of their lobes. 



4.3.1 Spectral index 



Parameters rs value 


t value 


Parameters 


rs value 


t value 


as z -0.12615 


-0.62 


da z 


0.15487 


0.77 


as T^saurce 0.20889 


1.05 


da ^source 


-0.20957 


-1.05 


as Pi5i 0.03932 


0.19 


da Pi5i 


0.26632 


1.35 



Table 4.2: Spearman rank values for all spectral index parameters where the sub- 
script s indicates that the spectral index is averaged over both lobes in a source. 

There are no strong correlations between spectral index, a^, and any of the fun- 
damental parameters as the Spearman Rank results in Table 4.2 demonstrates. The 
PCA analysis confirms these findings (Table 4.3). The difference in the spectral 
index between the lobes, da, also shows httle correlation with any of the funda- 
mental parameters. Veron et al. (1972) and Onuora (1989) find a strong spectral 
index-radio luminosity correlation, but no corresponding correlation with redshift. 
Interestingly Athreya & Kapahi (1999) find a strong correlation of spectral index 
with redshift, but no significant correlation with radio luminosity when they anal- 
ysed the MRC/ IJy sample. Although I do not find strong support for either trend, 
this could simply be clue to niy small sample size. 

*dDMj, = IDMj-j^^^j — DM^j^^^jI, see also the previous chapter. 
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Parameter 


1 


2 


3 


4 


z 


0.6354 


-0.2475 


0.0568 


-0.7292 


P 


0.5781 


-0.1595 


0.5304 


0.5992 


D 


-0.4690 


-0.1168 


0.8206 


-0.3051 


as 


-0.2053 


-0.9485 


-0.2050 


0.1271 


Eigenvalue 


44.6% 


24.6% 


20.5% 


10.2 % 


z 


0.5565 


-0.0810 


0.5422 


-0.6243 


P 


0.5487 


0.4404 


0.2629 


0.6602 


D 


-0.4038 


0.8418 


0.1920 


-0.3024 


da 


0.4756 


0.3014 


-0.7746 


-0.2879 


Eigenvalue 


51.4% 


21.3% 


17.7% 


9.5 % 



Table 4.3: Eigenvectors and Eigenvalues for spectral index. 



A recent study by Blundell, Rawlings & Willott (1999) using the 3CRR, 6CE 
and 7C samples found that a anti-correlated not only with radio-luminosity, but 
also with source size. They found that in general, larger sources had a steeper 
spectral index. Again the fact that I do not find this trend could simply be due to 
my much smaller sample size. 



4.3.2 Depolarisation 

In the previous chapter 1 motivated the adjustment of the measured depolarisation to 
a common redshift, z = 1, using the Burn^ theory of depolarisation (Burn, 1966), for 
all sources. There is httle difference between the redshift and the radio luminosity 
correlations with the average adjusted depolarisation, taken over the source as a 
whole, DM^, (see Table 4.4) as both correlations show significance levels greater 
than 99.9%. There is also an indication of a DMz-Dgource anti-correlation, but this 
is found not to be significant. The Partial Spearman Rank test is used to determine 
if the P151-DM2; correlation can be explained by the independent z-DM^ and Z-P151 
correlations. 

r.DM.,Pi5i = 0.40, D = 2.00 

rzP^,„DM. = 0.44, D = 2.21 

rp,,,DM.,. = 0.41, D = 2.04 

It is evident that there is still little difference between the redshift and radio- 
luminosity correlations. This is confirmed using the PCA results in Table 4.5 (top). 

t As noted in the previous chapter the Burn method for shifting the depolarisation may resuh in 
an over-estimation of the depolarisation. In the next chapter I discuss other methods for shifting 
the depolarisation and determine what effect this has on the correlations presented here. 
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Figure 4.2: (a - left) Redsliift against the average depolarisation of a source, (b - 
right) Radio- luminosity against the average depolarisation. The depolarisation has 
been shifted to a common redshift, z=l. Symbols are as in Figure 4.1. 



Parameters 


rs value 


t value 


Parameters 


rs value 


t value 


DM^ z 


0.68955 


4.66 


dDM, z 


0.74017 


5.39 


DM, Pi5i 


0.68821 


4.65 


dDM, Pi5i 


0.54462 


3.18 




-0.31145 


-1.61 


dDM, Bsource 


-0.31009 


-1.60 


DM, aRM^ 


0.45257 


2.33 


dDM, (TnM. 


0.55929 


3.09 


DM, dDM, 


0.89197 


9.67 


dDM, dRM, 


0.40514 


2.03 



Table 4.4: Spearman rank values for all shifted depolarisation parameters. 



Figure 4.2(a) shows that the higher redshift sources display, on average, a larger de- 
gree of depolarisation and also a much larger spread of depolarisation at any given 
redshift compared to their low redshift counterparts. There is also very little dif- 
ference between the two high redshift samples indicating that the difference in the 
radio-luminosity of the two samples is unimportant, see Figure 4.2(b). Interestingly 
Morris & Tabara (1973) found depolarisation to correlate with radio-luminosity 
whereas Kronberg et al. (1972) found the correlation was predominantly with red- 
shift. My results support the findings of both groups. 

The trend with redshift is more pronounced when the difference in the depolar- 
isation over the source, dDM,, is considered (Tables 4.4 and 4.5(bottom)). Using 
Partial Spearman Rank I find 

rzdDM,,Pi5i = 0-60, D = 3.25 

rzPi5i,dDM, = 0.53, D = 2.77 

rp,,,dDM.,z =0.05, D = 0.23 

This shows a definite correlation between redshift and dDM, at any given radio 
luminosity, but the corresponding Pisi-dDM, correlation is considerably weaker at 
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Parameter 


1 
i 


z 


Q 
O 


/I 

4 


z 


u.ooo ( 


noise; 

U.ZioO 


n nnm 

-U.UUUi 


n 7son 

-U. ( o/U 


P 


0.4699 


0.6727 


-0.1909 


0.5387 


D 


-0.4397 


0.6135 


0.6368 


-0.1570 


DM, 


0.4951 


-0.3511 


0.7470 


0.2713 


Eigenvalue 


52.7% 


22.0% 


15.5% 


9.8 % 


z 


0.5957 


0.1851 


0.0132 


0.7815 


P 


0.4760 


0.6755 


-0.2183 


-0.5192 


D 


-0.4281 


0.6314 


0.6249 


0.1662 


dDM, 


0.4851 


-0.3329 


0.7495 


-0.3036 


Eigenvalue 


52.0% 


21.6% 


16.9% 


9.5 % 



Table 4.5: Eigenvectors and Eigenvalues for depolarisation 



a given redsliift. 

It is not surprising that DM, and dDM, are highly correlated, with a significance 
greater than 99.9% (Figure 4.3). A source with a high average DM, will be located 
in a region with a dense Faraday medium, which also causes a significant difference 
between the lobes and hence the correlation. The Laing Garrington effect discussed 
in section 3.3.1 could also cause the observed correlation between DM, and dDM,, 
assuming that the sources are not angled fiat to the line-of-sight. 
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Figure 4.3: DM, against the difference in depolarisation, dDM,. dDM, is the differ- 
ence in the average depolarisation of each lobe. The depolarisation has been shifted 
to a common redshift, z = 1. Symbols are as in Figure 4.1. 

Table 4.4 shows a very weak DM,-D<jource anti-correlation. This detection is 
marginal and as the PGA results in Table 4.5 (top) demonstrate this anti-correlation 
is not very significant. Interestingly this trend has been found by Strom (1973); 
Strom & Jagers (1988); Pedelty et al. (1989); Best et al. (1999) and Ishwara-Ghandra 
et al. (1998), but was not found by Dennett-Thorpe (1996) [DT96]. The observations 
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by Strom (1973), Strom & Jagers (1988) and Ishwara- Chandra et ah (1998) are at a 
lower resolution than my measurements and thus their depolarisation measurements 
may be affected by this resolution difference. The observations by Best et al. (1999) 
are at a higher resolution, but also at higher frequencies (4.8 GHz to 8.4 GHz) and 
so I would expect any depolarisation trends to be stronger in this sample. Pedelty 
et al. (1989) (P89) however observe at the same frequency range and resolution as my 
samples and they find an anti-correlation with size. DT96 also observes at the same 
frequency and resolution as P89, but also does not find any trend with size. Even 
when my results from sample B, which are taken at the most similar redshift and 
radio-luminosity range to those of P89, are analysed separately, I still find no trend 
with size. There is no difference in the selection of the P89 sample and my sample 
B except that P89 chose their sample from sources which had strong emission lines. 
By preferentially selecting sources with strong emission lines P89 has chosen sources 
that would have strong depolarisation asymmetries and hence large depolarisation 
measurements overall. This could strengthen any D^ource-DM correlation found. 
However, there does not seem to be any convincing physical explanation why one 
set of sources should show such a strong size-depolarisation anti-correlation and 
another set of sources should show no trend. 

It is usually assumed that radio sources are located in stratified atmospheres. 
Therefore a small source will be embedded in denser gas which acts as a more 
efficient Faraday screen. As the source expands the lobes have a higher probability 
of extending beyond the denser inner atmosphere, thus reducing the amount of 
depolarisation observed. In section 4.2 it was shown that there is no significant 
z-Dgource anti-correlation in my sample. Therefore a lack of any DM^-D^o^^ce anti- 
correlation suggests that there is no significant difference in the source environments 
at any redshift, of my samples, as the sources become larger. This may be evidence 
that on scales of up to a Mpc the environment is relatively homogeneous and does 
not show evidence of stratification. The lack of any D^o^^ce-dDM^ could simply be 
caused by the Laing-Garrington effect (see section 3.3.1). 



4.3.3 Rotation measure 



The Spearman Rank results in Table 4.6 show that there are no strong correlations 
of the fundamental parameters with rotation measure, RM. In every case the corre- 
lation is below 95% significant. The PGA results in Table 4.7(top) also demonstrate 
that the RM of a source does not depend on the size, redshift or the radio luminos- 
ity of the source. This is consistent with the idea that the observed RM is due to 
Galactic variations and not to a local Faraday screen (see section 3.3) and would 
not be expected to vary in response to any change in the source properties. 

The difference in the rotation measure, dRM^, between the two lobes of a given 
source and the rms variation in the rotation measure, auM^ , are local to the source 
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Parameters 


IS value 


t value 


Parameters 


rs value 


t value 


RMa z 


0.16996 


0.79 




0.67885 


4.24 


RM, Pi5l 


0.27470 


1.31 


O'RM^ Pl51 


0.34881 


1.71 


RM(j ^source 


-0.05040 


-0.23 


O'RM^ ^source 


-0.30435 


-1.46 


dRM, z 


0.50000 


2.65 


Grm, dDM, 


0.55929 


3.09 


dRM, Pi5i 


0.22530 


1.06 


aRM, DM, 


0.45257 


2.33 


ciR-lVI^ ^source 


-0.22332 


-1.05 


aRM, dRM, 


0.68874 


4.35 


dRM, dDM, 


0.40514 


2.03 









Table 4.6: Spearman rank values for all rotation measure parameters. The subscript 
a indicates that the RM has not been shifted to the sources' frame of reference. 



and have been corrected to the sources' frame of reference. There is a dRM,-z 
correlation present in the Spearman Rank results, (see Table 4.6), but it is weaker in 
the PCA results where the correlation is strongly reversed in the 3rd eigenvector, see 
Table 4.7 (middle). Figure 4.4 shows that in principal I should find a strong dRM^-z 
correlation. However, the presence of the low redshift - high dRM, source (3C 299) 
weakens this trend. 3C 299 has a large dRM, produced by a lack of polarisation 
information in one lobe. Removing this source strengthens the correlation giving 
TdHM^ z = 0.632 which corresponds to a significance of 99.98%. 




500 



dRM, 



Figure 4.4: Redshift against difference in rotation measure, dRM, across the lobes 
of a source. dRM, has been shifted to the sources' frame of reference. Symbols are 
as in Figure 4.1. 

As previously noted, (Trm^ measures variations of the Faraday screen on scales 
smaller than dRM,, which may still be influenced by the Galactic Faraday screen. 
However, it is interesting to note that dRM, correlates strongly with ctrm^, see 
Figure 4.5. So the two parameters may sample the same, local Faraday screen 
even if dRM, may still contain some Galactic contribution, ctrm^ shows a strong 
correlation with redshift, > 99.9% significance, and a weak anti-correlation with 

■Dsoiirce* The 

c rm^^"^ correlation also strengthens with the removal of 3C 299, be- 
coming TaRM^^z = 0.767 which corresponds to a significance exceeding 99.99%. These 
correlations are more obvious when the PCA results are considered, see Table 4.7 
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Parameter 


1 

1 


o 

L 


Q 
O 


4 


z 


u.ooyo 


n OQon 
U.zoyU 


A n/ioi 
-u.U4y i 


U. i Zoo 


r 


U.59d8 


n /I An /I 

0.4094 


-0.1678 


A /^/^ A /I 

-0.6694 


U 


-U.3781 


0.4402 


A OAO A 

-0.8034 


A 1 O O C 

0.1335 


KM 




-0.7625 


A r A o 

-0.5693 


A A r /I A 

-0.0540 


Eigenvalue 


4o.o/o 


^4.i /o 


Zi.O/o 


1 n ro/ 
iU.o/o 


z 




-U.UZUi 




r\ 7AA7 
U. ( 44 ( 


P 


0.5820 


-0.1197 


0.4874 


-0.6399 


D 


-0.3224 


-0.8641 


0.3594 


0.1423 


dRM^ 


0.4017 


-0.4885 


-0.7644 


-0.1255 


Eigenvalue 


45.9% 


23.3% 


20.3% 


10.5 % 


z 


0.6215 


0.1765 


0.0623 


0.7607 


P 


0.4923 


0.4210 


-0.6151 


-0.4495 


D 


-0.3212 


0.8876 


0.3288 


0.0296 




0.5178 


-0.0615 


0.7139 


-0.4673 


Eigenvalue 


51.3% 


22.4% 


18.6% 


7.6% 



Table 4.7: Eigenvectors and Eigenvalues for rotation measure. 



(bottom). The first eigenvector contains 51.3% of the variation in the data and 
shows a strong 'l-cfrMz correlation and a weaker ^source-^^RM^ anti-correlation. The 
D source- o'ijAf^ anticorrclation reverses in the third eigenvector indicating that it is a 
considerably weaker trend compared to the z-cfrm^ correlation which only reverses 
in the last eigenvector. The Partial Spearman Rank results below confirm that there 
are no significant correlations with radio-luminosity for both dRM;^ and aRM^, 

r^dRM.P,,, = 0.49, D = 2.40 

rzPi5i,dRM, = 0.69, D = 3.79 

rp,5idi?M.,. - -0.19, 1^ = 0.86 



rzanM^Pwi = 0.65, D = 3.47 

^^Pi5i,<TflM, = 0.67, D = 3.63 

^Pi5i.HM.,.= -0-24, D=1.09 



In a survey of 27 high redshift sources Pentericci et al. (2000) [P02] found that the 
Faraday rotation was independent of size or radio-luminosity but they also found 
that the number of sources with high levels of Faraday rotation increased with red- 
shift. The results from P02 use the RM of a source in the sources' frame of reference. 
As noted in section 3.2 the RM from my sources is taken to be dominated by contri- 
butions from the Galaxy thus 1 use ctrm^ and dRM^ as indicators of the environment 
instead of RM. ctrm^, dRM^ and RM sample the same medium and both sets of 
results find a trend with redshift. The P02 trend indicates that the strength of 
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dRM, 



100 



Figure 4.5: The rms variation in the rotation measure, ct/jm^ against the difference in 
rotation measure across the lobes of a source. Symbols are as in Figure 4.1. 



the magnetic field or the density of the environment (or both) is increasing with 
redshift but by using only RM they are unable to determine how the disorder of 
the environment changes with redshift. Conversely by using only (Trm^ and dRM^ 
I have no information on the density of the environment or the magnetic field but I 
find that the disorder in the environment increases with redshift. 

The Spearman Rank results in Table 4.6 show that there is a strong (T/ja/^-DM^ 
correlation, indicating that (TrM:, is sampling the same medium as the depolarisation 
measurements and thus either parameter can be used as a test of the Faraday screen. 



4.3.4 Summary 



Depolarisation, DM^, and the differential rotation measure properties, dRM^ and 
ajiM^, correlate with redshift, indicating that there is a change in the source en- 
vironments as redshift increases. Only depolarisation shows a trend with radio- 
luminosity. By analysing the rotation measure down to its rms variations {ctrmJ 
and the depolarisation measure, I am probing smaller and smaller variations in the 
Faraday screen. I find that there is an increasingly strong correlation with redshift 
over radio-luminosity which suggests that the environment has no direct link to the 
radio-luminosity of a source. The fact that there are no spectral index correlations 
found also corroborates this view. 
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4.4 Source asymmetries 



The previous sections have allowed me to analyse the bulk trends of the source 

properties with the fundamental properties, but also with themselves. However, it 
is also interesting to look at the asymmetries of the sources and to determine their 
underlying physical cause. 

The angle at which a source is orientated to the line-of-sight, 9, can affect the 
projected lobe length. Sources orientated at a small angle will have lobes that 
appear more asymmetrical than sources at a large angle, assuming there are no 
environmental differences between the lobes. The orientation of a source will also 
affect the observed depolarisation (Garrington & Conway, 1991), the spectral index 
(Liu & Pooley, 1991b) and the rotation measure. In a simple orientation model 
emission from the lobe pointing away from the observer will have a longer path 
length through the local Faraday screen and thus any depolarisation measurements 
of this lobe would be larger than the lobe pointing towards the observer. This is 
known as the Laing- Garrington effect, see section 3.3.1. Liu & Pooley (1991b) also 
find that the least depolarised lobe has a flatter spectrum. 

Beaming can also cause asymmetries in a source, independent of path length 
through the Faraday screen. Assuming a source is angled at 6 to the line-of-sight, 
then the ratio of lobe lengths, D1/D2, is given by (Longair & Riley, 1979) 

Di 1 + Vo/c cos6 
D2 1 — Vo/c cos6 

where Vo is the velocity at which the hotspot (Longair & Riley, 1979) is moving away 
from the nucleus and c is the speed of light. As the hotspot becomes increasingly 
more beamed (i.e. ^ — > 0) the lobe ratio increases. The hotspot will begin to 
dominate the flux and spectral index in the beamed lobe. Thus as beaming becomes 
more dominant in a source, the beamed lobe will become brighter, longer and its 
spectrum flatter than the receding lobe. 

Half of my sources do not contain any detection of a core at either 4.8 GHz or 
1.4 GHz. This was a problem as I needed to deflne where the core was located to 

determine an accurate estimation of the lobe length. In the cases were there was no 
core detection 1 used previously published maps to determine where a core is likely 
to be located. However in several cases, most notably 3C 16, there was no core to be 
found. In these cases I estimated the location of the core position by using the most 
likely core position from the hterature (Leahy & Per ley, 1991). This may seem a 
rather drastic approach, but as the maps are not of very high resolution, the values 
calculated for the lobe length can only be an approximation to the true length. The 
lobe volume is calculated by assuming cylindrical symmetry and using the measured 
lobe length, Diobe, but also the parameter R, which determines the ratio of the lobe 
length to the width at half the lobe length. 
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V,o,e = ^ (4.6) 

In most cases R > 1 but in extreme cases like the southern lobe of 7C 1745+642 
(see Figure ??) R is less than unity, which causes the large volume ratio (> 10) for 
this source apparent in Figure 4.18(a). 

In the following sections the subscripts b and / are used to denote the value from 
the brighter lobe and the fainter lobe of an individual source, respectively. 

In the previous section I used the corrected z = 1 depolarisation in my analy- 
sis of the trend of depolarisation with the fundamental parameters. However, by 
correcting the depolarisation to a common redshift I am assuming a specific under- 
lying structure to the Faraday screen which may mask any underlying asymmetries 
between lobes. In previous studies (e.g. Pedelty et al., 1989; Liu & Pooley, 1991a; 
Gregory & Condon, 1991, etc.) the observed depolarisation is used as an asymme- 
try indicator. To ensure my results are comparable to those published I use only 
the observed depolarisation values given in Tables 3.1 to 3.3 in my analysis of any 
depolarisation asymmetries. 



4.4.1 Flux asymmetries 
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Figure 4.6: Plot of the flux of the brighter lobe against the fainter lobe at 4.8 GHz. 
Symbols are as in Figure 4.1. 

Figures 4.6 to 4.8 use the 4.8 GHz flux data as a test for flux asymmetry over 
the source. Tables 3.1 to 3.3 show that the 1.4 GHz flux is much larger compared 
to the 4.8 GHz flux for each source. Although this suggests that the 1.4 GHz flux 
data would be a better indicator of any source asymmetries than the 4.8 GHz flux 
data, since fainter lobes such as the northern lobe of 7C 1745+642 are much more 
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prominent, but the noise level in each source is also much higher. To ensure that 
my analysis was not unduly affected by high noise levels I chose to use the 4.8 GHz 
ffux instead of the 1.4 GHz flux, in my analysis of the asymmetries. 



Figure 4.6 demonstrates that in all my sources there is a tight correlation between 
the flux of the brightest lobe when compared to the flux of the fainter lobe. In several 
cases the fainter flux is much smaller than the brighter flux. 3 of these sources are 
quasars, 3C 351, 3C 68.1 and 3C 16 but interestingly 3C 299 and 3C 268.1 also 
display this atypical behaviour. McCarthy et al. (1995) noted that 3C 299 showed 
a strong asymmetry and was more similar to a high redshift {z >1) radio galaxy 
than a low redshift galaxy. Harvanek & Stocke (2002) also noted that this source 
has lobes of very different lengths. I find 3C 299 to have a ctjim^ comparable to the 
high redshift samples. 
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Figure 4.7: Plot of the flux ratio at 4.8 GHz against the lobe ratio, (a - left) with 
all sources and (b - right) zoomed in on the more 'normal' sources. Symbols are as 
in Figure 4.1. 



Figure 4.7(a) shows that 3C 299, 3C 268.1 and 3C 16 do show larger (>1) lobe 
ratios but 3C 351 and 3C 68.1 do not. However 3C 351's position in this plot is 
uncertain. The estimation of the lobe length requires an estimation of where the 
lobe terminates but in the case of 3C 351 's northern lobe it extends far to the west 
of the hotspots whereas the southern lobe is much more well defined (see Figure ??). 
Thus the estimation of the length of the northern lobe may be incorrect and hence 
the calculated lobe ratio may be smaller than it really is. 3C 68.1 and 3C 16 have no 
detected cores either by me or previously published data^. This means that the lobe 
length estimation in these sources may also be flawed. To further aggravate this 
problem both of these sources have very faint lobes (see Figures ?? & ??) and so it 
is hard to determine where exactly the lobe terminates. However, it is worth noting 
that when these 5 sources are removed from the analysis (Figure 4.7(b)) there is no 
significant trend for the brighter flux ratio sources to have a larger lobe ratio. 

■^Bridle et al. (1994) do marginally detect a core in deeper observations of 3C 68.1. 
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Figure 4.8: (a - left) Plot of the flux ratio at 4.8 GHz against redsliift, (b - right) 
and against radio-luminosity. Symbols are as in Figure 4.1. 



Figures 4.8(a) and 4.8(b) show that the trends noted are intrinsic to the sources 
and not due to any redshift and radio-luminosity differences in the samples. It is 
worth noting that sources in the high-redshift, low radio-luminosity sample, sample 
A ('x'), always seem to have a smaller flux ratio than sources in the other two 
samples. 



4.4.2 Asymmetries in spectral index 
4.4.2.1 Liu-Pooley effect 

The Liu-Pooley (LP) effect predicts that the lobe with the flatter spectrum in a 
source is also the least depolarised lobe (Liu & Pooley, 1991a; Ishwara-Chandra 
et al., 2001). Liu & Pooley (1991a) found that 12 out of their 13 sources demon- 
strated this effect. As noted in section 4.4, to allow a fair comparison between my 
results and those of Liu & Pooley (1991a) I use only the observed depolarisation in 
the following analysis. 

As Figure 4.9 shows I find no significant Liu-Pooley (LP) effect in my sample. As 
another comparison I used the rms variations in the rotation measure in the sources' 
reference frame, (Jrm^, instead of depolarisation measure. Clearly Figure 4.10 shows 
that with ajiM^ I fi^id no significant evidence of any LP effect. Even by plotting 
the difference in depolarisation or the difference in a^M^ against the spectral index 
difference, insuring that the difference in DM or ctrm^ were always positive, there 
was no significant LP effect. In all cases I find only a maximum of 14 sources out 
of the 26 sources (12 out of 23 in the case of ctrmJ to show definite signs of the LP 
effect. 
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Figure 4.9: The difference in the spectral index of the bright lobe to the faint lobe 
against the difference in the depolarisation of the bright lobe to the faint lobe. 
Depolarisation is taken to be the observed (unshifted) value. Symbols are as in 
Figure 4.1. 
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Figure 4.10: The difference in the spectral index of the bright lobe to the faint lobe 
against the difference in the ajiM^ of the bright lobe to the faint lobe. (Trm^ is the 
value local to the source. Symbols are as in Figure 4.1. 



Ishwar a- Chandra et al. (2001) [ICO] have shown that the LP effect is stronger for 
smaller sources and shows no dependence on redshift. They also find that the LP 
effect is not affected by the choice of quasars or radio galaxies in a sample and 
is equally significant in both species. They also find that for their large sample 
(comprised of the Liu & Pooley (1991a); Pedelty et al. (1989); Garrington et al. 
(1991) samples and the Molongo Reference Catalogue/1 Jy sample) only 58% of 
radio galaxies and 59% of quasars show the LP effect which is consistent with my 
results. The inherent difference between sources chosen in samples as well as orien- 
tation effects has been suggested by ICO as an explanation to the number of sources 
showing the LP effect, in any given sample. 
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4.4.2.2 Effect of changes in the lobe length 
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Figure 4.11: The ratio of spectral indices against the lobe ratio. Symbols are as in 
Figure 4.1. 

Figure 4.11 shows that the ratio of spectral indices of the sources' lobes is insen- 
sitive to changes in the lobe volume. A source with a large lobe volume ratio (> 
1) is equally likely to have a low spectral index ratio as a high spectral index ratio. 
This suggests that the environment, which affects how large a volume a source can 
grow to, has no direct effect on the spectral index of a lobe. Thus spectral index 
should not be used as a tracer for the environment. 
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Figure 4.12: The spectral index ratio against the flux ratio, (a - left) all sources and 
(b - right) zoomed in on sources with a flux ratio < 5. Symbols are as in Figure 4.1. 



Figure 4.12(a) shows that in general, a source with a higher flux ratio has a smaller 
spectral index ratio. This means that the brighter lobe becomes progressively flatter 
(i.e. ab/df — >■ 0) as the difference in the flux ratio's increases. This is consistent 
with the theory that beaming is present in these sources. 3C 16 and 3C 299 are 
notable exceptions to this trend. 3C 299 has already been discussed above and 
possess exceptional properties when compared to other radio galaxies. 3C 16 is a 
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quasar and I would expect this source to have shown more evidence for beaming but 
it has almost similar spectral indices in both lobes. Harvanck & Hardcastle (1998) 
found evidence that 3C 16 maybe a radio source that is just beginning to restart 
and hence the fainter lobe will not have a well defined hotspot, therefore beaming 
is currently not important in this source. Figure 4.12(b) shows that when the more 
extreme sources are removed the anti-correlation is even more evident. 
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Figure 4.13: The ratio of spectral indices against the lobe ratio. Symbols are as in 
Figure 4.1. 

Figure 4.13 shows that as the lobe length ratio increases the spectral index ratio 
decreases. This means that in the more length asymmetric sources the larger lobe 
has the flatter spectral index. Once again 3C 299 again contradicts this result. It has 
a high ratio of spectral indices and a high lobe length ratio. The trend between lobe 
ratio and the spectral index ratio adds further evidence to the presence of beaming 
in many of my sources. 

Although Figures 4.12 and 4.13 show evidence for beaming, another possibility 
would be that the observed asymmetries are caused by a difference in the density of 
the environment between both lobes. A denser environment would mean that the 
lobe has a higher minimum magnetic field. Thus any asymmetries in the minimum 
energy magnetic field would indicate that it is the environment and not beaming 
that is causing the observed asymmetries with spectral index. 

4.4.2.3 Minimum energy of the magnetic field 

A simple method to estimate the magnetic field of a radio source is to use the 
minimum energy argument (e.g Longair, 1994). For any given radio source of a 
given volume and radio-luminosity the minimum energy magnetic field argument 
calculates the minimum energy that must be present in the emitting volume. The 
minimum energy argument assumes that VFmag = |M/particies, i-e. a source has ap- 
proximately equal energy in its relativistic particles as it has in its magnetic field. 
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The minimum energy approach assumes a perfectly uniform magnetic field through- 
out the lobe. There are no energy losses from the electrons. Therefore, the energy 
distribution of the relativistic electrons is a perfect power-law throughout the life- 
time of the source. This is a simphstic model but it can be used to model the 
magnetic field, B^m, in the lobe with few assumptions required. For full details see 
Longair (1994) which gives. 



B„ 



where 



i_ 



V 262 



ml K 



(Tesla) 




(4.7) 
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and X = \ — la. The volume of a lobe, V, is assumed to be cylindrical with 
V = ^Dfobe/ before, ui = 10 MHz and 1^2 = 100 GHz are taken to be the lower 
and upper limits of the synchrotron spectrum and the luminosity, L^, is calculated 
using the 4.8 GHz flux of each lobe. 
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Figure 4.14: (a - left) The minimum energy magnetic field ratio against the fiux 
ratio and (b - right) against the spectral index ratio. Symbols are as in Figure 4.1. 

As Figure 4.14(a) demonstrates, there is no relationship between the minimum 
energy magnetic field ratio and the fiux ratio. This indicates that there is no 
significant environmental difference causing the fiux differences between two lobes 
of a source. 



Correspondingly Figure 4.14(b) shows that there is a weak correlation in the 
minimum energy magnetic field ratio when compared with the spectral index ratio. 
A beamed, fiat spectrum hotspot will have a stronger magnetic field which could 
explain the weak trend in Figure 4.14(b). However, the asymmetries in the minimum 
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energy magnetic field could also be caused by differences in the density of the 
environmental across the lobes. Figure 4.15 shows that there is no trend with 
redshift and so the minimum energy magnetic field- a ratio is intrinsic to the local 
source environment and is not affected by the redshift of a source. 
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Figure 4.15: The minimum energy magnetic field ratio against the redshift of the 
source. Symbols are as in Figure 4.1. 
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Figure 4.16: The minimum energy magnetic field ratio against (a - left) the depo- 
larisation ratio and (b - right) ctrm^ ratio. Symbols are as in Figure 4.1. 



It is interesting to note that the asymmetries in depolarisation and ojiM:, are 
not related to asymmetries in the minimum energy magnetic field (see Figure 4.16). 
Depolarisation and orm^ are both thought to be observational indicators of density^ , 
in the source environment but neither show any relationship with the minimum 

^Rotation measure is a direct function of density and the magnetic field of a source (section 3.2) 
and since (Jrm^ and the shifted depolarisation have been shown to be correlated (section 4.3.3) 
then both can be said to be tracers of the density of the environment. 
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energy magnetic field ratio. This suggests that the environment has no direct impact 
on either depolarisation or (Trm^ asymmetries. 



4.4.3 Asymmetries in depolarisation 

Unhke the flux asymmetry (Figure 4.6) the depolarisation shows no obvious correla- 
tion between the depolarisation of the brighter and fainter lobes, see Figure 4.17. It 
is worth noting that in most cases the depolarisation of the brighter lobe is smaller 
than the depolarisation of the fainter lobe. This is again indicative of beaming in 
these sources. In the case of beaming, the fainter lobe is assumed to be pointing 
away from the observer. The radiation of the fainter lobe will then have a larger 
path length to traverse through the local Faraday screen compared to the radiation 
of the brighter lobe and will appear more depolarised. 

ro I— I , . , , . , , . , 1 . , , . 1 
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Figure 4.17: Depolarisation of the brighter lobe against the fainter lobe. Symbols 
are as in Figure 4.1. 



The flux asymmetry is independent of any depolarisation asymmetry as 3C 299, 
3C 16, 3C 68.1, 3C 351 and 3C 268.1 (the high flux ratio sources) show similar 
spreads in their depolarisation ratios when compared to all the other sources. Figure 
4.18(a) at flrst sight indicates that the volume of the lobe also plays no part in the 
magnitude of the depolarisation ratio. However, on closer inspection it is worth 
noting that in 16 out of the 26 sources, the lobe with the larger volume is the more 
depolarised. In fact, only 4C 14.27, 3C 68.1, 3C 351, 6C 1256+36 and 7C 1745+642 
vary signiflcantly from this trend. 

By plotting the depolarisation ratio against the lobe volume ratio in Figure 
4.18(a) I assume that all the lobes in the three samples can be described by a 
cylinder. This is not a good approximation for extremely distorted sources like 3C 
351 and 7C 1745+642. As Figure 4.18(b) shows the trend still exists when the dif- 
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Figure 4.18: (a - left) Depolarisation ratio against the lobe volume ratio, (b - right) 
Difference in depolarisation against the difference in size of the lobes. Plotted so 
that Diobei-^iobe2> 1 cases. Symbols are as in Figure 4.1. 



ference in the lobe lengths is compared with the difference in the depolarisations of 
each lobe. In this case the trend is much more pronounced with only 6C 0943+39 
and 7C 1745+642 showing significant disagreement with the findings. Both of these 
sources show the same trends in (a) and (b). This disagrees with the results of 
Pedelty et al. (1989) and Liu & Pooley (1991a), who find that the shorter lobe is 
the more depolarised lobe. However DT96 also finds this reverse correlation with 
only 3C 351 in common between her sample and my sample. The DT96 correlation 
was found to be strong in the quasars, but the depolarisation of the radio galaxies 
was too low to be able to determine if the correlation was significant. This may 
indicate that there is an inherent difference in the Pedelty et al. (1989) and Liu & 
Pooley (1991a) samples and the DT96 sample and mine. This may be due to the 
fact that the Pedelty et al. (1989) sample was selected from sources with strong 
emission lines and the Liu & Pooley (1991a) sample only contains sources smaller 
than 23". 



4.4.4 Asymmetries in auM, 



Figure 4.19 shows there is no relationship between urm^ of the brighter lobe when 
compared to the fainter lobe. I find that only half the sources show that the fainter 
lobe has the larger a^M^ ■ This is unlike the depolarisation result where in almost all 
cases the fainter lobe was more depolarised. However, the depolarisation is affected 
by changes in the path length through the Faraday screen but <Jrm^ is insensitive 
to these changes and so would not be expected to show any trend with the lobe 
asymmetries. 

As with depolarisation there is little trend with lobe volume (Figure 4.20(a)). 
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Figure 4.19: ctjim^ of the bright lobe against the fainter lobe. Symbols are as in 
Figure 4.1. 



This confirms the idea that the extent that a lobe can grow to (i.e the volume) has 
little effect on the measured parameters such as spectral index, depolarisation and 
aRM^- Most of the sources show that the larger a^M^ is in the longer lobe. In only 
3 out of the 23 sources is this not the case (6C 1129+37, 6C 0943+39 and 3C 280) 
which is similar to the trend noted with depolarisation. 
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Figure 4.20: (a - left) ctrm^ ratio against the lobe volume ratio, (b - right) Difference 
in aRM^ against the difference in size of the lobes. Plotted so that Diobei-^iobe2> 1 
in all cases. Symbols are as in Figure 4.1. 



4.4.5 Summary 

I found that there is no significant LP effect, but this lack of trend could be attributed 
to the larger size distribution of my sample compared to the sample of Liu & Pooley 
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(1991a). This dependence on size has also been noted by ICOl. From the spectral 
and depolarisation asymmetries 1 found evidence that many of the sources were 
beamed. The minimum energy magnetic field asymmetry results could be explained 
by both a change in the environment or by the presence of beaming. The minimum 
energy magnetic field model is not an accurate picture of the magnetic field of a 
source but does give a rough estimate of the size of the field. 

The depolarisation and a^M^ asymmetries are similar, confirming the assumption 
that they sample similar regions of the Faraday screen. In general, the depolarisation 
and (TuM;, are smaller in the shorter lobe and almost all the fainter lobes show larger 
depolarisations. The a^M^ variations are insensitive to which lobe is brighter. 



4.5 Conclusions 

A number of source properties are derived directly from the radio observations. 
These include the spectral index, the rotation measure and the depolarisation mea- 
sure. I find no correlation of spectral index with any of the fundamental parameters 
(redshift, radio- luminosity and radio-size). I can therefore not confirm previous 
findings of a correlation of spectral index with redshift (Athreya & Kapahi, 1999), 
a trend of spectral index with radio luminosity (Veron et al., 1972; Onuora, 1989) 
nor a trend of spectral index with size (Blundell et al., 1999) 

In the previous chapter I argued that the overall rotation measure for a given 
source is mainly Galactic in origin. This interpretation is consistent with the ab- 
sence of any correlation of rotation measure with the fundamental parameters. Vari- 
ations in the rotation measure are quantified in order of decreasing angular scales by 
dRM^. (difference of rotation measure between the two lobes of an individual source) , 
apiMz (rms variation of rotation measure) and DM^ (depolarisation measure). 1 find 
that the significance of a correlation with redshift increases from dRM^ to DM^. 
Therefore I can conclude that the variation of the rotation measure on small an- 
gular scales is caused by a Faraday screen local to the sources. Furthermore, these 
Faraday screens give rise to greater variation of the rotation measure for sources at 
higher redshift compared to their low redshift counterparts. 

I have also studied the asymmetries of these derived properties with the funda- 
mental properties, but I have also looked at the asymmetries in the 4.8 GHz flux 
data, S, and the minimum energy magnetic fleld, B^j„. By deflnition, Sfe/S/ > 1 but 
several sources, of which 3 are quasars, show a much larger flux ratio. I flnd evidence 
from spectral asymmetries that many of the sources show indications of beaming. 
However, the results from the Bmin asymmetries are inconclusive and could also 
explain the asymmetries by implying a more environmental based explanation. In 
general all asymmetries were found to be insensitive to changes in the source red- 
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shift, or low frequency radio luminosity; suggesting that the environments of the 
high redshift and low redshift samples are similar and any asymmetries caused are 
due to localised changes in the environment. 

This interpretation is confirmed in the analysis of the depolarisation and (Jrm^ 
asymmetries with B,„i,j,„. I suggest that it is not the environment that is causing 
the depolarisation and (Jjim^ asymmetries but is actually the degree of disorder in 
the sources that is causing a change. This fits with my findings that the shifted 
depolarisation, DM^, dRM^^ and (Trm^ show trends with redshift but no trend with 
the minimum energy magnetic field. Thus it is the degree of disorder in the magnetic 
field that is changing with redshift and not the density of the source environment. 



Chapter 5 



Modelling of the Faraday screen 



To sample the environment around radio sources I use the observed polarisation 
and the observed RM variations. These parameters are used to determine how the 
environment evolves with redshift, radio-luminosity and radio size*. However, the 
observed polarisation, the observed RM variations and hence the observed depo- 
larisation measure are also effected by source based variations in the structure of 
the Faraday screen. The structure of the Faraday screen combined with the limited 
spatial resolution of any radio observation lead to the observed depolarisation of the 
synchrotron emission. 

A uniform Faraday screen comprised of large (> 50 kpc) cells will have little 
effect on the underlying polarisation and RM structure. A cell is defined as a region 
in the Faraday screen where the magnetic field has almost uniform strength and 
orientation. As these cells shrink (i.e the cells size becomes smaller than the size of 
the telescope beam) and the Faraday screen becomes increasingly non-uniform the 
effects of the underlying structure can become more pronounced. The effect of the 
Faraday screen on a source was first investigated by Burn (1966). I use the Tribble 
(1991), hereafter T91, models as a modern alternative to Burn's theory. The T91 
models assume a specific structure to the RM variations and calculate the effect 
of the structure function on the underlying RM distribution and hence the overall 
polarisation of a source at any given wavelength. These models can be used to 
correct the observed depolarisation of my sample sources to a common redshift and 
the results can be compared with the Burn shifted depolarisation, DM^, in section 
3.3. If the depolarisation measure is found to be correlated with redshift and radio- 
luminosity using the T91, models then this will confirm the findings in section 4.3.2 
where DM^ is found to be a function of both redshift and radio-luminosity. 



*The relationship between the environment and the depolarisation measure and uhm^ are 
discussed in detail in the previous chapters. 
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5.1 Critique of Burn's model 



The polarisation, p, at a given point in a radio map can be defined as a function of 
position within the source, x^, and wavelength. A, 



pixl, A) = / Wixl - f)e2*^^(*)-^'(i2f (5.1) 



where W is the restoring beam of the telescope and RM^ is the intrinsic rotation 
measure of the Faraday screen in front of the source. 

By assuming that the RM can be defined as a Gaussian distribution with a dis- 
persion a (T91), equation 5.1 reduces to < p(A) >= e~'^'^ ^ which is more commonly 
known as Burn's law for a foreground depolarisation screen (Burn, 1966). 

Burn noted that this law only applies when the scale of the Faraday screen fiuc- 
tuations are much smaller than the physical size of the source. Observations of 
Cygnus A by CariUi & Barthel (1996) found that their RM fluctuations were co- 
herent over scales of order 10 kpc. This scale length is comparable to the results 
found by Pedelty et al. (1989), who find scale lengths between 10-50 kpc, using a 
sample of high redshift, high radio-luminosity sources. These coherence lengths are 
not significantly smaller than the source size, violating one of the main assumptions 
of Burns' theory. 

This is not the only problem associated with Burn's theory. In his analysis of the 
polarisation properties of a source. Burn used the vector < p{X) > as a test for the 
polarisation which falls off rapidly to zero as the direction of the polarisation field 
randomises. Burn also predicts that the observed polarisation decreases steeply as a 
function of wavelength, which is very rarely observed. Table 5.1 demonstrates than 
even in a small sample of radio galaxies the decrease of the polarisation fraction 
with wavelength varies widely from the Burn relationship, p ~ A^'^. This confirms 
the assumption that by using Burn's theory to shift the depolarisation to a common 
redshift in section 3.4.2 I am almost certainly over-estimating the depolarisation. 
Any observations of the polarisation of radio emission measures the length of the 
polarisation vector, thus < |p(A)p > is a better indicator of the polarisation com- 
pared to Burn's < p{\) >. In the next section 1 explain the main points of the T91 
models which use < |p(A)p > to test the Faraday medium around a source. For full 
details see T91. 



^See section 3.2 for a full explanation of the rotation measure of a source. 
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5.2 The Tribble models 



The polarisation distribution as a function of wavelength, independent of any source 
model, is given by, 

< |p(A)|^>= J Ris)U^£s (5.2) 

where R{s) is determined by the size of the telescope beam W on the sky and 
Cp{s) = Cf{s, A) Ce{s) Ce(s). Ce{s) determines how the emissivity of neighbouring re- 
gions interact and becomes important when a radio source has prominent filamentary 
structures (Fomalont ct al., 1989). Ce{s) determines how the intrinsic polarisation 
angle of neighbouring regions interact. For a uniform source both the emissivity, 
ie{s), and the angle, ^6i(s), are unity, therefore ip{s) depends only on the distribution 
of the external Faraday screen, ^j(s. A). 

C,f{s, A) is the critical component in the determination of the polarisation as this 
relates directly to the structure function of the RM fluctuations, D{s). 

e/(.-',A)=e-^^^^(^-^ (5.3) 

where D{s) determines the RM distribution in a source and how neighbouring re- 
gions of local RM are correlated, i.e. D{s) =< [RM{x + s) - RM{x)]^ >. The exact 
form of the structure function depends on whether or not the RM field is assumed 
to be dominated by Quadratic fluctuations, Gaussian auto-correlation fluctuations 
or Power-law fluctuations. 



The exact form of R{s) is dependent on the structure of the telescope beam, 

R{s)^ J W{x)W{x + s)d'^x (5.4) 

assuming a Gaussian beam this reduces to, 

m - eH.IV.= (5.5) 

where t is the physical size of the telescope beam. The Tribble models are calculated 
in the restframe of the source. The size of the RM fluctuations, So, t and the 
observing wavelength, A must all be calculated in the source frame of reference. 



5.2.1 Quadratic model 

For the quadratic model the structure function is given by, 

D(s)= 2a\s/sof s<So (5.6) 

D{s) = 2(7^ s > So 
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where So is the size of the RM fluctuations within the telescope beam. t. The 
number of independent RM cells within a telescope beam is given by = (t/so)^- 
A completely unresolved source has So ^ t whereas a completely resolved source 
has So ^ t. 



The structure function can be substituted into equation 5.2 to give an expression 
for the polarisation, 



1-, — ' 



where a = sl/2t^ and (5 = 4cr^ A"^ 

However, to calculate the observed polarisation at any given wavelength the in- 
trinsic a of the source must also be known and this is affected by the choice of the 
structure function. 



For the quadratic model, 

<^observed = 4m[1 " (1 " (5.8) 



where (Jobserved is taken from the observations of the RM variations but the exact 
value of So is not known. In theory it should be possible to estimate the value of So 
from high resolution RM maps as So should be comparable to the scale length over 
which the RM is roughly constant. This would require higher resolution maps than 
those presented here, thus it is not possible to directly estimate the value of So from 
my observations. In section 5.3 I discuss the effects of varying So on the observed 
polarisation and the RM variations. 



5.2.2 Gaussian auto-correlation model 

The Gaussian auto-correlation structure function is given by, 

D{s) = 2a^[l - e-''/''] (5.9) 
Again substituting the structure function into equation 5.2 gives, 

< |p (A ) |2 >= ae-^ f e^'dt (5.10) 
Jo 

and the dispersion in the RM is given by. 
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5.2.3 Power-law model 



The power law model is slightly more complicated and depends on the choice of m, 
the power- law exponent. The structure function is given by, 

D{s) = 2(7' [1 - (1 + 2s7ms2)]-"^/2 (5.12) 

which using equation 5.2 gives, 

/oo 
g-ma./2g/3.-'»/2^^_ (5.13) 

This form of the structure function results in the relation p ~ X~^^"^. 

The dispersion of RM in this model is much more complex than either of the 
previous models and is given by, 

<^observed = 4m[W2]"^/' r(l - m/2, ma/2) e^^^/' (5.14) 

/•oo 

where T{a,x) = / e^H^'^dt. 



In all 3 models when the source is completely resolved 0"^, 



observed "^RM 



2 



5.2.3.1 Power law exponent 



T91 uses m = 2 in his analysis, but a sample of radio galaxies may show a range of 
exponents. The only way to determine how polarisation decreases with wavelength 
is to take a sample of sources with at least 3 sets of polarisation information widely 
separated in frequency, but also at high enough resolution that the polarisation 
measurements are not affected by beam dcpolarisation. The sample must also be 
observed at the same resolution at all three frequencies. Although there is a large 
volume of polarisation information in the archives, it proved to be difficult to find a 
large enough sample of sources with these requirements. Table 5.1 shows the best fit 
to p ~ \-^/^ using sources taken from Akujor & Garrington (1995) with observations 
at 1.4, 4.8 and 8.4 GHz. In general I find that the polarisation decreases with l^m< 
4. The north lobe of 3C 266 shows a much steeper decrease in its polarisation, but 
as Figure 5.3 demonstrates the effects of changing m are small. I use m — 2 and 
m = 4 in the following sections as limiting cases. 



5.3 Fitting the T91 models to the observations 

The T91 models have 2 free parameters, So and a^M and two parameters taken 
from observations, t and the observed variations in the RM, (Johservedi plus a third 
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Source 


Component 


— 4/m 


3C 67 


T 


1.03 


3C 173 


T 


1.79 


3C 213.1 


1 


1.72 


3C 266 


N 


0.26 


3C 266 


S 


1.72 


3C 299 


s 


2.34 


3C 305.1 


T 


3.17 


3C 309.1 


T 


1.29 


3C 346 


E 


3.94 


3C 380 


E 


3.16 


3C 380 


W 


1.62 


3C 455 


T 


2.56 



Table 5.1: Calculating the range of exponents ,m, seen in sources observed by Akujor 
& Garrington (1995) 

parameter, m, for the power-law model. There are three basic assumptions that can 
be made about the values of So and auM' 

i. ) The observations of cXobserved and z can be fitted assuming a constant So and 

a constant aRM with respect to redshift. If So is constant at all redshifts 
then this assumes that cells in the Faraday screen are always the same size 
independent of the source morphology or location and will have the same effect 
on the underlying RM distribution independent of redshift. By also setting 
CRM to be constant with redshift I am assuming that all sources have the same 
underlying RM distribution. This is a rather naive picture as Table 4.6 shows 
a strong evolution in the observed variations in the RM with redshift. 

ii. ) A more realistic assumption is that the observations of cTobserved and z can be 

fitted assuming a constant ctrm and varying So as a function of redshift i.e. 

So= 7(1 + ^)^. 

iii. ) Conversely, the observations may also be fitted assuming a constant So and 

varying aRM as a function of redshift, i.e. aRM — 7(1 + z)^- 

I use the Numerical Recipes subroutine AMOEBA (Press et al., 1997) to estimate 
the best fits for the three assumptions presented above. This is done by calculating 
the least squares difference between the predicted o"obscrvcd and the true (Jobserved 
for each source. The best-fit is then given by the values of So and aRM which 
gives the smallest least squares value*. These results are presented in Table 5.2. 
Unfortunately, no best-fit solutions can be given for either Power-law model when 
both free paiaiiictcis arc ccnistant with redshift. This is simply due to the fact that 

■tLeast squares value = (model - observed)^ 
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the Gamma function in equation 5.14 is very sensitive to the exact choice of So and 
and no reasonable fit can be achieved. 

constant auM, constant So 



Gaussian 


CRM 


= 89.1 


So 


= 15.9 


Quadratic 


<^RM 


= 85.2 


So 


= 3.5x10^ 




constant 


CRM, varying So 






Gaussian 




= 523.4 


So 


= 0.03(l+z)5-5 


Quadratic 


C^RM 


= 173.3 


So 


= 1.20(l+z)3-2 


Power-law with m 


= 2 aRM 


= 529.6 


So 


= 0.01(l+z)6i 


Power-law with m 


= 4 aRM 


= 504.6 


So 


= 0.03(l+zf-2 




constant 


So, varying ctrm 






Gaussian 


CRM 


= 6.1(l+z)4-^ 


So 


= 5.5 


Quadratic 


CRM 


= 4.9(l+z)4-6 


So 


= 2.0 


Power-law with m 


= 2 aRM 


= 6.9(l+z)4-8 


So 


= 2.4 


Power-law with m 


= 4 aRM 


= 5.4(l+z)4-4 


So 


= 12.0 



Table 5.2: Parameters of best fit. So in units of kpc and aRM in units of rad m 
All results are in the sources' frame of reference. 



5.3.1 Constant So, constant aRM 

Figure 5.1 demonstrates that the models with constant So and constant aRM do not 
adequately describe the data. Both models, despite having shghtly different aRM 

and drastically different So, have similar solutions. The only discernible difference 
occurs at very low rcdshift where the Quadratic model quickly falls to zero. It is 
obvious then that there must be some evolution with redshift in either So or aRM 
to account for the observed correlation between (Jobserved and redshift. 



5.3.2 Varying So and grm 

Figure 5.2(a) and 5.2(b) present the best fit solutions to the Gaussian auto-correlation 
model and the Quadratic model when So is varied keeping urm constant and also 
when aRM is varied keeping So constant. Figures 5.2(a) and 5.2(b) use the RM 
dispersion of the lobe and that averaged over the entire source, respectively when 
calculating the best-fit solutions. It is obvious that the solutions, although using 
subtly different data, produce almost identical fits. For this reason I use only the 
average source values in all further fits. 

The best-fit solutions to the Power-law model with exponents m = 2 and m = 4 
are shown in Figure 5.3. Gomparing all 8 solutions presented in Figures 5.2(b) and 
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Figure 5.1: Best fit solution of the Faraday screen keeping So and a^M fixed for the 
Quadratic and Gaussian models. The cXobserved values for all lobes are plotted as 'x'. 
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Figure 5.2: Best fit solutions of the Faraday screen for the Quadratic and Gaussian 
auto-correlation models. Line styles indicate whether So or ct^a/ are being varied 
as a function of (1+z). (a - left) Solution uses all lobe (Tobserved values in the and all 
lobe (Jobserved values are plotted, (b - right) Solution uses only the average (Jobserved 
values for each source and the average source cXobserved values are plotted. 



5.3 demonstrates that by using the (Xobserved from observations it is impossible to 
determine whether or not it is So varying with redshift or o^u varying with redshift 
that is causing the observed trend of (Tobserved with redshift. All 8 solutions are 
statistically identical in terms of their least squares fit value. The only noticeable 
difference occurs at very low redshifts when So varies as a function of redshift and 
grm is kept constant. Each of the different models shows a slight upturn as z — >-0. 
This is simply due to the fact that as z — >0, the beam size approaches zero and 
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hence a — ^ oo and so o"observed — ^rm- This is similar to the case when a source is 
completely resolved. 
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Figure 5.3: Best fit solutions of the Tribble models for the Faraday screen for a 
Power-law model with exponents m = 2 and m = 4 plotted against the average 
^observed values for each source. 



5.3.3 Depolarisation 

Although there is no significant difference between the solutions with varying So as 
a function of redshift compared to varying ajiM, I folded the results presented in 
Table 5.2 into the corresponding expressions for polarisation given in equations 5.7 
to 5.13. This allows the fractional polarisation that would be observed at 4.8 GHz 
and also at 1.4 GHz to be calculated. These values can then be used to calculate 
the depolarisation parameter DMf;| (see section 3.3). 



5.3.3.1 Constant Sq, constant ajiM 



Figure 5.4 shows that the calculated Quadratic depolarisation keeping So and aRM 
constant with respect to redshift, is always 1. It is obvious that this model fails to 
describe the data. When the Gaussian auto-correlation model is used to calculate the 
depolarisation, it predicts more depolarisation at lower redshifts, underestimating 
the higher redshift depolarisation measurements. Both models fail to describe the 
data confirming my earlier finding that there must be some evolution in either So or 

CRM- 
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z 

Figure 5.4: Model prediction for the depolarisation with constant So and ct/jm ob- 
tained from the best-fit solutions using the Quadratic and Gaussian auto-correlation 
models. The data points represent the average observed source depolarisation. 



5.3.3.2 Varying So and a^M 




Figure 5.5: Model predictions for depolarisation using all models with a constant 
aRM and a varying So with redshift, obtained from the best-fit solutions in Table 
5.2. Data points are as in Figure 5.4. 



All solutions presented in Figures 5.2 and 5.3 are equally significant independent 
of whether or not So or a^M is varied. The solutions are also independent of which 
model is used to calculate the observed RM variations. However, there is consider- 
able difference between the calculated depolarisations as shown in Figures 5.5 and 
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5.6. Figure 5.5 uses only the solutions which have a constant ctrm and varying Sq. 
In all cases the depolarisation is widely over-estimated. In the most extreme cases 
predicting depolarisations in excess of 30. Such high depolarisations have never been 
observed, even at high resolution (e.g Best et al., 1997; Carilli & Barthel, 1996; Best 
et al., 1998; Dennett-Thorpe et al., 1997). In the case of the Quadratic model the 
solution even predicts a strong anti-correlation of depolarisation with redshift which 
has never been observed. In striking contrast Figure 5.6, which shows the results 




I I I I I I 

0.5 



Figure 5.6: As Figure 5.5 but with a constant So and varying ctrm with redshift. 

for constant So and varying cr/jjv/? is a much better fit to the observations. At lower 
redshifts {z < 0.5) all models are equally likely and the scatter in the depolarisation 
encompasses all the solutions. It is only when the higher redshift [z >0.5) sources 
are considered that I am able to determine a significant difference between the mod- 
els. Both the Power-law model (with an exponent m = 2) and the Quadratic model 
overestimate the observed depolarisation. However, the m = 4 Power-law solution 
and the Gaussian auto-correlation solution give results that are consistent with the 
observations. The fact that there is some degree of scatter at all redshifts suggests 
that a fixed So is a simplification and So may vary slightly from source to source. 
However, the dominant effect is the variation of ctrm- 

In the next section I consider the effects of varying So slightly using the Gaussian 
auto-correlation model and the Power-law model with m = 4. 



5.3.4 Scatter in So 

Figure 5.7(a) shows the Gaussian auto-correlation model with a^M = 6.1(1 + z)'^-'^ 
rad m~^ with So varying around the best-fit value, So=5.5 kpc. By allowing So to 
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vary by a maximum of ±90% from the best fit value I can determine the effects 
of changes in So on the data. Figure 5.7(b) is the m = 4 Power-law model with 
c^RM = 5.4(1 + zY'^ rad m~^, again So is allowed to vary around the corresponding 
best-fit value of So=12.0 kpc. In both cases the low redshift results {z <0.5) are 



^o^^' / ^o^ 

s =0.5s ; ^' s =0.5s 

s =0.9s / ^ _ _ _ s„ = 0,9s. 




z z 



Figure 5.7: (a - left) Varying the value of So using the Gaussian auto-correlation 
model and (b - right) the Power-law model with m = 4. In both plots the bold line 
indicates the initial So taken from Table 5.2. 

largely insensitive to changes in Sq. This is not the case at the higher redshifts 
{z > 0.5). As So^ ±90%So the depolarisation becomes increasingly unrealistic. 
In both models using So=0.1so produces large over-estimates in the high redshift 
depolarisations. This is also true when So=0.5so using the Gaussian auto-correlation 
model. This suggests that although So doesn't vary with redshift, a small amount 
of variation can describe the scatter in the observed data seen in the DM-z plots of 
the previous chapters. Figure 5.7(b) is the best fit to the data suggesting that the 
Power-law model with m = 4 is perhaps the most accurate and the small variations 
in So do aptly describe the data. This model predicts p ~ A^^ compared to the 
Burn-law p ~ A~^. 

It must be noted that the fitting procedure only uses 23 out of a possible 26 
sources^. To gain a better estimate of which model is the most accurate a much 
larger sample, preferably with higher resolution, would have to be used. By using 
a sample of radio galaxies at high resolution the variations in So would be less 
significant and the sources would be well resolved and hence cr/?A/~ o"observed- 



5.3.5 Using So=10 kpc 

As I have already noted my RM maps are not of sufficient quality to determine how 
the RM structure changes over each lobe of a source. In several sources (e.g. 3C 299 

§As noted in the previous chapters 3 sources have no RM information in at least one lobe. 
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and 6C 0943+39) I have only a few pixels of RM data in each lobe. In section 5.3.3 
I demonstrated that keeping So constant and varying cfrm as a function of redshift 
is the most realistic of all the solutions presented. In all cases So was found to be 
less than 12 kpc. 

The method described in the previous sections does not include any information 
on how the depolarisation differs with changes in radio-luminosity. Figures 5.2 to 
5.3 show that the two high redshift samples ('x' = low radio-power sample A and 
'o' = high radio-power sample B) are equally distributed around the solutions. This 
method does not enable any determination of how depolarisation is affected by 
changes in P151 and I^iobe- To overcome this problem, I set So=10 kpc and calculate 
o'rm using the observed rms variation in the rotation measure, (Tobservcd; and the 
known beam size. This calculated grm is then used to work out the depolarisation 
for the source if it was located at 2; = 1. Only by doing this can I determine if any 
depolarisation trends are really due to redshift effects and not simply due to the 
fact that a third of my sample is located at z — 0.4. This is the same method as 
used in section 3.4.2 to shift the depolarisation using Burn's law. 

Quadratic model results 

rz,DM = 0.67857 t = 6.13 
rp^uM = 0.49405 t = 3.77 
^D,DM = -0.21899 t = -1.49 

Gaussian model results 
r^.DM = 0.67894 t = 6.13 
rp,DM = 0.49541 t = 3.78 
rD,DM = -0.22294 t = -1.52 

Power-law model results 
rz,DM = 0.67771 t = 6.1 
rp,DM = 0.49171 t = 3.75 
^D,DM = -0.23059 t = -1.57 

Table 5.3: Spearman rank results for So=10 kpc 

As Table 5.3 demonstrates, I always find a strong DM^^i-z correlation with a 
weaker correlation with radio-luminosity, independent of the model used. There is 
never a significant anti-correlation with size. These results are consistent with the 
results found in the previous chapters, obtained from Burn's law. 



5.4 Conclusions 

I have used the T91 models to determine how the size of the cells in an external 
Faraday screen, So, varies with redshift and also to determine if there is any evolution 
in the intrinsic RM variations, (Trm, with redshift. 
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I have compared the observed RM variations with the predicted RM variations 
using a constant Sq and cfrm with redshift and I find that this model is inconsistent 
with the data. There is no statistical difference between the models when So is varied 
compared to when (Trm is varied. Only by comparing the predicted depolarisation 
with the observed depolarisation is there any distinction between the models. In 
all cases when So is varied, keeping grm constant, the predicted depolarisation is 
over-estimated. Only by keeping So constant and varying Grm is the predicted 
depolarisation similar to the observed depolarisation. 

The Gaussian auto-correlation model and the Power-law model with exponent 

m = 4 are good fits to the data. When So is allowed to vary slightly in these models 
I find that it has little effect on the low redshift sources, but as redshift increases the 
size of So becomes important. A small amount of variation (typically < 50%) can 
be introduced to the models to allow for the obvious scatter in the depolarisation- 
redshift plots. In all cases it is the higher redshift sources that constrain the models. 

A fixed cell size and varying ctrm with redshift gives further evidence that the 
observed depolarisation is due to the increasing disorder with redshift and not due 
changes in the density of the environment with redshift (see previous chapter) . 

By setting So=10 kpc and calculating (Jrm using the observed (^observed and the 
known beam size 1 am able to determine how the depolarisation is affected by 
changes in radio-luminosity and radio-size. Irrespective of the model I use 1 always 
find that depolarisation correlates strongly with redshift and with radio-luminosity. 
In no cases do I find an anti-correlation with size. These findings agree with the 
results presented in sections 3.4.2 and 4.3.2. Although the Burn's law over-estimates 
the degree of depolarisation, 1 find similar results using the T91 models, which are 
more realistic than the Burn model. This proves that the simplistic Burn model 
can still be used to find general trends with depolarisation in a sample, irrespective 
of the fact that it may overestimate the depolarisation. Interestingly I find that 
p ~ best describes the data, compared to the widely used Burn p ~ result. 



Chapter 6 



Modelling FRII sources 



In the previous chapters I have used depolarisation and the variations in the rotation 
measure to sample the environment around powerful FRII radio galaxies. However, 
In chapter 4 I found evidence that these quantities do not sample the density of the 
source environment (e.g. Gregory & Condon, 1991), but are indicators of the degree 
of disorder of the Faraday screen. 

In the following I use the model for the evolution of the radio luminosity of FRII- 
type sources as a function of source age by Kaiser et al. (1997, hereafter KDA) to 
model the source spectra. This allows a more theoretical approach to determine 
source properties such as the age of the jets and their energy transport rate or jet 
power. The model also provides an estimate for the density of the gaseous source 
environment. By comparing observations of the source depolarisation and variations 
in the RM with model estimates of the density, jet-power and lobe pressure allows 
an independent check on the trends of the observational properties of the source 
environment derived in the previous chapters. 

I have obtained depolarisation and rotation measure observations for only a small 
sample of 26 sources. Because the KDA model can be used to model the radio 
spectra of a source I can apply it to the complete 3CRR, 6CE and 7C III samples. 
In practice I use flux measurements at 3 different observing frequencies, for each 
source, to constrain the KDA model. I then analyse how the model parameters (jet 
power, density of the environment, source age and lobe pressure) evolve as a function 
of radio-luminosity, size and cosmological epoch. This gives a more statistically 
significant result on the overall evolution of the environments of powerful FRII 
sources. 
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6.1 The model for FRII sources 

The large scale structure of FRII sources is formed from twin jets emerging from 
a central AGN buried inside the nucleus of the host galaxy. Each jet propagates 
outward from the core of the source and terminates in a strong shock giving rise 
to the so-called radio hotspots. The jet material inflates a lobe surrounding the 
jet, which drives a bow shock into the surrounding medium. This model was first 
proposed by Scheuer (1974). Kaiser & Alexander (1997), hereafter KA, showed that 
in a purely dynamical model of FRII evolution, the bow shock and lobe grow self- 
similarly, which has also been inferred from observations (e.g Leahy & Williams, 
1984; Leahy et al., 1989). 

The dynamical model of KA assumes a simple power-law for the external density 
distribution around the radio source, p^,. 



where d is the distance from the central AGN, Oo is a scale height, po is the density 
at d = Oo and < /? < 2. 

The KA model also assumes that the rate of injection of energy into the lobe, 
defined as the jet power, Qo, is constant over the lifetime of the source. Due to the 
high sound speed in the lobes, the lobe pressure is uniform throughout except for 
the regions close to the hotspots. It is the lobe pressure, piobe , that confines the jets. 

The KA model predicts that the length of a single jet, Diobe, grows with time as 



where Ci is a constant which depends on the geometry of the lobe, the thermody- 
namic properties of the jet material and the gas in the source environment (see KA 
for details), t is the age of the source. 

KDA added a prescription for calculating the synchrotron radio emission of the 
lobes to the dynamical model of KA. For this the lobe is split into small volume 
elements, SV. These are injected into the lobe at time ti containing a population of 
relativistic electrons with a power-law energy distribution with exponent p ranging 
from 7 = 1 to some 'jmax- Here, 7 is the Lorentz factor of the relativistic electrons. 
This electron distribution is then subject to energy losses due to adiabatic expansion, 
synchrotron radiation and inverse Compton scattering of the CMB. Since ti will 
be different for different volume elements, the resulting energy distribution of the 
electrons in each SY at time t is a function of the injection time t^. 




(6.1) 




(6.2) 
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The total radio emission of the lobe at a given frequency, i^, is the sum of all the 
contributions from the volume elements, 

P.= I ^(JTCUB—n{-i)5V{t) (6.3) 
J on V 

where ut is the Thomson cross section, is the energy density of the magnetic 
field and ^(7) is the density of relativistic electrons with Lorentz factor 7. 



In the model of KDA equation 6.3 becomes. 



67rz/(r + l) (i-«i/3-a2(t,ti)7)'"'' \^ 



dti, (6.4) 



where Ho is the normalisation of the electron energy distribution, r is the ratio 
of the energy density of the magnetic field and that of all particles at the time of 
injection, ti. R is the ratio of the length of the lobe and its width at half the lobe 
length. The various F are the adiabatic indices of the material in the lobe, Fc, and 
the magnetic field, which is assumed to be tangled on small scales, F^. a2{t,ti) — 
(t-«i/3)/^_ (tpi/-^)/^,. where ai is the exponent of the power-law time dependence 
of the volume elements, i.e. dV oc t"^, which KDA find to be Oi = (4-|-/3)/[Fc(5 — /?)]. 
Finally, 7j is the Lorentz factor at ti of electrons which have a Lorentz factor of 7 
at the current time, t. 

There is no analytical solution for equation 6.4 for general p. However, for p — 2 
the integral simplifies considerably to give 



7/4 

p ^ Plobe 



1 - 



mm 



(6.5) 



where s = ai (i + ^) + (|^), and 

<7T r^l' (r. + 1) (2i?) Vrc Bl^^ [/n(7,™a.) + 7^. - l]"^ / 27r^ \ 



G — 

1087rz/V2 c^^-P) cmT (F, - 1)^/^ {k' + 1) (r + 1)V^ (5 - /3)-2 

(6.6) 

where p;o6e is the pressure in the lobe at the current time t. a\ determines 03 = 
1 — ai(F5 + 1/3) and 04 = 1 — ai/3. F^ is the adiabatic index of the gas in the 
source environment. The lobe may contain non- relativistic particles, the ratio of 
their energy density and that of the relativistic electrons is given by k' . Finally, 
is the electron mass and Ho is the magnetic permeability of the vacuum.* 



Older parts of the lobe will no longer radiate at the observing frequency v due to 
severe radiative losses, so it is not realistic to integrate equation 6.4 over the entire 
age of the smucc. t„i,„ is the minimum injection time for volume elements to still be 

*In the derivation of equations 6.4 to 6.6 I have used the results of KA and KDA. 
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emitting radiation at the observing frequency v. I assume that the electrons only 
emit at their critical frequency 'y'^{t,ti)i'L(t,ti) where is the Lamour frequency, 
and use the condition that u = '^'^{t,tmin)^L{t-,imin) to give an implicit equation for 



V — 



g^-(2ai/3) 



(6.7) 



where X = (r piobet^''"'^)/[{i^ + 1) i^c — l)]; Uc is the energy density of the cosmic 
microwave background, Uc oc (1+2;)'^, (see also KDA). There is no analytical solution 
for train, but it is possible to numerically find a tmin, for any given source age and 
frequency, that solves equation 6.7 to a high degree of accuracy (i.e. within 1 x 10~^ 
or less.) 



The model calculates all the properties in the sources' frame of reference. This 
means that the observing frequency and Pjy have to be transformed using the chosen 
cosmology, see appendix D for details. 



6.1.1 Application of the model 



The model for an individual radio lobe depends on a number of parameters, some 
of which are determined from the radio observations. These are the observing fre- 
quency, u, the length of the lobe, Diobe, the aspect ratio of the lobe, R, and the 
monochromatic radio luminosity, Pj,. The source redshift is taken from archival ob- 
servations, see section 2.1 and appendices E to G. Another set of parameters are 
not directly accessible by observations, but their values are either well-constrained 
in general or are simply set to reasonable values. The adiabatic indices are set to 
Ta; = Tc = 5/3, assuming that the material in the lobe is mainly non-relativistic. 
The magnetic field in the lobe is assumed to be tangled on scales much smaller 
than Diohe: therefore = 4/3. Any particles in the lobe other than the electrons 
responsible for the emission of synchrotron radiation are neglected, so k' = 0. I 
further assume that the relativistic electrons and the magnetic field are initially, 
at time ti, given by their values appropriate for minimum energy conditions (e.g. 
Longair, 1994). Thus r = 3/4 for the chosen slope of the electron energy distribu- 
tion with p — 2. I set the high-energy cut-off of this distribution to ^i^max — 10^^ 
and '-fi^min = 1- It is clear from equation 6.6 that the exact value of ^i^max does not 
significantly infiuence the results. Finally, the gaseous atmospheres around powerful 
radio galaxies imply 1 < f3 < 2 and so 1 set /?= 1.5. This leaves three free model 
parameters: The current pressure in the lobe, piobe, the age of the source, t, and the 
minimum injection time, tmin- 
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Diobe, the lobe ratio and the rest frame frequency of the observations are used 
to determine the constant G (equation 6.6). Although piobe and t are the principal 
free parameters of the model, in the fitting process afpo is used instead of the lobe 
pressure. This is due to the fact that piobe and ty (the age of a source in years) are 
highly anti-correlated, see Table 6.4, and applying the model using these parameters 
could create a bias in the data. As Table 6.2 shows there is no significant a^po-ty 
correlation, making these two parameters a good choice for the input parameters 
for the fitting process. Initial 'guesses' for a^po and ty are used to calculate the lobe 
pressure, which depends on the jet length and the lobe ratio. 

^'"'^ 4i?2(r. + 1)(5 - (5Y \Qf-^y^ J ^^-^^ 



2(re+l)(5-/3)2i?2ty2 

The age of a source and its density, afpo, are used to calculate the minimum 
injection time, tmin, (equation 6.7) and the monochromatic luminosity, P^, using 
equation (6.5). This monochromatic luminosity is converted to a specific flux, 8,^ 
using equation (D.4). Comparing S^^ with the observed flux it is possible to obtain 
a xl deviation of the model at frequency u. This gives three measurements of xt at 
the three different frequencies. The Numerical Recipes subroutine AMOEBA (Press 
et al., 1997), is used to find the best estimate for ty and afpo for a lobe by minimising 

X2 _ ^2 _j_ ^2 _j_ ^2 
A.i'i I Az/'2 ' At/3 ' 

The jet power, Qo, is then calculated from the final values of ty, a^po and hence 

Plobe, 



Qo = "^'""^^^^^^^.^"^^ (6.10) 



For full details of the above model see KA and KDA. 



6.2 Modelling the small sample of 26 sources 



The model uses the three flux measurements around 1.4 GHz, 1.6 GHz and 4.8 GHz 
and their corresponding errors to constrain the free parameters of the model. Diobe, 
is calculated from the observed angle of the source on the sky, and the measured 
redshift using the chosen cosmology. The lobe ratio is taken from the total intensity 
maps. 
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6.2.1 Angle to the line of sight 

It is not possible to determine from the radio data whether or not the jets in a 
given source lie in the plane of the sky or whether the measured length of the lobe 
is projected. The length of the jet and the lobe ratio are both affected by this angle 
dependence, however it is the jet-length that is most affected. 

To test how sensitive the model is to a change of the viewing angle, each source 
is modelled with a number of viewing angles ranging from 10° to 90° to the line-of- 
sight. I found that in most cases the source is most probably angled around 57.3° 
to the line-of-sight, which is the average angle expected for a sample of randomly 
oriented sources. In these cases there is no significant difference in the results for 
any orientation unless the source is angled to less than 23° to the line-of-sight. 
These very small angles produced highly unlikely results, i.e. extremely large sizes, 
improbable jet-powers and/or extremely high densities. Sources which are already 
large, e.g. 3C 46, are unlikely to be orientated at angles much less than 57° as this 
makes the source larger than a few Mpc across. In several cases when the sources 
became large (small orientation angle) the source age was not consistent with the 
source size, i.e. less than 10^ years for sources of order of a Mpc across. 

There are a few sources that all angles seemed equally probable (e.g. 6C 0943+39 
and 3C 280), but these are in the minority. To assume a different angle for each 
source is impractical, especially in section 6.4 where 211 sources taken from the 
combined 3CRR, 6CE and 7C 111 samples are modelled. The model does not directly 
depend on the line-of-sight so it is not possible to determine an orientation angle 
from the fitting procedure. 

For a sample of randomly orientated sources the most likely orientation is 57.3° 
to the line-of-sight and only a few of the sources will be orientated at angles smaller 
than 23°. Thus these sources will not have a large effect on the statistics. In all 
the following sections the results have been calculated using an orientation angle of 
57.3° for all sources. 



6.3 Results 

The KDA model provides an estimate for the jet-power, density, lobe pressure and 
source age for each lobe in the small sample. However, in the next section I only 
consider how the density parameter, afpo j relates to changes in depolarisation and 
variations in the rotation measure. How the model parameters are affected by 
changes in the 'fundamental' parameters (i.e. redshift, radio-luminosity and source 
size) is presented in detail in section 6.4 using the combined 3CRR, 6CE and 7C HI 
samples. By using a large sample any results found will be more significant than 
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with only 26 sources. However, the large samples do not have depolarisation and 
rotation measure information for each source. Therefore I use my smaller sample to 
determine if depolarisation and density are related. 



6.3.1 Comparing a^po with DM^ and cthm^ 

In section 4.3.2 I found that the Burn corrected depolarisation measure, DM^, cor- 
related strongly with redshift. By analysing the minimum energy magnetic field 
asymmetries in section 4.4.2.3 I found evidence that the observed depolarisation of 
a source was unaffected by changes in the density of the environment. The KDA 
model gives an estimate of the gas density around each lobe which can be compared 
with the source depolarisation. 



C) o o 




Figure 6.1: (a - left) Unshifted depolarisation against a^po and (b - right) depolar- 
isation shifted using Burn's law. Symbols are as in Figure 4.1. 



Figure 6.1 shows that there is no relationship between depolarisation and a^po, 
independent of whether I use the observed depolarisation (Figure 6.1(a)) or I use 
the shifted depolarisation from section 3.3 (Figure 6.1(b)). Figure 6.2 shows that 
there is also no trend between the variations in the RM in the sources' frame of 
reference, (Jrm^, and a'^Po- 

Figure 6.3 shows a weak correlation between radio luminosity and afpo , but no 
corresponding correlation with redshift. The relationships between, z, P151 and 
a^Po are discussed in detail in section 6.4.6 using the large sample. However, it 
is worth noting that in the small sample there is no indication of any evolution of 
the environment density with redshift. Thus the evolution of depolarisation with 
redshift cannot be attributed to changes in density. In section 4.3.2 I also found 
a strong P151-DM2 correlation that was comparable to the z-DM^ correlation. As 
the Partial Spearman results show, at any given radio-luminosity there is still no 
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Figure 6.2: (T/jm^ against a^po- Symbols are as in Figure 4.1. 




Figure 6.3: (a - left) Radio-luminosity against afpo and (b - right) redshift against 
afpo- Symbols are as in Figure 4.1. 

significant DM^-afpo correlation. 



r 13 n,, D = —0.119 

agpo DNh, Pi5i 

r fj p r,., = 0.526 

doPo Pl51, DM^ 

Tp /5 = 0.504 

Pl5l DMz, aZpo 



D = 0.83 
D = 4.05 
D = 3.84 



At any given density the Pisi-DM^ correlation still exists, thus any changes in 
DM2 with either z or P151 are independent of changes in the density of the source 
environment. The lack of correlation between afpo and the observed polarisation 
parameters gives further evidence that depolarisation and ctrm^ should not be used 
as tracers of the density of the environment (e.g Burn, 1966; Garrington & Conway, 
1991; Blundell et al., 1999) and should only be used as a measure for the disorder 
in the environment. 
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6.3.2 Minimum energy, a"po and p;o6e 



Figure 6.4(a) shows that there is a tight correlation between the ratio of pressures in 
the lobes of a given source and the minimum energy magnetic field ratio. The only 
notable exception being 3C 16, which has a high minimum energy magnetic field 
ratio, but a relatively low ratio of lobe pressures. This agrees with the suggestion 
that 3C 16 is in the process of switching back on (Harvanck & Hardcastle, 1998, also 
see section 4.4.2.2). Figure 6.4(a) demonstrates that the simple minimum energy 
magnetic field model can be used to give an estimate of the lobe pressure. For the 
sources in the large sample, discussed in section 6.4, there is no information on the 
individual lobes of each source so it is impossible to estimate the minimum energy 
magnetic field for individual lobes. However, as the minimum energy magnetic 
field and Y>iobe are tightly correlated, I can be confident that the results using the 
KDA model give a good estimation of the magnetic field in a source. 
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Figure 6.4: (a - left) Ratio of minimum energy magnetic field against the ratio of 
lobe pressures and (b - right) against the ratio in a^po- Symbols are as in Figure 
4.1. 



Figure 6.4(b) shows that there is a weaker correlation between the density ratio 
and the minimum energy magnetic field ratio which corresponds to a significance of 
~99.5%. By removing 3C 351 and 3C 16 the correlation becomes >99.9% significant. 
In section 4.4.1 I noted that 3C 351 has a very distorted northern lobe, significantly 
extending to the west of the double hotspots. This could be attributed to a very 
low density on the western side of the north lobe, allowing the source to propagate 
westward. The southern lobe is much more constrained, implying a larger density, 
which could be the cause of the large ratio of densities between the two lobes. The 
fact that the minimum energy magnetic field ratio and a^po ratio are correlated 
suggests that it is the difference in the environment of the two lobes of a source that 
is causing the difference in the minimum energy magnetic field and not beaming as 
suggested in section 4.4.2.3. 
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6.4 Modelling the complete 3CRR, 6CE and 
7C III FRII samples 

6.4.1 Observations 

The large sample contains a sample of 211 FRII sources taken from the complete 
3CRR, 6CE and 7C III samples, see Figure 6.5. Only sources larger than 10 kpc 
are included in the analysis as below this limit the model may not be applicable 
(Alexander, 2000). For consistency all sources smaller than 10 kpc are excluded 
using the stated cosmology. 

Appendices E to G give the archival redshift, size, flux measurements, lobe ratio 
and their associated references for each source used in the analysis. In section 6.2 I 
stated that the model uses three flux measurements and their corresponding errors 
to constrain the free parameters. I apply the model to the flux measurements at 
178 MHzt, 365 MHz and 1400 MHz for the 3CRR sample; 151 MHz, 365 MHz and 
1400 MHz for the 6CE sample and 151 MHz, 327 MHz and 1490 MHz for the 7C HI 
sample, aiow is defined as the spectral index between the lower frequency and 365 
MHz and ahigh is defined as the spectral index between 365 MHz and 1400 MHz*. 

In section 6.1 I use p = 2 in the determination of the monochromatic radio 
luminosity, P^. This means that only sources with a spectral index, aboundary < —0.5 
(including errors) can be accurately modelled where aboundary is determined by the 
choice of p, i.e. p = 2q;+ 1 (Longair, 1994). For this reason 3C 382 and 6C 0922+36 
are excluded, since they have ahigh = —0.22 and ahigh — —0.30 respectively, which 
even when the errors are accounted for have ahigh > Oiboundary 



6.4.2 Application of the model to the combined sample 

The source size, G source, is calculated from the observed angle of the source on the 
sky, and the source redshift. The lobe ratio for each lobe (i?iobci and i?iobe2) is taken 
from high resolution maps, where they exist, in the literature (see appendices E to 
G). In the few cases where there are no high resolution maps, particularly the faint 
7C III sources, the lobe ratio is set to be the average sample lobe ratio, using the 
sample that the source is taken from. 

'^The different lower frequency of the 3CRR sample is simply due to the fact that it is defined 
at 178 MHz instead of 151 MHz. 

tPor the 7C III sample I use 327 MHz and 1490 MHz instead of 365 MHz and 1490 MHz. This 
is due to the fact that is only at 327 MHz and 1490 MHz that there is flux information for each 
source in the sample. 
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Figure 6.5: Redshift against the radio luminosity at 151 MHz. The 3CRR sample is 
represented by 'x', the 6CE sample by 'o' and the 7C III sample by '*'. A spectral 
index of -0.75 is used to shift the 3CRR luminosity data at 178 MHz to 151 MHz. 



For every source the flux of the lobe, at each frequency, is taken to be half the total 
flux. This a simplification, but there are no flux measurements for individual lobes, 
at all 3 frequencies, for each of the 211 sources in the literature. For this reason the 
lobe length is also taken to be Diobe= DsoMrce/2 and the lobe ratio is averaged over 
the source i.e. -Raverage = (-Riobei + -Riobe2)/2. This cusured that the results are not 
biased towards sources that are well observed and have high resolution maps (e.g. 
the majority of the 3CRR sample). This is particularly important when dealing with 
the sources of the 7C HI sample for which very little archival information exists. 



6.4.3 Fundamental Parameters 



In the following sections I investigate the relationship of the model properties of 
my sources with the 'fundamental' properties redshift, z, radio luminosity at 151 
MHz, Pi5i, and physical size, Diobe- Before I can understand how the model param- 
eters evolve with the fundamental properties I need to understand the relationship 
between the fundamental properties. This is similar to the analysis in section 4.2. 
Table 6.1 contains the associated Spearman Rank results of the fundamental param- 
eters. As expected, I find redshift and radio luminosity to be highly correlated (see 
Figure 6.5), to a significance of > 99.99%, but this is simply due to the flux limits 
of the samples and cannot be avoided. Using the statistical techniques described in 
chapter 4.1 this correlation can be isolated from any other correlations. 

By using the Partial Spearman Rank test I find that even with 3 large complete 
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Parameters 


rs value 


t value 


Z Pl51 


0.69302 


13.93 


Z Djobe 


-0.40682 


6.45 


Pl51 D«obe 


-0.29855 


4.53 



Table 6.1: Spearman Rank values for z, P151 and D/obe in the combined samples. 



samples there is no significant Pi5i-D;of,e correlation and that the size-redshift anti- 
correlation is weak when it is compared to the strength of the redshift - radio 
luminosity correlation. 



r.A..e,Pi5i= -0.288 D = 4.27 

r.P,,„A„,.= 0.655 D = 11.27 

^Pi5iA..e,.= -0.019 1^ = 0.27 



This is also evident in Figure 6.6(a), the weak Pisi-D^ofte ant i- correlation is almost in- 
distinguishable from the scatter. This confirms that even with such a large sample of 
sources there is only a weak relationship between the low frequency radio-luminosity 
and lobe size and agrees with the findings of Barthel & Miley (1988).^ 
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Figure 6.6: (a - left) Source size against radio-luminosity and (b - right) against 
redshift. Symbols are as in Figure 6.5. 



As already noted, redshift and lobe size are only weakly related, see Figure 6.6(b). 
The weak correlation can be see in Table 6.1 and also in the Partial Spearman 
Rank results above. However the z-Diote correlation is much weaker than the Z-P151 
correlation and will have little effect on the model results. 



The Pi5i-Djobg anti-correlation is discussed in detail in section 4.2. 



6. Modelling FR II sources 



-106- 



6.4.4 Source age, ty 



An important test of the model is to determine if the age of a lobe, ty, correlates 
with its size, Diobe- Clearly larger sources should be older. As Table 6.2 and Figure 
6.7 demonstrate, Diobe and ty correlate to a high degree, > 99.99% significant. The 
PCA results shown in Table 6.3 are dominated by this strong correlation and do not 
show any other significant trends. A strong D;o;,e-ty correlation is expected from the 
model, but it is the scatter caused by redshift and radio-luminosity that is important 
in determining how the source environment affects its age. 
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Figure 6.7: Source age against source size. Symbols are as in Figure 6.5. 



Parameters rs value t value 



ty Biobe 


0.80452 


19.63 


ty z 


-0.40583 


6.43 


ty Pi5i 


-0.29238 


4.43 


ty afpo 


0.01446 


0.21 


ty aiow 


-0.24527 


3.67 



Table 6.2: Spearman Rank values for the age of a lobe 



Table 6.2 shows an anti-correlation of the source age with redshift and also a 
slightly weaker anti-correlation with radio-luminosity. Figure 6.8(a) shows that the 
Pi5i-ty correlation is weak. The trend with redshift could be explained by the 
inclusion of low-redshift giants. In section 6.4.3 I have shown that there is a weak 
z-Diobe anti-correlation which combined with the strong evolution of source size with 
age could cause the apparent z-ty anti-correlation seen in Figure 6.8(b). 

The Partial Spearman Rank analysis shown below demonstrates, that for any 
given radio- size the z-ty correlation is barely significant. Interestingly at any given 
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Parameter 


1 


2 


3 


4 


z 


0.378 


0.575 


0.724 


0.054 


Pl51 


0.312 


0.658 


-0.685 


0.001 




-0.625 


0.318 


0.021 


0.713 


ty 


-0.608 


0.368 


0.078 


-0.699 


Eigenvalue 


53.09% 


30.96 % 


13.70% 


2.25 % 



Table 6.3: Eigenvectors and Eigenvalues for the age of the lobe 
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Figure 6.8: (a - left) Source age against radio-luminosity and (b - right) against 
redshift. Symbols are as in Figure 6.5. 



source age, the z-D/obe also becomes barely significant. This suggests that the ob- 
served z-Diobe correlation (and hence also the z-ty correlation) is simply a by product 
of the Malmquist bias arising from the flux-limited nature of the samples. Only by 
observing at fainter flux-limits, e.g TOOT (Hill & Rawlings, 2003), can this obser- 
vational bias be removed. 

r.A„,„t,= -0.148 D = 2.15 

r.t,,A„,. = -0.145 D = 2.10 

TD,^,^ ty, . = 0.766 D = 14.57 

It is also worth noting that when I use Partial Spearman rank to determine the 
relationship between redshift, radio-luminosity and age the Pi5i-ty anti-correlation 
vanishes at a fixed source age. 

ty, Pi5i = -0.295 D = 4.38 

rzp,,„ty= 0.657 D = 11.36 

rp,,,ty,z= -0.017 D = 0.2A 

This suggests that both the z-ty and Pi5i-ty anti-correlations are due to the strong 
Dio6e-ty and z-Pisi correlations. Sources at higher redshifts and/or higher radio- 
luminosities are not younger than their lower redshifts/lower radio- luminosity cousins. 
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Table 6.2 also shows that there is a weak anti-correlation of the low-frequency 
spectral index with the age of the source. However, I would expect ahigh to show 
the stronger correlation over aiow due to spectral ageing. Although Table 6.2 shows 
that there is a weak anti-correlation, Figure 6.9(a) indicates that the a/o«)-ty anti- 
correlation is insignificant. Figure 6.9(b) also shows that there is no ahigh-^J anti- 
correlation. Thus, between 151 MHz and 1400 MHz spectral ageing appears to be 
insignificant. 
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Figure 6.9: (a - left) Source age against aiow and (b - right) against ahigh- Symbols 
are as in Figure 6.5. 



6.4.5 Lobe pressure, piobe 



Table 6.4 shows that there is a strong Diobe-Piobe anti-correlation, at a significance 
exceeding 99.99%. This is not surprising as a small source will have a high initial 
lobe pressure, driving the source to expand which in turn will cause the lobe pressure 
to fall. Table 6.4 also shows a somewhat weaker tj-piobe anti-correlation, but this 
is simply due to the strong Diote-ij correlation from section 6.4.4 and the strong 
^lobe-piobe auti-correlatiou already discussed. 



Parameters 


rs value 


t value 


Plobe z 


0.70061 


14.23 


Piobe Pl51 


0.75192 


16.53 


Plobe Diofte 


-0.81074 


-20.07 


Piobe ty 


-0.64674 


-12.29 


Plobe a^Po 


0.67819 


13.37 


Piobe Othigh 


-0.33934 


-5.23 



Table 6.4: Spearman Rank values for the lobe pressure. 



Table 6.4 also reveals a strong Pi5i-p;o6e correlation and a slightly weaker p;o6e-z 
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correlation. Both of these strong correlations can be seen in Figures 6.10(a) and 
6.10(b) respectively. This suggests that the lobe pressure is dependent on both the 
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Figure 6.10: (a - left) Lobe pressure against radio-luminosity and (b - right) against 
redshift. Symbols are as in Figure 6.5. 

redshift and radio-luminosity of a source. By using Partial Spearman Rank on piobe , 
Pi5i and z it is possible to disentangle the relationship with lobe pressure. 
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I find that neither the p/obe-z nor the Pisi-p/obe correlation vanish at any given 
radio-luminosity or redshift, respectively. However, it is worth noting that in both 
the Spearman Rank results (Table 6.4) and the Partial Spearman Rank results the 
Pi5i-p«o6e correlation is always the stronger correlation. 

The Partial Spearman Rank results below show that at any given radio-size, the 
redshift - lobe pressure correlation still exists and is not due to any other independent 
correlation with the fundamental parameters. This suggests that although the P151- 
P/ofee correlation is stronger than the z-piobe correlation, the z-piobe correlation is 
significant and indicates that there is some evolution in the lobe pressure of a source 
with redshift. 



bei Plobe 



Plobe^ Dlobe 



lobe Plobei ^ 



0.386 
0.693 
-0.807 



D = 5.86 
D = 12.28 
D = 16.09 



The result is consistent with Bremer et al. (1992) who find that the pressure around 
^ ~ 1 radio-loud quasars is an order of magnitude higher than the pressures around 
a lower luminosity set of quasars at 2; ~ 0.3. Unfortunately the study by Bremer 
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et al. (1992) cannot determine whether it is the redshift or the radio-luminosity of 
the quasars that is connected to the change in the lobe pressure. 

The PCA results in Table 6.5 show that the first eigenvector holds 48.2% of 
the variation in the data. Suprisingly I find that this eigenvector is dominated by 

the Pi5i-pio6e and z-piobe correlations, with D/o6e-P/o6e the weakest. In the second 
eigenvector, only the Pi5i-p/o6e correlation does not reverse. However, it is the 
Diobe-Piobe Correlation that dominates this eigen-vector. The PCA results indicate 
that the Pi5i-p/o6e correlation is the dominant correlation and that the z-pio^e is not 
highly significant. 



Parameter 


1 


2 


3 


4 


z 


0.518 


-0.182 


0.716 


-0.431 


Pl51 


0.587 


0.322 


0.100 


0.736 


^lobe 


-0.347 


0.845 


0.380 


-0.144 


Plobc 


0.516 


0.386 


-0.577 


-0.501 


Eigenvalue 


48.2% 


22.9% 


18.1% 


10.8 % 



Table 6.5: Eigenvectors and Eigenvalues for the lobe pressure 

The complicated relationship between P151, piobe and Diobe is revealed by using 
the Partial Spearman Rank statistic. I find that the Pi5i-pzobe correlation is almost 
as significant as the T)iobe-Piobe anti-correlation. The fact that both relationships are 
almost equally significant means that PCA is unable to distinguish between the two 
separate relationships. 

m5iA„.e,P.o.e= 0.805 L> = 16.06 

^Pi5i P^o,., A..e = 0.912 D = 22.22 

^Ao.eP.<.6e,Pi5i= -0.931 L> = 24.08 

Using the minimum energy magnetic field argument or the KDA model it is 
possible to show that the density of the lobe environment directly infiuences the 
lobe pressure of a source (e.g Longair, 1994). It is not surprising then that I find 
Plobe and afpo to be correlated, see Table 6.4. The KDA model results demonstrate 
that by deriving the pressure in a lobe, we can obtain an accurate estimate of the 
gas density in the surroundings of a lobe. 

Interestingly Table 6.4 shows a weak anti-correlation between ahigh and piobe , but 
no corresponding relationship with aiow In section 6.3.2 I stated that a stronger 
magnetic field should cause stronger radiative losses in a source and would lead to 
a steeper spectral index. This explains the weak ahigh-Piobe anti-correlation seen in 
Table 6.4. The fact that 1 do not find any relationship of aiow with piobe indicates 
a lower level of radiative losses between 151 MHz and 365 MHz. 
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6.4.6 Density of the source environment, a'^po 



In section 6.3.1 I found that in the small sample of 26 sources there was a weak 
Pisi-^oPo correlation and no significant trend with redshift. Table 6.6 demonstrates 
that with the large sample containing 211 sources I find a much stronger Pisi-afpo 
correlation (see Figure 6.5(a)) and a weak z-afpo correlation (see Figure 6.5(b)). 



Parameters 


rs value 


t value 




0.47835 


7.89 


a-oPo Pi5i 


0.66294 


12.83 


O-oPo Globe 


-0.33283 


-5.11 


O'oPo C^high 


-0.31872 


-4.87 



Table 6.6: Spearman Rank values for a^p, 
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Figure 6.11: (a - left) Source density against radio-luminosity and (b- right) against 
redshift. Symbols are as in Figure 6.5. 

However, as the Partial Spearman Rank results below show, the z-afpo correla- 
tion vanishes at any given radio-luminosity. This demonstrates that the apparent 
correlation of density with redshift is simply due to the strong Z-P151 correlation 
and the independent afpo-Pi5i correlation present in the samples. 

r ,3 p = 0.035 D = 0.51 

r p = 0.572 D = 9.38 

Tp = 0.524 D = 8.39 

The PCA results in Table 6.7 shows only a weak z-afpo correlation which reverses 
in the second eigenvector. This adds further evidence to the lack of any significant 
evolution of the source density with redshift. The Pisi-afpo correlation means 
that a more powerful source resides in a denser atmosphere. By comparing sources 
with a range of radio-luminosities, any bias with respect to density of the source 
environment in a sample is effectively removed. However, studies which only select 
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Parameter 


1 


2 


3 


4 


z 


0.515 


-0.213 


0.677 


-0.481 


Pl51 


0.587 


0.349 


0.168 


0.711 


Djobe 


-0.372 


0.813 


0.406 


-0.187 


a^Po 


0.501 


0.415 


-0.590 


-0.478 


Eigenvalue 


47.2% 


22.5% 


19.6% 


10.7 % 



Table 6.7: Eigenvectors and Eigenvalues for a^po 



radio sources with a narrow range of radio luminosities (e.g. Pedelty et al., 1989) 
may be affected by the Pi5i-afpo correlation. 

Hardcastle & Worrall (2000), Prestage & Peacock (1988) and Longair & Seldner 
(1979) found that low-z (hence low radio-luminosity) FRII sources do not live in 
rich environments. Yates, Miller & Peacock (1989); Crawford & Fabian (1989) and 
also Yee & Green (1987) find that the most powerful radio-galaxies inhabit rich 
environments, generally rich clusters. This is consistent with my results presented 
here. However, it is worth noting that all of the observational studies suffer from the 
redshift-radio luminosity degeneracy. Thus it is hard to determine whether the den- 
sity of the source environment changes with redshift or radio-luminosity. A recent 
study by Blundell et al. (1999) found that when they modelled the radio-luminosity, 
size, redshift and spectral index plane of the 3CRR, 6C and 7C sample their models 
did not require any evolution of the source density with redshift to match the obser- 
vations. This may indicate that any evolution of density with redshift is relatively 
unimportant. 

Table 6.6 also shows that there is a weak anti-correlation of density with a^igh. 
The Partial Spearman Rank results between ahigh, afpo and piobe below, demon- 
strate that the apparent a^Po-OiUgh correlation vanishes at any given lobe pressure 
and is simply due to the a^Po-piobe correlation coupled with the piobe-cxhigh anti- 
correlation. This proves that the spectral index of a source is insensitive to direct 
changes in the environment. However, since density and lobe pressure are corre- 
lated, any change in the source density is reflected in a change in the lobe pressure 
and hence atigh- 

r „a „ = -0.128 D = 1.86 

O^high 0,0 Po, Plobe 

r = -0.177 D = 2.58 

O^high Plobe, «o Po 

r /3 = 0.638 D = 10.92 

Plobe loPo, CXhigh, 



6.4.7 Jet power, Q 



Figure 6.12(a) shows that jet power, Qo, and radio-luminosity are highly correlated, 
with a significance exceeding 99.99% (see Table 6.8). This is expected from the 
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model as it predicts that the more luminous sources have a stronger jet driving 
them. However, the KDA model does not predict any relationship between jet- 
power and redshift. Table 6.8 shows that there is a strong correlation of jet-power 



10^ 



o 



ooo ,2*>*< 



*l*x ^ > 

ell' 



, XOO V 
Oq + * 



( X XX 



151 (MHz) 



10^ 



X ^x 
X tf+i* o 



50 , 



, o o o 



0.1 



10 



Figure 6.12: (a - left) Jet power against radio-luminosity and (b - right) against 
redshift. Symbols are as in Figure 6.5. 



Parameters 


IS value 


t value 


Qo PZobe 


0.64330 


12.18 


Qo F)zo6e 


-0.14354 


-2.10 


Qo z 


0.70195 


14.28 


Qo Pl51 


0.92486 


35.24 


Qo O^high 


-0.30324 


-4.61 



Table 6.8: Spearman Rank values for the jet power 

with redshift. This strong correlation is evident in Figure 6.12(b). The PCA results 
in Table 6.9 are dominated by the Q0-P151 correlation which only reverses in the 
final eigenvector. Conversely the correlation with redshift reverses in the second 
eigenvector and shows a much weaker correlation in comparison to the Q0-P151 
correlation. There is no significant trend with the size of the lobe. 



Parameter 


1 


2 


3 


4 


z 


0.501 


0.156 


0.829 


0.912 


P151 


0.578 


-0.227 


-0.455 


0.638 




-0.222 


-0.927 


0.291 


0.079 


Qo 


0.604 


-0.254 


-0.146 


-0.741 


Eigenvalue 


58.3% 


24.3 % 


14.0% 


3.4% 



Table 6.9: Eigenvectors and Eigenvalues for jet power 



By using the Partial Spearman Rank test, it is possible to disentangle the rela- 
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tionship between redshift and radio-luminosity and jet-power. 



P^5, = 0.223 D = 3.26 



151 

^2 Pi51, Qo 

^Pi5iQo,.= 0.854 r> = 18.32 



0.162 D = 2.35 



It is obvious that although the z-Qo correlation is significantly weaker than the 
P151-Q0 correlation it does not completely vanish at any given radio-luminosity. 
FRII radio sources are powered by jets created in the vicinity of the most massive 
black holes known in the universe. It is commonly assumed that the mass of the 
black hole determines the jet power. Lacy et al. (2001) find that for their sample 
of sources taken in the redshift range 0.3 < z < 0.5 the black hole mass - 5 GHz 
radio luminosity relation can be described by L5GHZ oc M^^^^^'"^ whereas Dunlop 
et al. (2003) find Lsghz oc M^^ between 0.1 < 2; < 0.25. A recent study by McLure 
et al. (2004), at a narrow range of redshift {z ^ 0.5), find that to Sa there is a trend 
between the low frequency luminosity, Lisimhz, and black-hole mass. The fact that 
I also find a correlation between Qo and z suggests that the mass of the black-hole 
powering the radio sources changes with redshift as well as with the low-frequency 
radio-luminosity. This is evidence that the mass of a black-hole powering FRII 
sources is not constant with redshift. The difference in the findings of Lacy et al. 
(2001) and Dunlop et al. (2003) may be attributed to the slight difference in the 
cosmological epoch of their samples. 

Table 6.8 shows that there is a weak ahig^-Qo anti-correlation but, as the Par- 
tial Spearman Rank results show, this is simply due to the independent Qo-Piobe 
correlation and piobe-cxhigh anti-correlation. 

rpiobe Qo, = 0.603 D = 10.04 

W««..,Qo= -0.198 D = 2.88 

^a.,..Qo,P..e= -0-118 D^1.71 



6.5 Conclusions 



I find no trend of the polarisation parameters DM^ and ctrm^ with the density of the 
lobe environment, a^po- This provides further evidence that the polarisation trends 
do not map changes in the density of the source environment, but are indicators of 
the degree of disorder in the structure of the magnetic field in these environments 
evolving with redshift. 

It was suggested in chapter 4 that many of the sources from the small sample 
show evidence of beaming. However, the density estimate from the KDA model 
disagreed with this finding, providing evidence instead that it is changes in the lobe 
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environment that causes the observed changes in the minimum energy magnetic field. 
The lobe pressure and the minimum energy magnetic field are tightly correlated, 
since the pressure from the KDA models is directly related to estimates of the 
minimum energy magnetic field. 

I find little evidence of any spectral ageing between 151 MHz and 1400 MHz. 
There is no trend of the lobe pressure with the low frequency spectral index, but a 
weak trend with the higher frequency spectral index which indicates, as expected, 
that radiative losses become more important at higher observing frequencies. 

Unsurprisingly I find a strong correlation with lobe pressure and size. The lobe 
pressure of a source also shows a strong anti-correlation with radio-luminosity and 
a weaker correlation with redshift. The lobe pressure of a source was also found to 
correlate with the density of the source environment. However, density is found to 
show a strong trend with radio-luminosity but there is no corresponding trend with 
redshift. This indicates that although density and pressure are related, it is only 
the lobe pressure that is significantly affected by changes in redshift. There is no 
indication of any evolution of the source environment with cosmological epoch. 

Finally, I find a very strong correlation between jet power and radio-luminosity 
with a much weaker correlation with redshift. This provides evidence of the evolution 
of black hole masses with low-frequency luminosity, but more importantly with 
redshift. 



Chapter 7 



Conclusions 



I have presented a complete data set of three samples of radio galaxies and radio- 
loud quasars. The three samples were defined such that two of them overlap in 
redshift and two have similar radio luminosities, allowing a study into the trends 
of various properties with redshift and radio luminosity. The spectral index, the 
rotation measure and the depolarisation measure were derived directly from the 
radio observations. 

The spectral index was found to be insensitive to the fundamental parameters 
(redshift, radio-luminosity and radio-size) and thus 1 can not confirm previous find- 
ings of a correlation of spectral index with redshift (Athreya & Kapahi, 1999), a 
trend of spectral index with radio luminosity (Veron et al., 1972; Onuora, 1989) nor 
a trend of spectral index with size (Blundell et al., 1999). This is most probably 
due to the comparatively small size of my sample. 

All sources were found to have an external Faraday medium that was local to 
each source and is responsible for variations of RM on small angular scales. The ob- 
served rotation measure itself was insensitive to changes in any of the fundamental 
parameters which is consistent with a Galactic origin. This Galactic contribution 
dominates the RM properties of the sources on large scales. At small angular scales, 
e.g. (JRM^ (rms variation of rotation measure) and DM^ (depolarisation measure), 
the Galactic contribution is effectively removed. The higher redshift sources dis- 
play a greater variation in their rotation measures properties. Only the observed 
depolarisation shows any significant trend with radio-luminosity. My results are 
consistent with both Morris & Tabara (1973) and Kronbcrg et al. (1972), as I find 
both redshift and radio-luminosity to correlate with depolarisation. 

The lack of any asymmetries in grm^ and DM with the minimum energy magnetic 
field suggests that the environment has no direct impact on the depolarisation or 
CflM^ in a source, but it is the degree of disorder in the structure of the magnetic field 
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that is changing with redshift. This is consistent with my findings that there is no 
strong correlation between the projected sizes of the sources and their depolarisation 
measure which agrees with the findings by Dennett-Thorpe (1996), but disagrees 
with findings by Strom (1973); Strom & Jagers (1988); Pedelty et al. (1989) and 
Ishwara-Chandra et al. (1998). The lack of any trend between (Trm^ and DM^ and 
the density of the source environment (from the KDA models) adds further evidence 
against the assumption that depolarisation and density are related. 

Asymmetries in spectral index show evidence of beaming in the majority of my 
sources. The results from the minimum energy magnetic field asymmetries are 
inconclusive and could also be explained by an environmental effect. The density 
estimates from the KDA model suggest that it is more likely that environmental 
differences and not beaming that causes the observed asymmetries. 

In general all asymmetries were found to be insensitive to changes in the source 
redshift, or low frequency radio luminosity; suggesting that the environments of the 
high redshift and low redshift samples are similar and any asymmetries caused are 
due to localised changes in the environment. 

I did not find any significant spectral ageing between 151 MHz and 1400 MHz 

for the sources in the 3CRR, 6CE and 7C III samples. The lobe pressure shows 
only a weak trend with the spectral index between 365 MHz and 1400 MHz. This 
indicates, as expected, that radiative losses are more important at higher observing 
frequencies. 

The 3CRR, 6CE and 7C HI sources also show that the density of the source 
environment correlates with radio-luminosity but seemed to be insensitive to changes 
in redshift. This is consistent with the findings of Longair & Seldner (1979); Prestage 
& Peacock (1988); Yates et al. (1989) and Yee & Green (1987) that there is no 
indication of any evolution of the source environment with cosmological epoch. 

Finally, I find that jet power and radio-luminosity are tightly correlated. This 
is expected as the more luminous sources have larger black-holes and thus have 
stronger jets emerging from their accretion disk (Urry & Padovani, 1995; McLure 
et al., 2004; Wang et al., 2003; Lacy et al., 2001). However, I also find evidence of 
jet-power and hence black-hole mass evolution with cosmic epoch. 



7.1 Evolution of the Faraday screen 

The observational trends noted above are further complicated by the presence of an 
external Faraday screen. Using my observations of variations in the rotation measure 
and the Tribble (1991) models I was able to determine the effects and evolution of 
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the Faraday screen with redshift, radio-luminosity and source size. 

The observed depolarisation-redshift trend was found to be best explained in 
terms of the intrinsic rotation measure, urm , varying as a function of redshift and 
a constant cell size, Sq. This is consistent with my earlier findings that it is the 
increasing disorder in the magnetic field and not density, that changes with redshift. 
A slight variation in the cell size, typically no more than ±50% around the best fit 
value, describes the scatter in the depolarisation results at any given redshift. In 
fact, it is generally the high redshift results that constrain the model. This suggests 
that to map the evolution of the Faraday screen, for all epochs, we must use samples 
containing sources from a large range of redshifts. 

Although the Burn (1966) law over-estimates the degree of depolarisation, I find 
similar results using the Tribble (1991) models, which are more reahstic. The sim- 
plistic Burn model however, can still be used to find the basic trends with depolari- 
sation in a sample, with fewer assumptions about the underlying distribution. I find 
that linear polarisation is proportional to best describes the data, compared to 
the widely used Burn result p ~ A"''. 

This work indicates that there is no significant evolution of the density of the 

source environment with cosmic epoch. Almost all of the observed trends with 
redshift are directly related to the changes in the degree of disorder in the magnetic 
field with redshift. The weak trend noted between jet-power and redshift could be 
attributed to a difference in the formation of black-holes at higher redshifts, i.e. a 
higher rate of mergers at higher redshifts. 

The significant parameter in determining source characteristics seems to be the 
radio-luminosity of the source. This is evident in both the small and large samples. 
This agrees with the findings of many other authors in different wavebands (e.g 
Alexander & Leahy, 1987; Yates et al., 1989; Prestage & Peacock, 1988; Inskip 
et al., 2003). 

Thus in conclusion, I find that there is evidence for a relationship between radio- 
luminosity and the environment in which a given radio source lives, but there is no 
significant evolution of the source environment with redshift. 



7.2 Further Work 

As with any significant body of work there are many avenues in which to extend 
the research. Ideally 1 would like to make the same observations between 1.4 GHz 
and 4.8 GHz for all sources in the 3CRR, 6CE and 7C HI. This is highly unrealistic 
and would also be too time consuming. A compromise would be to extend my small 
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sample to lower-luminosities and higher redshifts. A sample of around 50-60 sources 
would be feasible and would represent a large cross section of the FRIZ population. 

It would also be interesting to do a multi-waveband study of the sources already 
in my small sample. A X-ray study of each source could determine if the gas halo 

around each source was contributing to the depolarisation and if there was any 
interaction between the radio lobes and the gas detected due to its X-ray emission. 

Several of the high redshift 7C sources show rather twisted morphologies and 
the signal-to-noise level was too low to get detailed depolarisation and rotation 
measure maps. Higher resolution and longer observations would show if there is any 
underlying physical reason for their distorted structure. 



Appendix A 
Sample A maps 

Please see http://www.astro.soton.ac.uk/jag/Thesis.zip for full images. 



Appendix B 
Sample B maps 

Please see http://www.astro.soton.ac.uk/jag/Thesis.zip for full images. 



Appendix C 
Sample C maps 

Please see http://www.astro.soton.ac.uk/jag/Thesis.zip for full images. 
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Cosmology 



D.l Basic assumptions used 

The frequency at which I observe my sources is redshift dependent, 

'-'emit = (1 + ^) '-'observed (-D-l) 

where z/emit is the emitted frequency, ^'observed is the observed frequency and z is the 
redshift of the source. Although I observe at 1.4 GHz and 4.8 GHz the emitted 
frequency will be slightly different for each source. This is the reason why the po- 
larisation observations are shifted to a common redshift, thus removing any residual 
redshift effects. 

To convert from the angular size of a source on the sky, ^, to a physical size, 
Dsource) I usc the angular diameter distance formula. 



_ RoSr 6 , . 

^source ^ \-*-^"^/ 



where, R.s, = (-) ^^^^^ . (D.3) 

where ^2 is the density parameter of the Universe, c is the speed of light and Ho is 
the Hubble constant. 

The flux density of a source is also dependent on the redshift of a source, 
where is the radio-luminosity at a specific frequency u. 
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3CRR sample 
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2.33" 


0.93 


PH74 


3C 184 


0.99 


5 


14.4 


9.08 


2.58^ 


2.06 


LP91 


3C 184.1 


0.12 


182 


14.2 


9.57^ 


3.30" 


2.84 


LP91 


3C 186 


1.06 


3 


15.4 


6.58 


1.24^ 


2.00 


AG95 


DA 240 


0.04 


2109 


23.2 


130.00^ 


0.19^ 


1.48 


T82 


3C 190 


1.20 


7 


16.4 


9.09 


2.73^ 


1.21 


AG95 


3C 191 


1.95 


5 


14.2 


7.49 


1.85^ 


2.62 


AG95 


3C 192 


0.06 


200 


23.0 


14.39" 


4.80" 


3.18 


BHB88 


3C 196 


0.87 


6 


74.3 


49.02 


15.01^ 


0.75 


ASZ91 


3C 200 


0.46 


25 


12.3 


6.41 


2.04^ 


1.90 


DRAGN 


4C 14.27 


0.39 


36 


11.2 


4.41 


1.03^ 


2.30 


G04 


3C 204 


1.11 


37 


11.4 


5.49 


1.38^ 


2.92 


HLL94 


3C 205 


1.53 


18 


13.7 


9.27 


2.26^ 


1.29 


BMS88 


3C 207 


0.69 


13 


14.8 


8.92 


2.61^ 


1.37 


PH74 


3C 208 


1.11 


14 


18.3 


10.33 


2.36" 


1.91 


HLL94 


3C 212 


1.05 


9 


16.5 


8.34 


2.37^ 


1.50 


ASZ91 


3C 215 


0.41 


56 


12.4 


6.30^ 


1.59" 


1.52 


HLL94 


3C 216 


0.67 


30 


22.0 


13.93 


4.23^ 


1.70 


AG95 


3C 217 


0.90 


12 


12.3 


8.37 


2.16" 


2.88 


BLR95 


3C 219 


0.17 


190 


44.9 


14.00^ 


8.46" 


1.50 


NRH99 


3C 220.1 


0.61 


35 


17.2 


9.23^ 


2.24^ 


2.29 


JPR77 


3C 220.3 


0.69 


10 


17.1 


11.59 


2.89" 


1.19 


JPR77 


3C 223 


0.14 


306 


16.0 


8.51" 


3.58" 


2.95 


LP91 


3C 225B 


0.58 


5 


23.2 


10.17 


3.34^ 


1.98 




3C 226 


0.82 


35 


16.4 


9.34 


2.39" 


1.80 


BLR97 


4C 73.08 


0.06 


1004 


15.6 


24.20'^ 


2.65^ 


1.33 


DRAGN 


3C 228 


0.55 


46 


23.8 


10.44 


3.69^ 


2.40 


JLG95 


3C 234 


0.18 


112 


34.2 


15.78 


5.38" 


3.29 


HPPR97 


3C 236 


0.10 


2478 


15.7 


7.60 


3.28^ 


1.13 


ALB94 


3C 239 


1.78 


11 


14.4 


7.37 


1.46" 


1.38 


BLR94 


4C 74.16 


0.81 


40 


12.8 


25.90-^ 


2.30-^ 


1.98 




3C 241 


1.62 


1 


12.6 


7.18 


1.69" 


0.80 


PPR85 


3C 244.1 


0.43 


51 


22.1 


13.89 


3.94" 


4.32 


F02 
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Source 


z 


e 


Sl78MHz 


SseSMHz 


S 1400MHz 




R 








(arc sec.) 


Jy 


Jy 


Jy 




REF 


3C 


245 


1.03 


9 


15.7 


9.45 


3.31^ 


1.98 


- 


3C 


247 


0.74 


15 


11.6 


7.61 


2.88^ 


4.15 


ALB94 


3C 


249.1 


0.31 


47 


11.7 


7.55-^ 


2.34^ 


1.59 


HLL94 


3C 


252 


1.11 


60 


12.0 


6.56 


1.21^ 


1.12 


G04 


3C 


254 


0.73 


15 


21.7 


12.57 


3.13^ 


1.98 


- 


3C 


263 


0.66 


51 


16.6 


11.12 


3.00^ 


1.98 


- 


3C 


263.1 


0.82 


7 


19.8 


16.96 


3.13^ 


0.98 


NRH99 


3C 


265 


0.81 


78 


21.3 


12.18 


2.84^ 


2.95 


G04 


3C 


266 


1.28 


5 


12.1 


6.91 


1.43^ 


1.98 


AG95 


3C 


267 


1.14 


38 


15.9 


9.07 


2.16^ 


1.77 


G04 


3C 


268.1 


0.97 


46 


23.3 


10.02^ 


6.80^ 


1.87 


G04 


3C 


268.3 


0.37 


1 


11.7 


11.27 


3.72^ 


1.56 


AG95 


3C 


268.4 


1.40 


10 


11.2 


6.59 


1.98^ 


1.98 


- 


3C 


270.1 


1.52 


12 


14.8 


9.74 


2.85 ^ 


1.69 


ALB94 


3C 


274.1 


0.42 


150 


18.0 


5.12^ 


2.78^ 


3.56 


AL87 


3C 


275.1 


0.56 


16 


19.9 


9.31^ 


2.90^ 


1.29 


ALB94 


3C 


277.2 


0.77 


55 


13.1 


7.09 


1.85^ 


1.67 


PRL97 


3C 


280 


1.00 


15 


25.8 


14.94 


5.10^ 


0.98 


G04 


3C 


284 


0.24 


178 


12.3 


7.07 


1.93^ 


3.14 


HPPR97 


3C 


285 


0.08 


180 


12.3 


3.14^ 


2.04^ 


1.25 


AL87 


3C 


289 


0.97 


10 


13.1 


8.10 


2.40^ 


2.34 


BLR97 


3C 


292 


0.71 


140 


11.0 


5.34 


2.07^ 


3.82 


AL87 


3C 


294 


1.78 


15 


11.2 


6.12 


1.32 ^ 


1.98 


- 


3C 


295 


0.46 


6 


91.0 


65.44^ 


22.82^ 


1.73 


AL87 


3C 


299 


0.37 


11 


12.9 


8.17 


3.15^ 


1.17 


G04 


3C 


300 


0.27 


101 


19.5 


10.47 


3.70^ 


3.23 


HPPR98 


3C 


303 


0.14 


47 


12.2 


7.24 


2.67^ 


1.14 


HPPR98 


3C 


318 


1.57 


1 


13.4 


9.21 


2.69^ 


1.25 


AG95 


3C 


319 


0.19 


105 


16.7 


7.66 


2.50^ 


2.49 


HPPR97 


3C 


321 


0.10 


307 


14.7 


9.31^ 


3.50^ 


2.04 


BHB88 


3C 


322 


1.68 


33 


11.0 


6.94 


1.98^ 


1.53 


LLA95 


3C 


324 


1.21 


10 


17.2 


8.91 


2.44 ^ 


1.07 


G04 


3C 


326 


0.09 


1206 


22.2 


6.60^ 


3.10^ 


3.23 


DRAGN 


3C 


325 


0.86 


16 


17.0 


12.12 


3.56^ 


3.70 


FBP97 


3C 


330 


0.55 


60 


30.3 


21.88 


7.10^ 


1.20 


F02 


3C 


334 


0.56 


58 


11.9 


6.80 


2.15^ 


1.58 


HLL94 


3C 


336 


0.93 


28 


12.5 


7.88 


2.74^ 


1.65 


HLL94 


3C 


341 


0.45 


74 


11.8 


6.35 


2.15^ 


3.14 


HLL94 


3C 


340 


0.78 


46 


11.0 


7.79 


2.49^ 


2.10 


JPR97 


3C 


337 


0.63 


46 


12.9 


9.98 


2.94^ 


2.07 


BLR97 


3C 


343.1 


0.75 


0.3 


12.5 


8.00^ 


4.83^ 


0.89 


AG95 
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Source 


z 


e 

(arc sec.) 


Sl78MHz 

Jy 


SseSMHz 

Jy 


S 1400MHz 

Jy 






R 

REF 


3C 349 


0.21 


86 


14.5 


9.95 


3.30^' 


2 


.77 


HPPR97 


3C 351 


0.37 


74 


14.9 


8.37 


3.27^ 


1 


.34 


G04 


3C 352 


0.81 


13 


12.3 


7.85 


1.87^ 


1 


.57 


ALB94 


3C 356 


1.08 


75 


12.3 


6.39 


1.41^ 


3 


.00 


FBB93 


4C 16.49 


1.30 


16 


11.4 


6.27 


1.46 ^ 


1 


.10 


G04 


4C 13.66 


1.45 


6 


12.3 


6.95 


1.70^ 


1 


.49 


LGS98 


3C 368 


1.13 


8 


15.0 


7.40 


1.11^^ 


2 


.67 


BLR95 


3C 381 


0.16 


73 


18.1 


8.71 


3.910^ 


2 


.11 


HPPR97 


3C 382 


0.06 


185 


21.7 


6.55 


5.83^ 





.86 


F02 


3C 388 


0.09 


50 


26.8 


16.04 


5.74^ 


1, 


.21 


PH74 


3C 390.3 


0.06 


229 


51.8 


30.19^ 


11.90^ 


1 


.40 


F02 


3C 401 


0.20 


24 


22.8 


14.34 


5.07^ 


1 


.49 


HPPR97 


3C 427.1 


0.57 


27 


29.0 


15.52^ 


3.80^ 


2 


.78 


NRH95 


3C 432 


1.81 


15 


12.0 


6.48 


1.58^ 


1 


.33 


BHL94 


3C 433 


0.10 


68 


61.3 


31.27 


11.80^ 


1 


.98 




3C 436 


0.21 


109 


19.4 


5.74^ 


3.40^ 


2 


.75 


HPPR97 


3C 437 


1.48 


34 


15.9 


8.68 


2.76^ 


1 


.94 


BLR97 


3C 438 


0.29 


23 


48.7 


26.40 


6.37^ 


1 


.81 


HPPR97 


3C 441 


0.71 


37 


13.7 


8.34 


2.52^ 


1 


.58 


FBP97 


3C 452 


0.08 


280 


59.3 


32.67^ 


10.54^ 


2 


.16 


DRAGN 


3C 455 


0.54 


4 


14.0 


9.18 


2.89^ 


1 


.46 


AG95 


3C 457 


0.43 


210 


14.3 


5.68 


1.75^ 


2 


.45 


G04 


3C 469.1 


1.34 


74 


12.1 


6.46^ 


1.76^ 


1 


.98 




3C 470 


1.65 


24 


11.0 


5.98 


1.94^ 


3 


.35 


BLR97 



Table E.l: Parameters used in the fitting on the KDA 
model in chapter 6. Column 1: 3CRR source name. Col- 
umn 2: redshifts. Column 3: Angular size in arcseconds. 
Column 4: Flux density at 178 MHz. Column 5: Flux 
density at 365 MHz taken from Douglas et al. (1996) 
unless otherwise noted. Column 6: Flux density at 1400 
MHz and Column 7: Average source lobe ratio. Columns 
2, 3 and 4 use data taken from Willot (2003). Sources 
with no high quality maps have lobe ratios which are 
the average for the 3CRR sample. Sec Table E.2 for an 
explanation of the reference codes used. 
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Listing Full Reference 



A 


\ A / In X r \-< o/^lyoT" 

VVlllL(:; OC JI>6CK61 


^^iyyz J 


13 


ijdige et di. 


l^iyoi ) 


P 


J-ld/lIlg 06 red/COCK 


(-1 QQ.r\\ 


u 


WPN<^<^ POT T ARDR ATTDN 


^ 1 yyo J 




OOllQOll et dl. 


^^iyyo ) 


r 


iveiiei iiidiiii ei di. 


(A QP.Q\ 

^^iyoy j 


P 


xvueiii et dl. 


(^ QS1 
l^iyoi j 


tr 

1 1 


r dUllIiy- ±(jLi1 cL dl. 


\^iyoo j 


T 
i 


V lliei o6 H/I ICKbOIl 


^^iy ( j 


T 
J 


ndieb ei dl. 


^^iyyo j 


K 

IV 


ijecKei et di. 


^^iyyo ) 


T PQ1 


ijediiy o£ It ei ley 


1 iyy i 1 




VjUUtlieL ei dl. 


yZ,\j\j'± j 




J-idW vji cell cL dl. 






INLii dl. 


^ 1 yyo j 


r uz 


r ei mill 




A T RQ/I 


/VKUJOI et dl. 


('[ QQA\ 

1 iyy^t 1 


r>R A P 1\T 


Lieany ex ai. 


[ZUUU ) 


IJPPRQ7 


Hardcastle et al. 


M oo7^ 
l^iyy / j 


WPRQS 
rlr^riyo 


Hardcastle et al. 


(^lyyo ) 


nijijy^ 


1— cviiH \ OT o 1 

J311Qlt; cL dl. 


1 lyy^ ) 




AKUJOI et dl. 


n QQ1 

1 iyy i 1 


BLR97 


Roof of q1 


(1997) 


LCFOl 


Lara et al. 


(2001) 


PPR85 


Pearson et al. 


(1985) 


FBB93 


Fernini et al. 


(1993) 


AG95 


Akujor & Garrington 


(1995) 


SC85 


Strom & Conway 


(1985) 


PH74 


Pooley & Henbest 


(1974) 


T82 


Tsien 


(1982) 


BHB88 


Baum et al. 


(1988) 


JPR77 


Jenkins et al. 


(1977) 


BMS88 


Barthel et al. 


(1988) 


JLG95 


Johnson et al. 


(1995) 


AL87 


Alexander & Leahy 


(1987) 


FBP97 


Fernini et al. 


(1997) 


LGS98 


Ludke et al. 


(1998) 



Table E.2: Reference table 



Appendix F 
6CE sample 



Source z 9 Si5imhz Ssssmhz Si4oomhz R 









( PI TP *?PP 1 


Tv 




Tv 


TTl Tv 




REF 






i.oO 


94 


z.oy 


i , 


1 

. iO 


91 "^-4 

Z id 


i .Oo 




6C 


0822+341 


0.41 


18 


3.06 


1, 


.77 


552-^ 


1.83 




6C 


0822+343 


0.77 


21 


2.93 


1, 


.03 


111^ 


1.64 




6C 


0823+375 


0.21 


81 


3.35 


1, 


.42 


425^ 


1.41 




6C 


0824+353 


2.25 


8 


2.42 


1, 


.83 


958 ^ 


1.80 


NAR92 


6C 


0825+345 


1.47 


7 


2.10 


1, 


.22 


284^ 


1.95 




6C 


0847+375 


0.41 


33 


3.07 


1, 


.43 


614^ 


2.67 


MMC97 


6C 


0857+390 


0.23 


24 


2.71 


1, 


.51 


503^ 


1.83 




6C 


0901+355 


1.90 


4 


2.07 


1, 


.08 


219^ 


1.24 


NAR92 


6C 


0902+341 


3.40 


5 


2.14 


1, 


.12 


298^ 


2.42 


COM94 


6C 


0905+395 


1.88 


5 


2.82 


0, 


.94^ 


233^ 


2.85 


LEL95 


6C 


0908+373 


0.11 


39 


2.33 


1, 


.33 


614^ 


1.08 


PRF86 


6C 


0913+390 


1.25 


9 


2.27 


1, 


.80 


1005-^ 


0.80 


HBW95 


6C 


0919+380 


1.65 


10 


2.72 


1, 


.05^ 


65^ 


1.84 


NAR92 


6C 


0922+364 


0.11 


17 


3.27 


1, 


.09 


725^ 


1.85 


MPC97 


6C 


0930+385 


2.40 


5 


2.21 


1 


.06 


278^ 


1.84 


FPDOl 


6C 


0943+395 


1.04 


12 


2.31 


1, 


,20 


323^ 


2.01 


G04 


6C 


0955+384 


1.41 


22 


3.45 


1, 


.75 


388^ 


1.66 




6C 


1011+363 


1.04 


66 


2.10 


1, 


.07 


374^ 


2.07 




6C 


1016+363 


1.89 


31 


2.28 


1, 


.39 


528 ^ 


1.61 




6C 


1017+371 


1.05 


9 


2.68 


1 


.44 


359^ 


1.83 




6C 


1018+372 


0.81 


83 


2.52 


0, 


.97 


271^ 


2.07 




6C 


1019+392 


0.92 


9 


2.99 


L 


.76 


392^ 


2.05 


BEL99 


6C 


1025+390 


0.36 


1 


2.97 


1, 


.98 


659^ 


2. 6 


NAR92 


6C 


1031+340 


1.83 


3 


2.33 


1, 


.55 


478^ 


2.72 




6C 


1042+391 


1.77 


11 


2.68 


1, 


.96 


655-^ 


1.57 


NAR92 


6C 


1043+371 


0.79 


5 


2.62 


1, 


.18 


249^ 


1.07 


NAR92 
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Source 


z 


9 

(arc sec.) 


SlSlMHz 

Jy 


SseSMHz 

Jy 


S 1400MHz 

mJy 






R 

REF 


6C 


1045+340 


1.83 


22 


2.00 


1.20 


266^ 


1 


.98 




6C 


1045+355 


0.85 


9 


2.07 


1.26 


267 ^ 


1 


.44 


NAR92 


6C 


1045+351 


1.60 


0.1 


3.03 


2.44 


1035-^ 


1 


.83 


- 


6C 


1100+350 


1.44 


14 


2.26 


1.18 


304^ 


4 


.25 




6C 


1108+395 


0.59 


16 


2.10 


2.01 


390^ 


1 


.83 


- 


6C 


1113+345 


2.41 


17 


2.33 


1.39 


440^ 


2 


.80 




6C 


1123+340 


1.25 


0.2 


3.40 


3.80 


1378^ 





.98 


NAR92 


6C 


1125+374 


1.23 


18 


2.07 


1.06 


325^ 


1 


.83 


- 


6C 


1129+371 


1.06 


19 


2.36 


1.42 


416^ 


1 


.72 




6C 


1130+345 


0.51 


78 


3.20 


2.10^ 


524^ 


2, 


.31 




6C 


1134+365 


2.13 


17 


2.07 


1.10 


235^ 


1 


.27 




6C 


1141+352 


1.78 


12 


2.40 


1.13 


290^ 


1 


.65 




6C 


1143+370 


1.96 


0.1 


2.06 


1.51 




1 


.83 


- 


6C 


1148+363 


0.14 


27 


3.21 


1.41 


550^ 


1 


.83 


- 


6C 


1148+384 


1.30 


10 


3.83 


2.20 


652^ 


2, 


.22 


NAR92 


6C 


1158+343 


0.53 


40 


2.12 


1.30 


314^ 


1 


.83 


- 


6C 


1159+365 


1.40 


2 


2.20 


1.40 


352^ 


1 


.28 




6C 


1204+351 


1.38 


63 


3.43 


1.96 


505^ 


1 


.49 




6C 


1205+391 


0.24 


24 


3.83 


1.61 


634^ 


1 


.55 


MMC82 


6C 


1212+380 


0.95 


0.6 


2.14 


1.27 


287^ 


1 


.20 


NAR92 


6C 


1213+350 


0.86 


0.1 


2.39 


2.32 


1507^ 


1 


.25 


XRP95 


6C 


1217+364 


1.09 


0.5 


2.40 


1.32 


396^ 





.93 




6C 


1220+372 


0.49 


36 


2.52 


1.44 


463^ 


1 


.83 




6C 


1230+345 


1.53 


12 


2.90 


1.89 


424^ 


2 


.30 




6C 


1232+394 


3.22 


51 


3.27 


1.47 


1130^ 


1 


.70 


NAR92 


6C 


1255+370 


0.71 


0.6 


3.66 


2.33 


759^ 


1 


.83 




6C 


1256+364 


1.13 


18 


2.88 


1.61 


527^ 


1 


.53 




6C 


1257+363 


1.00 


40 


2.40 


1.07 


196^ 


1 


.80 




6C 


1301+381 


0.47 


28 


3.46 


2.01 


536^ 


1 


.83 





Table F.l: Parameters used in the fitting on the KDA 
model in chapter 6. Column 1: 6CE source name. Col- 
umn 2: redshifts. Column 3: Angular size in arcseconds. 
Column 4: Flux density at 151 MHz. Column 5: Flux 
density at 365 MHz taken from Douglas et al. (1996) un- 
less otherwise stated. Column 6: Flux density at 1400 
MHz and Column 7: Average source lobe ratio, high res- 
olution maps taken from Law-Green et al. (1995) unless 
otherwise stated. Columns 2, 3 and 4 use data taken from 
Rawlings et al. (2001). See Table F.2 for explanation of 
reference codes. 



Listing Full Reference 



A 


White & Becker 


(1992) 


E 


Condon et al. 


(1998) 


I 


Viner & Erickson 


(1975) 


J 


Hales et al. 


(1995) 


K 


Becker et al. 


(1995) 


L 


Ficarra et al. 


(1985) 


M 


Pilkington & Scott 


(1965) 


N 


Becker et al. 


(1991) 


G04 


Goodlet et al. 


(2004) 


HBW95 


Henstock et al. 


(1995) 


PRF86 


Parma et al. 


(1986) 


COH94 


Carilli et al. 


(1994) 


LEL95 


Law-Green et al. 


(1995) 


MPC97 


Morganti et al. 


(1997) 


FPDOl 


Fanti et al. 


(2001) 


BEL99 


Best et al. 


(1999) 


NAR92 


Naundorf et al. 


(1992) 


MMC82 


Machalski et al. 


(1982) 


XRP95 


Xu et al. 


(1995) 



Table F.2: 6CE Reference table 
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Table G.l: Parameters used in the fitting on the KDA 
model in chapter 6. Column 1: 7C III source name. 
Column 2: redshifts taken from Lacy (1999). Column 3: 
Angular size in arcseconds taken from Lacy et al. (1992). 
Column 4: Flux density of 151 MHz. Column 5: Flux 
density of 327 MHz. Column 6: Flux density of 1490 
MHz and Column 7: Average source lobe ratio. Sources 
with a lobe ratio of 1.8 do not have high enough quality 
maps. All other lobe ratios use maps taken from Lacy 
et al. (1992). Columns 4-6 are taken from Lacy (1999). 
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